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Cover photo: Hurricane Hugo, which snapped trees and dumped salt water on 
South Carolina forests in September 1989, suddenly and dramatically altered a 
pime ecosystem. W.T. Swank photo was taken in the vicinity of the North Inlet 
LT ER Site. 
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INTRODUCTION TO LTER WORKSHOP 


ON CLIMATE VARIABILITY AND ECOSYSTEM RESPONSE 


David Greenland and Lloyd W. Swift Jr.? 


The Intersite Climate Committee of the Long- 
Term Ecological Research (LTER) Program, which is 
sponsored by the National Science Foundation, has 
the mission of facilitating investigations of the 
etmospheric environment in LTER ecosystems. The 
Committee has developed standards for 
meteorological measurements at LTER sites 
(Greenland 1986; Swift and Ragsdale 1985) and has 
summarized the climates at the first 11 LTER sites 
(Greenland 1987). This climate summary 
demonstrated the obvious: very different 
ecosystems have very different climates. This 
report, and discussions at the LTER Data 
Processing Workshop at Las Cruces in January 1986, 
suggested that moisture (including soil moisture) 
and climate variability were distinctive forcing 
variables at each site. The Climate Committee 
decided to defer investigations of water budgets 
and to concentrate first on climate variability. 
This decision recognized the LTER network’s 
potential importance to research on the global 
climate change question and the growing public 
interest in that question. 


Eleven of the 15 sites then in the LTER network 
attended a workshop on Climate Variability and 
Ecosystem Response. Ten sites are represented by 
papers in this volume. Each site was invited to 
examine its longest time series of climatic data 
for temporal variability and to comment on the 
relation of that variability to ecosystem 
responses. The variability of many data sets was 
characterized by multiyear climatic periods 
punctuated by strong and dramatic responses to 
specific weather events. All sites found duration 
of record and spatial representativeness to be 
limiting factors in assessing variability. 


The workshop was held in August 1988 at the 
Mountain Research Station of the University of 
Colorado, the field headquarters for the Niwot 
Ridge-Green Lakes Valley LTER site. The keynote 
address on Global Warming and Ecosytem Response 
was given by Dr. Stephen Schneider of the National 
Center for Atmospheric Research. Following this 
and formal presentations of papers, the authors 
and Dr. Gary Cunninghem from the Jornada LTER site 
met as the LTER Climate Committee to review the 


1 Professor of Geography, University of 
Colorado, Boulder, CO; and Research 
Meteorologist, Coweeta Hydrologic Laboratory, 
Forest Service, U.S. Department of Agriculture, 
Otto, NC. 


research and our understanding of climate 
variability and ecosystem responses of the LTER 
sites. The discussions reported in the overview 
chapter concluded that recognition and utilization 
of time and space scales are keys to understanding 
response phenomena. 


The papers in this volume represent a variety 
of ecosystems, environments, and approaches to our 
topic. This variety represents one of the riches 
of the LTER program, although it makes 
standardization and generalization difficult. The 
volume starts with information from three forest 
sites. Federer examines the record at the Hubbard 
Brook Experimental Forest in New Hampshire in 
order to determine whether real variation in the 
climate can be deduced from existing records. A 
powerful statistical technique, the Z-T extreme 
event analysis, is employed by Swift and others to 
determine the uniqueness, or return period, of 
extreme events in streamflow and precipitation 
data from the Coweeta Hydrologic Laboratory in 
North Carolina. Viereck and Adams infer effects 
of climate warming on vegetation patterns from 
data on spatial variation in microclimate and 
related plant successional development at the 
Bonanza Creek Experimental Forest. 


Aquatic ecosystems were represented by the 
North Inlet South Carolina and the Northern Lakes 
Wisconsin LTER sites. Michener and others show 
impacts of chronic and acute climate events upon 
the estuary ecosystem and how the scale of 
climatic variability affects productivity. This 
paper was written before the North Inlet site was 
severely impacted by Hurricane Hugo in September 
19898. At the Wisconsin site, Robertson 
demonstrates how historical data for fresh-water 
lakes can be used as a measure of longer term 
climatic change. Predominantly agricultural 
landscapes were addressed by Wendland and by Crum. 
Wendland summarizes the history and quality of 
Illinois weather observations that supported the 
former Illinois Rivers LTER site. Crum found 
little or no indication of climate change in 100- 
year temperature and precipitation records and was 
able to relate corn yield to midsummer 
precipitation at the Kellogg Michigan LTER site. 


Three sites represented landscapes where 
extreme climates limit vegetation cover. 
Greenland describes a marked variation in the 
temperature and precipitation record on the alpine 
tundra at the Niwot Ridge-Green Lakes Valley 
Colorado site, but this variation was not well 
correlated to obvious ecosystem responses. In 
contrast, Kittel reports that climate variability 





in the shortgrass steppe ecosystem at the Central 
Plains Experimental Range site in Colorado is well 
related to ecosystem function. The final site 
paper, by Hayden, describes how storm events at 
the Virginia Barrier Island site move and reform 
coastal terrain and influence vegetation 
distribution. The latter point was a good 
demonstration of the contrast between time and 
space scales, which must be recognized in any 
attempt to relate climate variability to ecosystem 
response. 
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CHANGE, PERSISTENCE, AND ERROR 


IN THIRTY YEARS 


OF HYDROMETEOROLOGICAL DATA 


AT HUBBARD BROOK— 


1/ 


1 2 
C. Anthony Federer— 





Abstract.--Daily precipitation, air temperature, and solar 
radiation data have been collected on the Hubbard Brook 
Experimental Forest, New Hampshire, for over 30 years. A 
tradeoff occurs between cost and accuracy. Various 
instrument errors can make real climatic variation difficult 
to detect. Periods of above or below "normal" temperature 
or precipitation persist for up to several years, but their 
ecosystem effect is probably slight. Rare hurricanes cause 
the greatest ecological response to any weather events. 


Keywords: Temperature, precipitation, solar radiation. 


INTRODUCTION 


The Hubbard Brook Experimental Forest (HBEF) 
in New Hampshire is a recent addition to the 
Long-term Ecosystem Research (LTER) network, 
but its weather records and its research history 
extend for more than 30 years. Only in the 
last couple of years has the Hubbard Brook weather 
record been organized on a computer in such a 
way that long-term analysis is relatively 
simple. 


In this paper I will examine possible 
long-term change or variation in the Hubbard 
Brook record, to see whether real variation can 
be separated from errors. The possible effects 
of climatic variation on the ecosystem are 
briefly discussed. 


The differences among weather, climatic 
variation, and climate change are differences of 
time scale. The climate of a place is usually 
defined in terms of means over a 30-year period. 
Climate change, therefore, takes place over 
decades; weather, on the other hand, over days 
to weeks. This leaves open the question of 
whether a drought or cold period of several 
to many months should be considered as weather 
or climate. Though I would prefer such phenomena 


1/ F . . : 
—/Presented at the LTER Climate Committee 


Workshop on Climate Variability and Ecosystem 
Response, Nederland, CO, August 22-23, 1988. 
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to be classed as weather variation, the term 
“climatic variation" is now commonly used for 
such phenomena that have time scales of months 
to several years. 


Variation or Error? 


Most of us know, when we read an instrument, 
that the reading may be slightly incorrect. With 
weather instruments and data analysis, there are 
several sources of error, some of which can be 
confused with climatic variation and even 
climate change. 


Mistakes are human-induced by misreading 
or miscopying. If they are large enough, they 
may be caught by eye or by a computer program, 
but smaller ones will always exist and will be 
merged in with the randomness inherent in 
weather data. Modern electronic data collection 
systems may be able to eliminate mistakes. The 
effect of random mistakes on analysis of climatic 
variation and change is probably negligible. 
However, persistent mistakes that bias data 
could be interpreted as real variation. 


Other sources of error are more important 
when evaluating climate variation. Sensor 
calibrations change, sensors are replaced, 
instrument exposure changes or instruments are 
moved, and processing procedures change. Each 
of these error sources causes apparent, but not 
real, long-term variation in the data. I will 
point out several possible instances in the 
Hubbard Brook data set. 


There is little question about how to 
minimize such errors in weather data. With 
enough equipment and technicians (i.e. 
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money), highly accurate and complete data could 
be produced. But at ecological research sites, 
our main objective is not data collection but 
research. We must decide on the tradeoff between 
accuracy and cost. We must decide what sort of 
error rate or amount is tolerable. At Hubbard 
Brook over 30 years, the demand for or use of 
weather data, aside from watershed precipitation, 
has been minimal, though not quite non-existent. 
We have put a great deal of effort into 
collecting some data that have never been used. 
Concern for error reduction is wasted on data 
that will not be used. 


HUBBARD BROOK HISTORY 


The Hubbard Brook Experimental Forest was 
established by the USDA Forest Service 
(Northeastern Forest Experiment Station) in 
1955. HBEF is a single, oval-shaped basin of 
3,160 ha located in the White Mountain National 
Forest of central New Hampshire (Fig. 1). Basin 
elevation ranges from 222 to 1,015 m. HBEF 
is drained from west to east by Hubbard Brook, 
which has many similar-sized tributaries on 
opposing south-facing and north-facing slopes. 
The area was selected specifically for small 
watershed research; stream gages were built on 
eight of these tributaries between 1956 and 1967. 
The gaged watersheds range in area from 12 te 
76 ha. Three of them have been treated by 
cutting the 70-year-old northern hardwood 
(beech-birch-maple) forest that covers virtually 
all of HBEF. 
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Instrumentation 


Collection of precipitation and air 
temperature data began immediately after 
establishment of HBEF. Openings were cut in 
the forest to eliminate tree crowns above a 45 
angle from the precipitation collectors. The 
network density of standard rain gages is about 
1 per 20 ha; 11 are on adjacent south-facing 
watersheds 1-6, and 9 are on the adjacent 
north-facing watersheds 7 and 8. Weighing 
recording gages are located at five of these 
locations, and hygrothermographs in standard 
weather shelters are at three of them (plus one 
at Headquarters). All gages are mounted high 
enough to clear up to 2 m of snow, and are 
equipped with Alter-type windshields. Each 
weather shelter also includes maximum, minimum, 
and standard thermometers. 


oO 


Electric line power and telephone 
communications have never been available anywhere 
at HBEF except at Headquarters. AI11 
precipitation, streamflow, and temperature 
recorders are spring-wound or battery operated. 
Spring-wound analog recorders have proven very 
reliable for year-round operation in the ~-30 to 
+40°C environment. Visitors to Hubbard Brook 
are usually surprised by our "primitive" 
approach to hydrometeorological data collection. 
The fundamental reasons for this are a philosophy 
that prefers to spend available money on 
research rather than on data collection, and a 
belief that inexpensive mechanical systems are 
more reliable than expensive electronic ones. 





~~ Hubbard Brook 
* rain gage 
+ weather stn. 
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Figure 1.--Map of the Hubbard Brook Experimental Forest 
showing rain gage and weather station locations. 
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Data Processing 


Data processing, on the other hand, has 
changed drastically over time. In the 50's and 
early 60's, all data were picked from charts by 
eye and tabulated manually using rotary 
calculaters. In 1964, we began digitizing 
streamflow charts and processing data by computer. 
Over the years, we have used six different mainframe 
and mini computers (GE-265, IBM 1620, IBM 360, 
DEC-10, PRIME~750, and DG-Mv4000).3/ Our Fortran 
programs migrated rather easily. Input and 
output data were stored for many years on punch 
ecards. Now, most of our streamflow and 
meteorological data are stored on magnetic 
tapes in ASCII files; this method has been 
chosen because it moves easily to new machines. 


Precipitation and temperature data continue 
to be read from charts by eye, but are 
processed by computer. Only daily total 
precipitation and daily maximum and minimum 
temperatures are obtained routinely. 


In early years, when the amount of data was 
less and labor was cheap, every data point and 
every calculation was checked by a second 
person. Now computer programs catch only 
the grossest errors. There is certainly a 
possibility that our random error rate now 
is higher than it used to be. 


Headquarters Weather Station 


At Hubbard Brook, we emphasize what happens 
on the gaged watersheds. The foot of the 
lowest and closest of these is 3 km away and 
200 m higher than our Headquarters (HQ) office. 
Though a weather station has been maintained 
at Headquarters, data from it have not always 
been processed. For instance, we have miles of 
unprocessed charts from a tipping bucket rain 
gage, and from an anemometer and wind vane. 
Headquarters data provide neither complete nor 
representative values for weather on the gaged 
watershed areas. This raises a question of 
which weather station or stations to use as a 
climatic standard for HBEF. 


A weather station meeting LTER Level 2 
standards was established at HQ (where power 
is available) in 1981. A variety of problems 





3/ 


~The use of trade, firm, or corporation 
names in this publication is for the 
information and convenience of the reader. 
Such use does not constitute an official 
endorsement or approval by the U.S. Department 
of Agriculture or the Forest Service of any 
product or service to the exclusion of others 
that may be suitable. 


have beset this system, and there are many gaps 
in the record. For example, a week was lost 

once when I forgot to push the ON button. 

Sensors, tape recorder, and data logger all fail 
occasionally. Lightning damage has been frequent, 
though we keep trying to improve protection. 
Repair and maintenance of this station has not 
had top priority. Although the daily data that 
it produces have sometimes been used, there have 
never been any requests for its hourly data. 


Data Amount and Availability 


Data piles up over the years. About one and 
a half million values are now contained in the 
routine hydrometeorological data set at 
Hubbard Brook (Table 1). And each of these 
has a time and location attached to it. Yet, 
we still get requests to "send us all your 
precipitation and streamflow data!" 


Our philosophy has been to provide this data 
to anyone who requests it. We have recently 
established the Hubbard Brook Bulletin Board. 
Anyone with an MS-DOS microcomputer and a 
modem (300, 1200, or 2400 baud) can obtain 
many of the files containing this data by 
calling 603-868-1006 outside of normal working 
hours. Routine or long-term data from many 
cooperating scientists in the Hubbard Brook 
Ecosystem Study will gradually be added to the 
Bulletin Board. Data are provided to users in 
exactly the same ASCII files that are produced 
by the scientists. 


SOLAR RADIATION 


Daily total solar radiation has been measured 
at Headquarters since 1960. Sengors have 
included successively a Belfort— pyranograph, 
a Weather Measure pyranograph, and two LiCor 
pyranometers. The pyranographs were calibrated 
occasionally against a Kipp pyranometer. 
Averaging the pyranograph charts by eye over 
2-hour intervals provides good daily totals. 
For the LiCor sensors, the manufacturer's 
calibration has been used and hourly integrals 
are obtained by the data logger. Comparisons 
among the various sensors generally have shown 
differences of less than 5Z. 


The scatter of data by day of the year shows 
two problems (Fig. 2). A few data points 
exceed the potential insolation but are not 
obvious errors and have not been eliminated. 
Second, the midsummer maximum daily values are 
only 80-85% of potential, though 90Z would be 
expected for clear days. The site horizon is 
elevated up to 15° in the northeast and northwest 
by nearby trees. These trees have been growing 
gradually taller, and presumably reducing daily 
solar radiation gradually, especially in summer. 
Should we cut down the trees or move the sensor? 
Or would that destroy any usefulness the 





Table 1.--Approximate amount of raw hydrometeorological data at HBEF 





Number of Days per 
Years Year 
Precipitation 30 x 365 
Temperature 30 x 365 
Streamflow 30 x 365 
HQ Weather Station 7 x 365 


mM POM 





Values per Number of Number of 





Day Stations Values 

1 x 22 = 241,000 

2 x 5 = 110,000 

8 x 8o/ = 701,000 

24 x 6— 368,000 
1,420,000 








ane values per hour 
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Figure 2.--Scatter plot of daily solar radiation. 
at HBEF versus day of the year, for 1960-87. 
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potential insolation for a flat horizon. 


long-term record might have for detecting 
real change in radiation? 


Missing values were estimated through the 
1960's using National Weather Service data from 
several stations in New England. That network 
no longer exists. Since 1981, the pyranograph 
provides a backup instrument. To obtain 
annual totals, values for a few missing days 
in each year had to be guessed. How much 
effort should be put into estimating such 
missing values? 


The annual radiation over 1960-87 has no 
long-term trend (Fig. 3), but several deviations 
may be instrument related. Are the 3 low 
years, 1967, 1968, 1969, real, or was the 
Belfort instrument calibration changing before 
its replacement? Are the low 1986, 1987 values 
due to the second LiCor sensor? Without 
considerably more effort in calibration and 
cross-checking of sensors, it is doubtful that 
such a record could detect long-term change, 
such as that caused by increasing haze or 
cloudiness. 


ATR TEMPERATURE 


Air temperature data begin in 1957 for HQ and 
Station 1, in 1961 for Station 6, and 1965 


for Station 14. Which station or stations 
should be used to look at long-term trends? 
Averaging over several stations reduces the 
effect of any errors and the effect of 
geography in any one station but shortens the 
record. The number of stations included in an 
average cannot change over time or a bias will 
be introduced. 


From an ecological viewpoint, daily maximum 
and minimum values contain more information 
than the daily mean. The average diurnal range 
and the normal extreme temperatures can be 
important to injury and survival of plants and 
animals. Similarly, seasonal differences are 
also important. The variation over years of 
average minimum and maximum temperatures for 
winter (December, January, February) and summer 
(June, July, August) provides a reasonably 
comprehensive picture in a small amount of space 
(Fig. 4). There is no evidence for gradual 
change and the year to year variations appear 
random. Only the sharp increase from winter 
1981-82 to 1982-83, and from summer 1982 to summer 
1983 suggests some abnormality such as an 
instrument problem. Because year-to-year 
variation is high, and errors are possible, 
evidence for climate change should not be 
looked for in records for a single location, 
even over a 30-year period. 
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Figure 3.--Annual solar radiation (3 jm”) at HBEF for 
1960-87. Arrows show times of sensor changes. 
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Figure 4.--Average summer (June, July, August) and winter 
(December, January, February) maximum and minimum daily 
temperatures for Stations 1, 6, and HQ at HBEF, 1961-87. 


stations, abrupt changes do not occur, but 
gradual shifts of stations with respect to 

The deviations of individual stations from each other do occur. For instance, winter 
the averages described above may be partly temperatures at Station 6 were about L°C 
due to instrument calibration (Fig. 5). The higher than Station 14 in the early record, 
hygrothermographs are currently readjusted and but were about the same after the mid 1970's. 
rotated annually in spring. In earlier times, Watershed 4, in which Station 6 is located, 
an instrument stayed in the same location for was strip clearcut in 1970-74. It is 
several to many years. Standard and max-min tempting to attribute some change to cutting, 
thermometers are read weekly at each weather but such change is not obvious. We have 
station, but the hygrothermograph data are not strongly discouraged interstation comparison 
corrected to agree closely with the thermometers. of temperatures, because of this lack of 
In general, we have expected that errors of accuracy. Yet the general order of temperatures 
+1°C were likely. The analysis for this paper from highest to lowest--HQ, Station 1, Station 
of year-to~year variation within a station 6, Station 14-~is in the expected inverse 
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shows a range of about 2°C (Fig. 5). The sudden relation to elevation. HQ and Station 14 differ 

sharp drop in Station 1 summer maximum from by about 3°C in maximum temperature, and by 

1970 to 1971 is suspicious, but is not evident 470 m in elevation, just what is expected from : 

in the summer minimum. The 1969 winter data for the lapse rate. The data may be better than I 

Station 1 are also suspect. For other years or we think. | 
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Figure 5.--Deviation by station of summer and winter maximum 
and minimum temperatures from the average values in Figure 


4. 


Persistence 


Ecologists customarily use monthly or 
annual averages to evaluate abnormality in 
weather, but this ignores or discards the 
important information contained in the daily 
values. If July was 1°C above normal, does 
that mean that each day was 1°C above normal 
or that a few days were many degrees above 
normal? If a summer is characterized as "hot", 
what does that really mean in terms of 
distribution of temperature over time? 


In work with tree rings in the Northeast, 
I first began looking at accumulated deviation 
of daily mean temperature from its normal. I 
was struck by the apparent persistence for 
months of below or above normal periods, and 
by the abrupt transitions from one "regime" 
to another. Station 1 at Hubbard Brook shows 
such behavior over its 31 years (Fig. 6). In 
such a plot, also used by Barry (1985), positive 
slope indicates above normal conditions and 
negative slope below normal. 


On a year-to-year time scale, below 
normal temperature persisted from mid-1964 to 
early 1968, and from early 1980 to late 1982. 


Conditions were generally above normal from 
early 1968 to early 1976, and from late 

1982 to early 1985. A comparison with 
regional weather data used in our tree-ring 
studies shows that these persistence regimes 
are regional in scope. 


On a month-to~month time scale, regimes 
occur with a persistence of several months 
(Fig. 7). November 1969 through April 1970 
was almost uniformly below normal. The use 
of daily data shows that transition between 
these regimes can be very abrupt, and can 
often be pinpointed to a particular day. The 
intensity of the persistence regimes is not 
great. Regimes of about 1°C above or below 
normal are most common, with few extended 
periods reaching more than 2°c difference 
from normal. 


Time series analysis of this data yields 
nothing more than a lag 1-day autocorrelation 
of about 0.5. The series is like taking two 
steps in one direction and sliding back one 
before deciding randomly on the next direction 
and step size. Persistence of the type shown 
here is characteristic of a random walk, or lag 
1 autocorrelation. Nonetheless, it is tempting 
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Figure 6.--Accumulated deviation of daily mean temperature 
for Station 1 from its normal for that day of the year, for 


1957-87. 
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Figure 7.-~Accumulated deviation of daily mean temperature 
for Station 1 from its normal for that day of the year, 


for 1967-71. 


to believe that the transition from one regime 
to another is forced by regional or even 
global weather patterns, such as jet stream 
movement (Blackmon et al. 1977), that may be 
produced by El Nino Southern Oscillation 
episodes, or North Pacific sea surface 
temperature anomalies (Namias et al. 1988), 
which have similar persistence. 


PRECIPITATION 
Persistence 


Persistence behavior occurs for 


precipitation as well as for temperature (Fig. 8). 


The "normal" precipitation for a given day of 
the year needs to be smoothed because the 
actual values have large day-to-day variation. 
We used a cubic smoothing spline with a cutoff 
of 50 days. The accumulated deviation tends 

to go upward rapidly and in big jumps, due to 
individual large storms, and come downward 

more slowly. Because of the skewed distribution 
of daily precipitation, normal time series 
modeling cannot be done. 


The 3l-year record for Watershed 1 shows 
a single below-normal period from late 1960 
through mid-1966 (the well-known northeastern 
drought), and above-normal conditions from early 
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Figure 8.--Accumulated deviation of daily precipitation for 
Watershed 1 from its smoothed normal for that day of the 


year, for 1957-87. 


1972 through mid-1978 (Fig. 8). The magnitude 
of deviation in the persistent regimes is on 
the order of 30 mm from normal per month, or 
about 307%. 


Station Differences 


Daily data from individual standard 
rain gages have not been entered onto the 
computer before 1965, so longer-term 
comparisons can only be made on a watershed 
basis. Both double mass comparison of one 
watershed with another and plotting of 
deviations from means show various, mostly 
unexplained, blips and kinks in the data. I 
have looked for one specific effect here. 
The cutting of Watershed 4 in 1970-74 altered 
the three rain gages in it so that they were 
in the open rather than in openings. Watershed 
4 precipitation does seem slightly reduced 
in the period 1972-77, and then drifts 
upward (Fig. 9). The reduction is only 
about 2% and may be meaningless. There is 
a slight tendency for Watersheds 1 and 3 
to drift downward while 6 and, recently, 4, 
drift upward. There is no obvious reason for 
such drift to exist. The data are dependent 
only on standard rain-gage catch, so no 
sensor or calibration error is involved. 
Exposure change by tree growth around the 
opening is possible; the openings are 
occasionally enlarged to compensate for 
ingrowth at the edges. Vegetation within the 
openings is cut only every several years. 
Should we put any more effort into examining 
these possible drifts and correcting for 
them? 
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ECOSYSTEM EFFECTS 
Weather Variation 


Weather variation has many effects on a 
forest ecosystem. These effects are often 
understood qualitatively, but are poorly 
quantified. We believe that drought reduces 
tree growth, but cannot specify how much. We 
know temperature affects photosynthesis and 
respiration, but cannot relate growing season 
temperature to net primary productivity. 
Recently I tried to relate ring widths of red 
spruce in New England to physiologically-based 
weather variables, with no success (Federer et 
al. 1988). Red spruce seems to be affected 
more by injurious winter conditions than by 
summer conditions (Johnson et al. 1986). But 
quantifying the weather conditions that cause 
winter conifer injury has proven elusive. 
Similarly, weather is known to affect outbreaks 
of pathogens, but quantification and modeling 
of the interactions are in their infancy. In 
general, we know that the variations in weather 
discussed in this paper affect the ecosystem, 
but we cannot say by how much. 


Most effects of weather and climate 
variation on ecosystem processes are complex 
and non-linear. But many analyses of the 
relationships assume simple linear additive 
responses to monthly precipitation and mean 
temperature. The question of whether the 
available data is appropriate to the research 
problem is often avoided. For instance, 
precipitation data, or even Palmer drought 
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Figure 9.--Ratio of annual precipitation on each watershed 
to the mean of watersheds 1, 3, 4, and 6, for 1964~-86. 


index, are not good substitutes for 
soil-water deficits in terms of drought effects. 


With respect to known Hubbard Brook 
history, the single weather event with the 
most impact on the ecosystem was probably the 
New England Hurricane of September 21, 1938. 
Enough trees were uprooted to induce the USDA 
Forest Service to carry out salvage logging in 
the Hubbard Brook basin. This was the only 
severe weather~related disturbance since 
the area was heavily cut prior to 1920. 
Hurricanes and local wind storms cause both 
regional and local disturbance at irregular 
intervals. However, treethrow mounds are 
common at Hubbard Brook; treethrow provides 
both soil mixing and seedbed for certain 
species. 


Fire is another weather-related event 
with potentially severe ecosystem impact. Fires 
that burn surface litter are fairly common in 
northeastern forests in early spring. These 
have some effect on nutrient cycling. Only 
one such fire, of a few ha, has occurred in 
the last 30 years at Hubbard Brook. More 
severe crown fires are very rare, with 
prehistoric return periods of hundreds of years, 
but ecosystem impacts obviously would be 
severe. 


Climatic Variation and Change 


Climatic variation, as evidenced by 
persistence in the temperature and 
precipitation records, occurs essentially 
continuously at Hubbard Brook. Periods of 
months to years occur with temperature 1°¢ 
above or below normal, or precipitation 302% 
above or below normal. Such persistence may 


favor certain species of plants or animals 
over competing species, but demonstration and 
quantification of such cause and effect would 
be quite difficult. 


The 30-year weather record at Hubbard Brook 
is too short to detect climatic change if 
climate is defined by a 30-year mean. Climatic 
change is very difficult to detect in any event; 
elaborate examination of many studies and data 
sets has not been able to prove the existence of 
a CO,-induced warming trend (Ellsaesser et al. 
1986). Looking for climate change in data 
sets from LTER sites is futile. Nevertheless, 
climate change has and will affect Hubbard 
Brook. Hubbard Brook is at the temperate 
forest-boreal forest ecotone. Climatic 
temperature changes of a couple of degrees C 
may greatly change the ratio of spruce-fir to 
northern hardwoods and, thereby, the 
structure and function of the whole ecosystem. 
Changes in the ecosystem itself may become 
evident long before the weather data can 
prove that the climate has changed. 


Chemical Climate 


Although excluded from this paper, changes 
in the chemistry of the atmosphere may have 
greater effects on the ecosystem than changes 
in physical weather and climate. Acid 
precipitation, ozone and SO, levels, and 
heavy metal deposition can be considered part 
of weather and climate in their larger sense. 
Acid precipitation and related air pollution 
have certainly increased considerably in 
the northeastern United States over the past 
100 years. More recently at Hubbard Brook, 
sulfate input in precipitation has decreased as 
regional sulfate emissions are reduced (Likens 
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et al. 1984). Lead deposition to the forest 
floor also has decreased recently. The 
ecosystem impacts of these pollutants therefore, 
have been partially reversed. The chemical 
climate at Hubbard Brook probably is changing. 
The impacts via acidification of soils and 
water and consequent ecosystem effects are 
receiving considerable attention at Hubbard 
Brook and elsewhere. 


CONCLUSIONS 


Climatic data collection and processing has 
changed almost incredibly over the past 30 
years. We cannot predict what the next 30 
years will bring, but need to expect more 
great changes. 


In spite of the electronic age, mechanical 
weather sensors may still be more reliable 
than hi-tech equipment unless funds are made 
available for a full-time electronic technician 
and data analyst, backup equipment, and careful 
calibration. 


For temperature and solar radiation, if 
not for precipitation, weather variation, 
persistence, and climatic change may be 
difficult to separate from error. Replication 
of sensors and systems can help. 


Where long-term records already exist, it 
may be better to continue with existing 
instruments, methods, and exposure than to make 
changes that will alter long-term mean values. 


There is a tradeoff between accuracy of 
weather data and costs. Where ecological 
research is the top priority, weather data 
should not be expected to detect climatic 
change. 


Ecological systems may be affected more by 
occasional extreme events than by 
long-term changes in mean weather or climate. 


Our knowledge of the quantitative relations 
between weather and ecosystem processes lags 
far behind our ability to collect reasonably 
good weather data. We should be spending 
much more money and time on research and 
maybe much less on instrumentation. 
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APPLICATION OF THE Z-T EXTREME EVENT ANALYSIS 


USING COWEETA STREAMFLOW AND PRECIPITATION DATA! 


Lloyd W. Swift, Jr., Jack B. Waide, and David L. White? 


Abstract.-—— A technique for drought or flood analysis, 
after Zelenhasi¢ and Todorovié, promises to improve the 
definition of both duration and magnitude of extreme flow 


events for river-sized basins. 


The technique has been 


applied to a smaller basin at the Coweeta LTER Site using 
both streamflow data and a longer precipitation record. 
This report illustrates the technique and describes needed 
adjustments to apply the method to stream-sized basins. 


Keywords: LTER, drought, duration, recurrence interval. 


INTRODUCTION 


Long-Term Ecological Research (LTER) sites are 
mandated to sponsor research in five core areas. 
One core area deals with natural and human-caused 
disturbances and their impacts on ecosystems. 

For most sites, common natural disturbances are 
driven either by short-term meteorological events 
such as storms and droughts or by long-term 
climate episodes. In either case, the researcher 
wants to be able to state how an apparently 
unique event or trend fits with past experience. 
In hydrology, traditional analyses of extreme 
value data are based upon flow magnitude and 
often use only the most extreme annual event 
(Chow 1964, Section 81). These analyses yield 
little information about duration of the event 
and ignore other unusal but less extreme events 
falling within the same year. Duration is 
incorporated in the Palmer Drought Index, but the 
ability to deal with probability of return 
interval is not. Zelenhasié and Salvai (1987), 
who extend work by Todorovié and Zelenhasié 
(1970), illustrate a method that considers the 
entire data record and allows statements to be 
made about the recurrence intervals of both the 
magnitude and duration of an extreme event. 


1 Presented at the LTER Workshop on Climate 
Variability and Ecosystem Response, August 21-24, 
1988 at Niwot Ridge-Green Lakes LTER Site, CO. 


2 Research Meteorologist, Coweeta Hydrologic 
Laboratory, Otto, NC., Principal Research 
Ecologist, Forest Hydrology Laboratory, Oxford, 
MS, and Ecologist, Southeastern Forest Experiment 
Station, Clemson, SC, all Forest Service, U.S. 
Department of Agriculture. 


Their illustration is based on a 8 800 000-ha 
river basin. We propose that the method can be 
applied to smaller basins if certain refinements 
in technique are made which allow for greater 
responsiveness of streams compared to rivers. 


Our application of the Z—-T method seeks to 
describe the southeastern drought of 1985-1988 and 
its significance to the Southern Appalachian 
Mountains and the Coweeta LTER site. Questions 
appropriate to the drought disturbance are: how 
dry was this period? how long did it last? and how 
does it rank with other droughts on record? In 
order to apply the Z-T method, we had to determine 
whether the technique could function with data 
from a smaller (approximately 1/10,000th) basin 
and whether the technique could be used with 
precipitation data. If the latter were feasible, 
the period of record could be doubled to over 100 
years and comparison of the 1925 and 1986 droughts 
would be possible. 


CONCEPT OF THE Z-T METHOD 


For droughts, the Z-T method analyzes deficit 
events created by dividing a continuous streamflow 
record into periods of unusually low flow 
alternating with periods of all higher flow rates. 
The separation is made around a derived flow 
reference value, Q,. A sequence of deficit and 
inter—deficit flow periods may contain intervals 
of normal flow that are separated by short, small 
deficits. Likewise, an extended low-flow sequence 
may contain periods when flow rates are 
temporarily above @,. These minor excursions 
above or below Q, are culled from the analysis by 
averaging them into the adjacent flow periods. 
Statistical tests for validity and serial 
correlation are applied to the resulting set of 
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largest deficits and curves are fitted to the 
cummulative relative frequencies of the deficits 
and their durations. From these fits, recurrence 
intervals are calculated for the most extreme 
events. 


THE Z-T PROCEDURE 


Mean daily streamflow values for Coweeta 
Watershed 8 (WSO8) from January 1936 through 
December 1988 provided the data for the 
streamflow test of the method. Watershed 8 is a 
forested headwater basin of 760 ha draining into 
the Little Tennessee River. Mean annual 
precipitation of 1988 mm is evenly distributed 
through the year, and flow averages 1163 mm per 
year. 


Selecting Q, 


All daily flow values for Coweeta WSO08 were 
ranked in descending order without regard to 
date; Q, is the value occurring at the "r" 
percent interval. For this work, Q, was selected 
as the flow rate ranked just under 90 percent of 
all larger daily streamflow values. Monthly Q,, 
values were also selected by ranking daily flows 
separately by months. Q,, may be appropriate for 
slowly changing data but responsive or flashy 
streams may require a Q,, or even larger reference 
value. The periods of consecutive days when flow 
was below Q, are the deficit events (Figure 1). 
The required data for the following analyses are 


TRUNCATION OF A FLOW SERIES 
INTG DISCRETE DEFICIT EVENTS 


STREAMFLOW 


Cumulative Deficit 
Flow Amount 






Reference 
ow 






1 t TIME. 


Figure 1.--Separation of a period of streamflow by 
Q, into flow deficit events. 


the durations of these deficits (T) and the 
period sums of the differences between Q, and 
recorded flow (D). 


The seasonal cycle of streamflow suggests that 
low flow is relative; 2 level that is unusually 
low in the spring wet season might be a moderate 
level for the fall season. Using a single Q,, for 
the whole year confines most deficits to the low- 
flow season (see days 270-310 in Figure 2). 
However, an atypically low flow, such as occurred 
in the spring of 1985, could have significant 
impact upon stream ecology. To include such 
cases, a separate Q,, was calculated and applied 
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to each month’s flow data; thereafter, deficits 
were identified in all seasons. Because a smooth 
transition between months was not possible when 
adjacent monthly Q,, were quite different, a curve 
was fitted to the monthly Q, and reference values 
calculated for each day of the year. Figure 3 
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Figure 2.—-Daily streamflow for Coweeta WSO8 for 
1985 showing reference values for annual and 
monthly Q,5.- 


compares the distribution of deficit periods in 
1986 derived by the three types of Q, and suggests 
that the extra effort needed to apply daily 
varying Q, was not justified by these data. 


Culling Minor Events 


The chronologic series of daily streamflow was 
transformed by subtraction into excess flows above 
Q, and deficit flows below Q,. Each unbroken 
series of deficit days was summed, defining a 
deficit event with duration T and cummulative 
deficit (deviation from Q,) D. Similarly, excess 
days were summed. Low-flow periods may be 
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Figure 3.—-Average flow for deficit and inter— 
deficit periods in 1986 as separated by annual, 
monthly, and daily varying Q,. 


punctuated by short events of flow above the 
reference value, usually caused by storms that 
briefly raise streamflow but have no lasting 
effect upon accumulating drought conditions. The 
next step was to remove these brief excess 
periods and to combine adjacent deficit periods 
that effectively function es a single event. The 
rules for combining were: (1) the magnitude of 
the excess had to be less than the absolute value 
of each of the adjacent deficits, and (2) the 
absolute sum of the two deficits had to be 
greater than three times the magnitude of the 
excess. The duration of the combined deficit 
event was the sum of the three durations and the 
cumulative deficit was the algebraic sum of the 
excess and two deficits. Iterative passes were 
made through the series of alternating deficit 
and excess events until no further minor excess 
events met the criteria for combining. A drought 
period extending through the end of the year was 
listed as falling in the year containing more 
than half of the duration. In the event of an 
even split, then the drought period was listed in 
the year having more than half of the deficit. 


We believe that the same iterative process 
also should be used to eliminate minor deficit 
events. The function fit to the cumulative 
relative frequency, to be described later, was 
sensitive to the number of very small events. 
The goal was to obtain a list of separate and 
independent events drawn from the lowest 10 
percent of flows on record. From this, lists of 
drought duration events and drought magnitude 
(cumulative deficit below the Q, value) events 
became the data points for the analytical steps. 


Statistical Tests 


Several tests were applied to ensure that the 
selected deficit events were independent, 
identically distributed random variables. First, 
a chi-square test determined whether the 
distribution of the number of droughts in each 
year differed from the distribution of a time~ 
dependent Poisson process. Then the list of 
deficits and the list of durations were each 
tested to ensure that the events were not 
serially correlated and that the lists did not 
include runs of consistently increasing or 
decreasing values. Finally, the drought deficits 
were ranked, the drought durations were ranked, 
and the correlation between rank numbers for each 
event was tested to confirm that both measures of 
drought severity ranked the seme events in 
essentially the same order. 


The number of droughts each year should fit 
the Poisson distribution. For this test, we 
calculated the overall mean annual number of 
deficits and counted the number of years having 
0,1,2,...k deficits each. For each class, k, the 
expected frequency was calculated as 


Feyp = Y * e(-Ny) * N,* / k!, where 


Y = total number of years, 
N, = mean annual number of deficits, and 
k = class interval of 0,1,2,... deficits/year. 


exp 


The sum of F,,, * k must equal the total number 
of deficits. If the sum of (Fy, ~ Fexp)?/Fexp 


across all k classes was less than the chi-square 
value for the number of classes minus 2, the list 
of selected deficits was judged to not differ from 
the Poisson distribution. 


Serial correlation was tested on the 
chronological list of deficit magnitudes and on 
the list of durations. Three correlations were 
calculated for each list: (1) each event with the 
event following, (2) each event with the second 
following, and (3) each event with the third 
following. For a serial correlation calculation 
of two offset sequences of the same data, we used 
the variance of all events in the denominator 
rather then the usual product-moment formula. If 
all six correlations were nonsignificant, then the 
lists of deficits and durations were judged to be 
acceptably free of serial correlation. If not, 
then we started over, revising our rules for 
culling minor events to effect more combining of 
adjacent drought periods. Alternatively, 
selecting a larger Q, may reduce serial 
correlation for some data sets. 


The runs test determined whether the 
chronological lists of deficits or durations 
contained trends or cycles. Each list was 
searched for groups of sequential values that 
formed runs of increasing or decreasing magnitude. 
Runs may consist of single values that reverse the 
previous and succeeding trends. Adjacent pairs of 
identical values were considered part of the same 
run. If r is the number of runs and n is the 
total number of values in the list, then the 
parameter for a t-test is approximated (Sokal and 
Rohlf 1969, pg. 628) by 


t, = {r — [(2n - 1)/3]} / {(16n - 29)/90}*?, 


where most lists would contain enough events to 
test against t,, = 1.96. If t, was greater than 
1.96 and the number of runs was low, the test 
suggested a systematic trend or bias in the list. | 
If the number of runs was high, cycles in the data 
were suspected. 


The two lists of deficits and durations were 
individually ranked and rank numbers assigned to 
the drought events in each list. Thus, any 
particular drought event might have different rank 
numbers for its deficit and its duration. The 
correlation between these pairs of rank numbers is 
a measure of the consistency of the data. 
Zelenhasié and Salvai (1987) found that reducing 
Q, increased the rank correlation. High 
correlation was not required for the frequency 
analysis but would be an important criteria if the 
results of this analysis were used to construct a 
synthetic drought series. 
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Recurrence Interval of Drought 


The end product of this analysis is the 
recurrence intervals of the largest drought events 
based upon the deficits selected and tested by the 
preceding steps. To calculate recurrence 
intervals we required: the lists of deficits (D) 
and durations (T), each ranked in ascending order; 
the total number of droughts (N); the mean annual 
number of droughts (N,); and the values of the 
mean deficit (D,) and mean duration (T,,). The 
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same calculation procedures, described below for 
the list of deficits, were applied also to the 
list of durations. 


The cumulative relative frequency (CRF) was 
calculated by dividing the rank number of each 
deficit by the total number of droughts+1; CRF 
therefore approached 1.00 for the very largest 
deficit. The expected distribution for CRF is 


CRF.,, = 1 - e(-D/D,); 


where D was each value from the list of deficits. 
The Kolmogorov-Smirnov goodness of fit test 
(Ostle 1963, pg. 471) was used to verify that the 
observed CRF was not significantly different from 
the expected distribution of CRF. For 
acceptance, the maximum difference between CRF,,, 
and CRF,,, had to be less than the K-S parameter 
of 1.36/(N)'” for the 0.05 level. Instead of 
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Figure 4.—-Observed points and expected curve for 
cumulative relative frequency of Coweeta WS08 
streamflow deficits. 


calculating all possible differences to find the 
maximum, one could select the most likely range 
of values from graphs drawn to visually verify 
the fit (Figure 4). If the distributions were 
not different, then the formula for the expected 
CRF could be used to represent the observed data. 


Because recurrence interval is usually stated 
in years, we calculated the distribution function 
of the largest annual deficit, which is an 
extension of the function for CRF,,,, as 


Fexp(D) = e[-Ny * e(-D/D,) ] 


where N, was the mean number of deficits per year 
and D was as in CRF,,,. The goodness of fit of 
this frequency distribution also could be 
verified by testing the maximum difference 
between the F,,, and a set of observed cumulative 
relative frequencies recalculated from an 
abbreviated list, omitting lesser deficits for 
years having more than one. The recurrence 
interval 


RI =1/ [1 - F,,,(D)] 


was calculated for the largest deficits. Also, 
the estimated deficit for a selected return 
interval was calculated as: 


Day = ~Dy * In[{ -In(1 -1/RI) / Ny]. 
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ILLUSTRATION OF THE Z-T METHOD 


Table 1 lists the 44 drought periods selected 
from Watershed 8 daily flow records for the 53 
years from 1936 through 1988. Droughts were not 
defined in 29 of the years. The distribution of 
droughts per year and two variations on the 
Poisson calculation are shown in Table 2. The 
first set of columns produce a chi-square sum much 
larger than the critical value of 9.49 for 6-2=4 
classes at 0.95 probability. Ostle (1963, 
pg. 127) recommends combining classes in cases 


Table 1.—-Coweeta Watershed 8 drought periods 
based upon a Q,, truncation applied to daily flows 
in area-~mm for 1936-1988. 


Year Julian Date Dur Def — Rank CRF Expected CRF — obs-exp 


Begin Mid End T 2D YD FF OD T 2D Pr. o8 


59 283 344 40 123 LAL 43 ONL 956.966.996.055 -.041 
40 4647 4B 53 15 1 & .022 .167 .079 .144 -.057 023 
40 4 HG 72 82 LT 867617 832.094 035 
a 8 8% 9 5 Le 3 13 056 .278 1128 .170 -.073 .108 
40 104 107 109 6 2.38 5 18 1M 400 2152 1212 -.041 188 
41 34108 181147 40.9 43 42.956 933.982 986 -.027 -.052 
41 284 292 299) 16 10 OD A833 355099078 -.065 
42 125) 129 134 DOL 3244067260 108-016 -.041 
43 Sda S52 SHB 132.2 1S 322 367.300 20412163 
4“ sO 6 26 AE 7B 373463105 1S 
7 8 43 6 3.2 2 22 433 489 355.282 «078.206 
4 177 216 25 078 BS 40 SL BBD 689 8B2 586 = 007.108 
S2. 204 2M 2i7 A RSNA .059 -. 013 
92 2% SO S23 OF 2.2 OF 7644 367.523.204.022 163 
98 278 (SOL 52345 GLH OOF 800 600.709 469091131 
4 255 293 Sl 76 10.7 39 34 B67 756 876.670 -.009  .085 
4 334 SSG 558 9 4 3 6 056 122 128 135 -.073 -.015 
4 355 «358 S61 7 2.0 7 14 156 311.175 187 -.019 .124 
S 8 3 2 0 8.9 24 32 533 711 422 603 11) 1108 
% Sl 4 8 8 13.3 3 FF 733 822 637.748.096.074 
% Soh 340 S46 137 AS WD 322.233.300.162 -.019 072 
3 MO 0 S60 4.2 6 25567556 438353129202 
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such as this where small counts and high chi- 
square values invalidate the test. The 
recalculation in the right side of Table 2, with 
3 or more deficits per year in the final class, 
yields a chi-square under the critical value of 
5.99. Thus, the distribution of droughts per 
year is judged to not differ from the Poisson 
distribution. 


Table 2.--Test for fit of number of droughts per 
year to the Poisson distribution. 


Class Expected Observed Obs-Exp Chi-Sq Observed Gbs-Exp © Chi-Sq 
k years years years 


Sum= 20.58 5.66 


Table 3 gives the results of the serial 
correlation tests on drought duration and 
deficit. All values of r are small. The runs 
test confirms the conclusion that neither the 
duration nor deficit lists are serially 
correlated. Each list contains 28 runs, thus the 
t-value by the approximate formula equals -0.37. 


Table 3.—-Serial correlation (r) to the third level 
for drought durations and deficits. 


Level N Duration Deficit 
ixitl 4 °°» - -.003t—«~ -.057 
ix i+2 42 —.097 ~.023 
ix i+3 41 206 . LO6 


The columns for duration in days (T) and 
deficit in mm (D) in Table 1 were each ranked in 
ascending order. The correlation between the two 
rank numbers for each drought was 0.84. At this 
point, the lists of selected droughts had met all 
the statistical tests and recurrence intervals 
could be determined. 


The cumulative relative frequencies (CRF) of 
durations and deficits were calculated by 
dividing rank numbers by 45 (Table 1). Based on 
the mean number of droughts, 0.83 per year, the 
expected CRF and the difference between observed 
and expected CRF were calculated. The largest 
deviation in durations was 0.129, well below the 
Kolmogorov-Smirnov criteria of 0.205 for 44 
observations. However, the largest deviation of 
deficits from the expected distribution was just 
barely above the K-S criteria. Figure 4 shows 
the result of this borderline fit in the mid- 
range of CRF for drought deficits. The penalty 
for a weak fit, as seen below, was an outrageous 
estimate of recurrence interval for the most 





extreme drought. Table 4 gives the estimated 
recurrence intervals for the eight most extreme 
droughts. Although we resisted accepting the 
recurrence intervals estimated for 1986, both the 
analyses and Figure 4 underscored the fact that 
277 consecutive days with streamflow in the lowest 
10 percent of all recorded daily flows was a 
highly unusual event. Figure 4 suggests that 
reducing the number of small droughts would pull 
the curve for observed data closer to the 
estimated distribution. A quick recalculation of 
selected CRFs, after we arbitrarily culled the 
four smallest droughts, decreased both the 
observed and expected CRF and their difference in 
midrange and cut the recurrence interval for the 
1986 deficit in half. Users of the Z-T method 
should seek to improve the fit of the CRF 
distribution by experimenting with several levels 
of combining small events or by trying other Q, 
values. 


Ten different years accounted for all the most 
extreme drought events in the Coweeta Watershed 8 
record. Table 4 lists the recurrence intervals 
for eight deficits and eight durations. Years 
1939, 1941, 1986, and 1988 ranked high in both 
lists. For the WSO8 data, a 100-year drought was 
estimated to last 161 days and have an accumulated 
deficit of 42.5 mm below the Q,, flow level. 


Table 4.—-Estimated recurrence intervals (RI) for 
the largest drought events from 1936 through 
1988. 


Year Dur RI Year Def RI 

days years 7 years 
86 277 2392 &6 92.0 16e13 
4} 147 68 39 54.1 330 
88 131 44 4] 40.9 84 
39 123 36 81 23.6 14 
5 78 il 8 21.7 12 
54 76 10 40 17.2 8 
85 70 9 él 16.3 7 
Gl 53 6 56 13.3 5 


PRECIPITATION COMPARISON 


Climate histories of the Southern Appalachian 
Mountains place the last major drought about 1925, 
a decade before the beginning of the climate 
record at Coweeta. A precipitation history began 
at Highlands, NC, in November 1877. This is a 
station in the U.S. Historical Climatology Network 
of the National Weather Service, 20.4 km east of 
Coweeta at the same elevation as the upper portion 
of Watershed 8. The Z-T method was applied to the 
lll-year Highlands record to determine whether 
precipitation could be used to compare the 
severity of the 1986 drought with that of 1925. 


Because precipitation is an event phenomenon, 
daily values could not be used as was done with 
continuous streamflow data. When monthly totals 
were analyzed, almost all durations were one month 
because individual storms would interrupt and mask 
otherwise dry periods. Streamflow is the 
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cumulative response of past precipitation input 
and evapotranspiration. As a test to simulate 
the timelag nature of streamflow, we used for 
each monthly value of precipitation the moving 
average of that month and the preceding three 
months. This technique identified 69 droughts 
with durations of 1 to 7 months. The 1925 
drought was identified by Highlands precipitation 
as the driest period (Table 5). Four of the 10 


Table 5.--Estimated recurrence intervals (RI) for 
the largest drought events identified by Highlands, 
NC, precipitation data, 1878-1989. 


Year Deficit RI 
mam years 
1925 168.1 334 
1954 154.7 21g 
1883 136.4 122 
1939 135.9 121 
1879 108.2 50 
194] 99.0 38 
1878 94.1 32 
1981 92.5 31 
1955 83.4 23 
1986 78.2 20 


most extreme events occurred before the Coweeta 
record began. Both the Highlands and Coweeta 
precipitation records identify the mid-1950’s as 
a more important drought period than did WS08 
streamflow. Precipitation data certainly does 
not rank the 1985-86 period as the most extreme. 


CONCLUSIONS 


The Z-T method of drought analysis seemed 
complex as presented by Zelenhasi¢é and Salvai 
(1987) with all the derivations. In application, 
the procedure is a series of logical steps and 
statistical tests to confirm that the data fit 
the assumed distributions. A purpose of our work 
was to determine whether the Z-T method could be 
applied to a fourth-order stream. We determined 
that the method could be used, but found that the 
responsiveness of our stream, compared to the 
large river, required careful selection of both 
the reference value Q, and the techniques for 
culling small events and combining them into 
adjacent extended wet or dry periods. An 
unsatisfactory fit to the cumulative relative 
frequency distribution was purposely included to 
illustrate the importance of obtaining a data 
series that fits the assumptions. 


A second purpose was to determine whether a 
longer precipitation record could be used to 
extend our interpretation of a 53-year streamflow 
record. Precipitation was not a perfect analog 
for streamflow, but this test showed that some 
form of time~averaging of precipitation can 
simulate the storage and release characteristics 
of the watershed. 


The third purpose of this application was to 
describe the 1986 drought at Coweeta and compare 
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it with previous events. This analysis shows 
conclusively that duration and deficit of 1986 low 
flows were unique in the 53-year streamflow record 
of Coweeta Watershed 8. The cumulative effects of 
drought extended into 1988. However, the 1939 
through 1941 period contained three separate 
drought periods that together must have had a 
severe impact on mountain watersheds and 
ecosystems. The extended precipitation record 
also identifies these same drought sequences but 
rates the 1925 drought as having the greatest 
rainfall deficit. Diaries and newspaper accounts 
report a dry period in the late 1800’s, and this 
is supported by the Highlands precipitation 
record. The low precipitation in 1954 does not 
appear as a major streamflow drought. Although 
the May—through-October period in 1954 established 
the all-time minimum growing season total for this 
station, average or above-average precipitation 
before and after these months apparently 
maintained soil moisture at levels sufficient to 
sustain streamflow. This may suggest that a 
longer moving average time base would be 
appropriate for precipitation data used to 
simulate streamflow drought periods. 
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VARIATION IN MICROCLIMATE AMONG SITES AND CHANGES OF CLIMATE 
WITH TIME IN BONANZA CREEK EXPERIMENTAL FOREST’ 


Leslie A. Viereck and Phyllis C. Adams 


Abstract.--The microclimate in Bonanza Creek Experimental 
Forest varies with slope and aspect and with the stage of 
forest succession. Examination of the climate record from 
nearby Fairbanks shows that there has been an increase in 
mean annual temperature since the mid-1970s. Climate 
warming would most likely result in a shift in distribution 
of permafrost and existing vegetation types with an 
expansion of the more productive forest types. 


Keywords: Boreal forest, taiga, permafrost, climate change, 
Alaska, forest succession. 


SITE DESCRIPTION AND VEGETATION 


The Bonanza Creek Experimental Forest (BCEF) is 
a 5045 ha research area located approximately 20 
km west of Fairbanks in interior Alaska. The 
area includes a section of the Tanana River 
floodplain at an elevation of approximately 120 
m and adjacent uplands rising to a ridge crest 
of 470 m. 


Upland forest types vary from highly productive 
aspen (Populus tremuloides Michx.), paper birch 
(Betula papyrifera Marsh.), and white spruce 
(Picea glauca (Moench) Voss) stands on 
south-facing, well drained slopes, to permafrost 
and moss-dominated black spruce (Picea mariana 
B.S.P.) forests of low productivity on north 
facing and lower toe slopes. Floodplain stands 
of balsam poplar (Populus balsamifera L.) and 
white spruce comprise productive forests on 
recently deposited river alluvium, where 
permafrost is absent; slow-growing black spruce 
stands and bogs occupy the older terraces, which 
sie underlain by permafrost (Viereck and others 
1983). 
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Successional communities are dominant within 
BCEF both in the uplands and on the floodplain. 
In the uplands there are three ages of forest 
resulting from wildfires in 1983, 1915, and 
approximately 1800. The youngest, early 
successional stands are primarily in the 
herbaceous, shrub, and sapling stages of 
recovery (Foote 1983); the middle aged stands 
(75-years-old), on upland slopes, are in a dense 
aspen and birch stage, often with an understory 
of white spruce; mature stands are of white 
spruce or a mixture of white spruce and paper 
birch. On poorly drained sites in the upland, 
and on older river terraces, fire in black 
spruce stands has usually resulted in the 
reoccurrence of black spruce stands, although on 
some sites black spruce has been replaced by 
paper birch (Foote 1983). 


On the floodplain of the Tanana River, erosion 
and silt deposition result in primary succession 
which begins with willows and alders and is 
followed first by balsam poplar and then by 
white spruce. Following one to several 
generations of white spruce, the floodplain 
sites develop permafrost and the productive 
white spruce stands are replaced by low 
productivity black spruce stands. 


SYNOPTIC REVIEW OF BCEF CLIMATE 


Climatic data reported here come from the 
National Weather Service observation station, at 
the Fairbanks International Airport (lat. 64 Onge 
N., long. 147°52' W.). The airport is on the 
floodplain of the Tanana River approximately 20 
to 25 km northeast of the BCEF and at an 
elevation of 132 m. Data collected at the 
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airport are most representative of the climate 
of the floodplain section of BCEF. 


The climate of BCEF is strongly continental and 
is characterized by temperature extremes from 
-50 to +35 ~C. The region lies within a rain 
shadow created by the Alaska Range. The 
physical barrier created by the mountains 
prevents the area from receiving precipitation 
from coastal storms and also results in rapid 
warming in winter as "chinook" type winds flow 
down the north slope of the mountains. The mean 
annual temperature of -3.3 ~C at Fairbanks 
results in the formation of permanently frozen 
soils (permafrost) on north-facing slopes and 
poorly drained lowlands. Thus the Forest is in 
a zone of discontinuous permafrost (Ferrians 
1965). At Fairbanks, strong temperature 
inversions occur 80 percent of the time during 
December and January, with gradients as steep as 
21 ~C/100 m of elevation during periods of 
extreme cold (Benson and Rizzo 1979). July is 
the warmest month with a mean daily temperature 
of 16.4 °C and January is the coldest with an 
average temperature of -24.9 °c. Because of its 
location at high latitude, BCEF experiences 
extremes of day length and sun angle which 
result in large differences in available solar 
radiation. At winter solstice, day length is 3 
hgurs, 42 minutes with a maximum sun angle of 
1-42', while at summer solstice there are 21 
hours, 50 minutes of sunlight and the maximum 
sun angle is 48°42', This results in average 
daily solar radiation of ,231 KJ/m /day in 
December and 22,375 KJ/m /day in June. 


The average annual precipitation at BCEF is 269 
mn. Most precipitation falls as rain in the 
summer months, a result of short-duration 
thunder storms and moist air masses that move in 
from the Bering Sea. Approximately 37 percent 
of the annual precipitation falls as snow from 
mid-October through April and remains as a 
permanent cover for 6 to 7 months each year. 
Maximum snow depths, averaging 75 cm, are 
commonly reached in February and March. The 
water equivalent at this time averages 11 cm. 
According to the Thornthwaite (1948) 
classification the climate of BCEF is semiarid, 
mesothermal, with little or no water surplus, 
and temperature efficiency normal to warm 
microthermal (D C'2de'2). 


VARIATION IN MICROCLIMATE AMONG SITES IN BONANZA 
CREEK EXPERIMENTAL FOREST 


Low sun angles, coupled with the continental 
climate, tend to make slope and aspect extremely 
important in the distribution of vegetation 
types. Permafrost also exerts strong controls 
over vegetation distribution by acting as a 
barrier to soil drainage, thereby creating wet 
or waterlogged soils. Presence or absence of 
permafrost is partially controlled by slope and 
aspect. These gradients of soil temperature and 
soil moisture are reflected in the distribution 
of plant communities and the productivity of 
forests and, in turn, result in a wide array of 


20 


microclimatic conditions within BCEF. The 
presence of permafrost on north-facing slopes, 
lower toe slopes, and old river terraces results 
in sharp contrasts with well-drained 
permafrost-free upland and floodplain soils. 
Consequently, sharp vegetation boundaries and a 
broad array of vegetation types are observed 
within close proximity to each other. Low 
productivity black spruce stands are formed on 
cold, wet permafrost-dominated soils and 
productive white spruce and successional 
hardwood stages occur on warm, mesic sites 
(Viereck and others 1983). The long hours of 
sunlight during the summer months offset to some 
extent the shortness of the growing season, 
which averages 111 frost-free days. 


Although sharp contrasts exist among vegetation 
types on different slopes and aspects, air 
temperatures on opposing north- and south-facing 
slopes at the same elevation do not differ 
significantly (Slaughter and Long 1975). Soil 
temperatures do differ appreciably and, in 
addition to resulting in the presence or absence 
of permafrost, greatly influence vegetation 
distribution. It has been shown that there is a 
good correlation between annual soil temperature 
sums and forest productivity (Viereck and Van 
Cleve 1984). 


In order to obtain more information on the 
effects of slope and aspect on microclimates, 
vegetation distribution, and productivity within 
BCEF we examined some microclimatic features of 
four stands along a gradient from hot and dry to 
cold and wet. These four stands were: (1) an 
open aspen-grassland with sage brush (Artemisia 
frigida Willd.) on a steep (75 percent) 
south-facing slope, (2) a mature white spruce 
stand on a 25 percent south-facing slope at an 
elevation of 396 m, (3) a mature paper birch on 
a 32 percent east-facing slope at 381 m, and (4) 
a black spruce stand on a 30 percent 
north-facing slope at an elevation of 400 m. 


The following parameters were measured in the 
four stands during the growing season: air 
temperature at 1.2 m in a standard weather 
shelter, soil temperature at depths of 5, 10, 
20, and 50 em, and soil moisture in the organic 
layer and at 5 and 10 cm in the mineral soil. 
For comparison of sites we used cumulative 
degree days for air and soil temperature. We 
present here a portion of the data collected, 
specifically for 1978 and 1979, which show some 
of the trends that were apparent along the 
environmental gradient. 


In figure 1, the monthly average air 
temperatures for the four sites are shown. 
Although there are some differences in extreme 
temperatures among the sites, the average 
monthly temperatures show very little 
difference. Mean annual temperatures for the 
four sites for the period April 1978 to May 1979 
were: aspen -0.08, white spruce -0,9, paper 
birch =-2.0, and black spruce -1.4 “C. During 
the same period the mean annual temperature at 
the Fairbanks Airport was -2.6 °c, Air 
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SITE e-c-e Aspen 
Figure 1.--Average monthly air temperatures for 
four stands in Bonanza Creek Experimental Forest 
for the period April 1978 through March 1979. 


temperature degree day sums (based on 5 7¢) also 
show no major differences among the four sites 
and Fairbanks. 


Soil temperatures did reflect slope and aspect. 
Temperature sums at 5 em (based on 0 ~C) for 
the period April 1978 through September 1979 
ranged from 3966 at the aspen stand to 922 in 
the black spruce stand. At 20 cm, temperature 
sums were 2393 at the aspen stand and only 539 
in the black spruce stand. Of the four stands, 
only the black spruce stand was underlain by 
permafrost. In figures 2 and 3 soil 
temperatures are shown for the four sites for 20 
em depth. 


Soil moisture also showed large differences 
among the sites. In the steep, south-facing 
aspen stand, during an extended dry period in 
July and August of 1978, soil moisture in the 
mineral soil layer remained below 10 percent for 
over a month and reached a minimum of 2 

percent. During the same period soil moisture 
was at 20 to 25 percent in the white spruce 
stand, and 25 to 60 percent in the birch stand. 
The mineral soil in the black spruce site 
remained frozen but the organic soils overlying 
the permafrost had a moisture content of between 
440 and 700 percent (based on oven-dry weight). 


It has also been shown that the snowpack and the 
time of snow melt differed considerably at these 
four sites (Slaughter and Viereck 1986). 
Snowpack was deepest on the black spruce site 
and remained nearly a month and a half longer 
than at the aspen site, which was snow-free by 
the end of March. 








SITE e-e-o Aspen +++ White Spruce ++ « Birch 66 Black Spruce 


Figure 2.--Monthly average soil temperatures at 
20 cm depth for four stands in Bonanza Creek 
Experimental Forest for the period April 1978 
through March 1979. 
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Figure 3.--Soil temperature sums (based on 0 “¢) 
for 20 em depth in four stands in Bonanza Creek 
Experimental Forest for the period April 1978 
through March 1979. 


CHANGES IN MICROCLIMATE WITH SUCCESSIONAL 
CHANGES 

1. Upland Secondary Succession Following Fire 
There are significant changes in soil 


temperatures with time in both upland and 
floodplain succession. In previous studies (Van 
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Figure 4.--A, Forest succession following fire 
on upland south-facing slopes in BCEF. Numbers 
in boxes refer to "turning points" to be studied 
under the LTER program. B. Soil degree day sums 
and organic layer thickness for the seven 
successional stages in 12a (modified from Van 
Cleve and Viereck 1983). 


Cleve and others 1980, Van Cleve and Viereck 
1983) we have shown that one result of fire in 
spruce stands in interior Alaska is warming of 
the soil profile. This increase in temperature 
results from a removal of the moss and dead 
organic layer and a change in the surface 
albedo. These changes in the soil surface allow 
for more absorption of heat. Figures 4A and 4B 
illustrate successional stages of the forest 
following fire on south aspect slopes with 
corresponding changes in organic layer thickness 
and soil degree day sums. As forest succession 
develops following fire, the organic layer on 
the forest floor also thickens, and soil 
temperatures gradually decrease. A significant 
turning point in the successional sequence is 
the establishment of feathermosses on the forest 
floor under spruce. The insulating effect of 
the feathermosses results in lower soil 
temperatures, and in black spruce sites, the 
reestablishment of a shallow thaw layer (active 
layer) above the permafrost. In our LTER study 
we will examine soil temperatures and soil 
moisture at three critical points in the 
successional sequence leading to the development 
of productive white spruce stands in the uplands 
of BCEF. These are (1) early succession, 
typified by the colonizing species which arrive 
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Figure 5.--A. Forest succession on the 
floodplain of the Tanana River. Numbers in 
boxes refer to "turning points" to be studied 
under the LTER program. B. Soil degree day sums 
and organic layer thickness for the twelve 
successional stages in 13a (modified from Van 
Cleve, Dyrness, and Viereck 1980). 


within 5-10 years following fire; (2) 
mid-succession turning point, characterized by 
the shift from aspen and birch to spruce 
dominance with the establishment of the 
feathermoss layer on the forest floor; and (3) 
late succession, represented by a mature white 
spruce stand with a well developed organic layer 
(figure 4A). 


2. Primary Succession of the Tanana Floodplain 


Similar trends in soil temperatures with time 
have been shown for the successional sequence 
that develops on the floodplain of the Tanana 
River following the deposition of new sand 
bars. Twelve stages of forest development on 
the floodplain have been described (Viereck 
1989). The sequence begins with the 
establishment of willows on recently deposited 
alluvium, proceeds through an alder shrub stage 
to a balsam poplar stage which, in turn, is 
Slowly replaced by stands of white spruce. As 
feathermosses develop beneath the spruce canopy, 
the organic layer increases in thickness and 
soil temperatures become lower (figures 5A and 
5B). In later stages the white spruce is 
replaced by black spruce and permafrost develops 
in the soils. At 10 em depth in the soil the 
degree day (DD) sums decrease from approximately 
1500 DD in early stages to 480 DD in the black 
spruce stands. 
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In the LTER study we will examine changes in 
soil temperature and other microclimatic factors 
at four critical points. These are (1) newly 
colonized alluvial deposits, with young willows 
and herbs; (2) dense willow and alder stands in 
transition to closed balsam poplar stands; (3) 
decadent balsam poplar intermixed with 
canopy-height white spruce and a developing 
feathermoss layer; and (4) mature white spruce 
stands with intermittent permafrost. 


LONG TERM RECORDS AND CLIMATIC VARIABILITY 


Climatic data used for examining long term 
trends were obtained from the National Weather 
Service station at the Fairbanks International 
Airport because there are no long term climatic 
records from BCEF. Climatic records for 
Fairbanks go back to 1907 but a reliable record 
extends back to only 1917. However, several site 
changes since 1917 make analysis of long term 
trends suspect. Two papers in a recent 
symposium on the effects of climate changes in 
Alaska (Bowling 1984, Juday 1984) analyzed long 
term climatic records from Alaska. Both authors 
showed that 1970-1980 has been warmer than 
previous decades for much of Alaska, but neither 
presented conclusive evidence about whether the 
changes are cyclic or indicate a significant, 
long-term warming trend. Most models of 
predicted climatic change resulting from the 
"greenhouse effect" predict greater warming at 
northern latitudes than at middle latitudes. 
Recently, Bowling has presented information 
which show that winters in Fairbanks have b¢en 
well above the 30-year average since 1975. 

She also has shown that the mean annual 
temperature of three cities in Alaska has been 
above average for the same period. 


The mean annual temperature for Fairbanks for 
the 70 years of the reliable record, 1917 to 
1987 is -3.1 ~C, close to he 30-year average 
from 1951 to 1980 of -3.3 ~C. In contrast the 
mean annual temperature for the decade 1979 
through 1988 is -1.8 -C. If the warming trend 
of the past decade continues, the new 30-year 
average, from 1961 to 1990, will be considerably 
higher than the previous 30-year average. 
Although the increase in mean annual temperature 
may be small, because it is close to O C it 
could have significant effects on permafrost 
distribution and subsequently on vegetation. 


We have analyzed the National Weather Service 
climate records from the Fairbanks International 
Airport from the year following the last site 
change, 1951 through 1988. Data sets of average 
annual temperature were created to examine 
temporal variation. Annual temperatures at the 
Fairbanks Airport from 1951 to 1988 are quite 
variable with a standard deviation of 16.2 
around a mean of -3.0 ~C. The linear trend over 
this time period indicates an increase of 0.06 


Fairbanks Daily News Miner, Sunday, January 
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Figure 6.-~Average yearly temperatures at the 
Fairbanks Weather Service station for the period 
1951 to 1987 with a fitted quadratic regression 
curve. 


“6 per year while the quadratic regression curve 
fit to the annual temperatures shows a general 
increase in mean annual temperatures (figure 

6). To further examine this trend, we subjected 
this data set to smoothing techniques. Unlike 
regression, smoothed data points are derived 
only from nearby points, and remove the 
influence of a single extreme value. These 
smoothing algorithms are described by Tukey 
(1977) and were applied as follows: repeated 
medians of 3 (with end value smoothing), 
splitting peaks and valleys, splitting peaks and 
valleys again, the Hanning algorithm, and 
medians of 3. Smoothing was done by means of a 
SAS macro program (SAS Institute 1985) written 
by Janet Nelson at the Forestry Sciences 
Laboratory, Portland, Oregon. The time series 
plot of this smoothed curve (figure 7) shows 
that the average yearly temperature has 
increased since the mid-1970s. Some of this 
increase may be due to building construction and 
other human activities around the airport which 
has resulted in the extension of the Fairbanks 
"heat island" to the airport environs in recent 
years (Bowling and Benson 1978). 


We also created data sets of mean temperature 
values for winter months and for summer months. 
Winter data consisted of temperatures for 
December, January, and February and summer data 
included June, July, and August. These data 
sets were subjected to the smoothing techniques 
described above and plotted as time series. 
Plots of the average winter temperatures show a 
steady increase from a low in the late 1960s 
(figure 8). The summer temperatures show less 
variability than winter temperatures, with a 
noticable increase beginning in the mid 1980s 
(figure 9). 


Another long term record of interest for the 
Bonanza Creek site is the record of ice break-up 
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Figure 7.--Average yearly temperatures 
(smoothed) at the Fairbanks Weather Service 
station for the period 1951 to 1987. 


on the Tanana River at Nenana, about 50 km 
down-river from the LTER floodplain site (figure 
10). The exact month, day, minute, and second 
of break-up has been determined for an annual 
gambling pool. The river break-up seems to 
integrate a number of climatic factors including 
severity of past winter and depth of winter snow 
accumulation. In addition, weather conditions 
during April have a strong influence, making 
predictive modeling of the break-up date 
difficult. Although we have not analyzed these 
data for trends, during the warmer period of the 
last decade there have been an above average 
number of "early" (on or before May 1) break-up 
dates. 


PREDICTIONS OF CLIMATIC VARIABILITY AND 
ECOSYSTEM RESPONSE 


The principal objective of this workshop was to 
discuss ways in which the LTER ecosystems might 
respond to climatic variability. In a previous 
paper, the senior author and Van Cleve (Viereck 
and Van Cleve 1984) discussed the possible 
effects of a predicted climatic warming of 5 ~C 
on the taiga ecosystems of interior Alaska. The 
following discussion is largely derived from 
that paper. 


There are a number of possible scenarios that 
might occur with a climatic warming in interior 
Alaska, but we discuss only two here: one in 
which a mean annual increase is spread 
throughout the year with no change in 
precipitation, and a second in which most of the 
increase in temperature is in the winter and and 
there is an increase in summer precipitation. 
Most models predict an increase in precipitation 
along with a temperature increase as open water 
replaces sea ice. Also, a warmer period in the 
northern Yukon Territory between 11,000 and 
8,900 BP was reported to be a wetter period than 
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Figure 8.--Average summer (June, July, and 
August) temperatures (smoothed) for the 
Fairbanks Weather Service station for the period 
1951 to 1987. 


present (Cwynar 1982). 


A general warming of 5 °c will probably have 
little overall effect on the species present or 
the composition of the existing plant 
communities in the taiga in interior Alaska. 
Pollen records show no discernible change in 
past warmer and colder periods (the warmer 
Hypsithermal and the colder Little Ice Age) for 
interior Alaska (Ager 1975). Increased 
temperatures would, however, most likely change 
the distribution of vegetation types. 


Higher summer temperatures with no change in 
precipitation would likely increase the 

frequency of naturally occurring wildfires. 

This in turn would increase the proportion of 

the vegetation in early stages of succession and 
decrease the area of mature stands of white and 
black spruce. Paper birch would invade many of 
the colder sites now dominated by black spruce. 
With a change toward a warmer and drier climate, 
an expansion of the steppe-like grass type and 
aspen types would be expected. The closed spruce 
forest would expand at the expense of the open 
forest types on wet sites. Open forest types, in 
turn, would invade some of the bog sites. 


One significant effect of climatic warming would 
be the thawing of permafrost on some sites and 
an increase in the active layer thickness on 
others. This would result in the expansion of 
forest types now found on permafrost-free sites 
and those with a thick active layer. An 
increase in soil temperatures, especially on the 
colder sites, would result in a more rapid 
turn-over of nutrients and increased ecosystem 
productivity. 


It is more difficult to predict the changes 
resulting from climatic warming if there is 
little change in summer temperatures but an 
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Figure 9.--Average winter (December, January, 
and February) temperatures (smoothed) for the 
Fairbanks Weather Service station for the period 
1951 to 1987. 


increase in precipitation. With an increase in 
precipitation the frequency of fire should 
decrease; early successional stages of 
vegetation would decrease while mature 
vegetation types would increase. Increased 
summer precipitation might also result in an 
expanded moss layer with a shallower depth of 
thaw on some permafrost sites, resulting, in 
turn, in an expansion of open black spruce 
forests and some bog and fen types. On the 
other hand, increased precipitation could result 
in higher nutrient turnover and more biological 
activity in some soils. The overall trend would 
be toward higher forest productivity and larger 
standing crops of biomass and nutrients. 
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Figure 10.--Break-up dates for the Tanana River 
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CLIMATIC VARIABILITY AND SALT MARSH ECOSYSTEM RESPONSE: 


RELATIONSHIP TO SCALE 
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Abstract. Small scale spatial variability in rainfall, seasonal 
changes in wind velocity and directionality, and hydrometeoro- 
logically induced intrusions of freshwater are characteristic 
features of the North Inlet Estuary. Precipitation resulting from 
hurricanes or tropical storms is an integral component of the 
regional water balance. The 1986 drought was a regional scale 
phenomenon which significantly altered estuarine productivity. 


Keywords: Crassostrea virginica settlement, drought, fish 
community structure, hurricanes, primary production, Spartina 
alterniflora, tropical storms, water budget 


INTRODUCTION 


The North Inlet marsh-estuarine system is 
located along the northeast-southwest oriented 
coastline, 70 km northeast of Charleston, South 
Carolina. Research facilities are located on 
Hobcaw Barony (a 7085 ha tract of maritime forest; 
Fig. 1) 8 km east of Georgetown, SC. The primary 
research areas are a 2,630 ha high-salinity 
Spartina alterniflora marsh and 715 ha of tidal 
creeks and intertidal flats which are separated 
from the Atlantic Ocean by sandy barrier islands. 
The estuary is bordered on the west by loblolly 
and long leaf pine forests. Hydrographic 
characteristics include an annual average seasonal 
salinity range of 30 to 34 ppt (monthly mean 
salinities average 19 to 36 ppt), average channel 
depth of 3 m, and a seasonal water temperature 
range of 3° to 33° C). Wetland habitats include 
exposed and sheltered sandy beaches; intertidal 
flats and oyster beds; submerged macroalgal mats; 
sand, shell, and mud benthic habitats; shell 
middens; and bird rookery islands. At mean tide, 
Spartina alterniflora comprises 73.0 percent, 
tidal creeks 20.6 percent, oyster reefs 1.0 
percent and exposed mud flats 5.4 percent of the 
marsh-estuarine zone (Dame et al., 1986). More 
than 1200 ha of brackish and freshwater marshes 
‘border the Winyah Bay side of Hobcaw Barony. 


lPresented at the LTER Workshop on Climate 
Variability and Ecosystem Response, Mountain 
Research Station, University of Colorado, August 
21-24, 1988. 
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University of South Carolina, Columbia, SC 
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Figure 1. Hobcaw Barony and the North Inlet 
Estuary (Rain gauges R1-R8; Haulover Creek, HC; No 
Man’s Friend, NMF; South Jones Creek, SJ; Town 
Creek, TC). 
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The North Inlet temperate-subtropical climate 
is greatly influenced by the low elevation of this 
coastal site. Temperatures are moderated by the 
proximity of the Atlantic Ocean and the Gulf 
Stream which produce both relatively higher 
temperature minima and lower temperature maxima 
than would be found farther inland. The mean 
annual temperature (1951-1980) is 18.06° C. 
Temperatures range from an average of 8.44° C in 
January to 26.94° G in July. Record high 
temperature was recorded July 10, 1977 at 40.56° ¢ 
and record low of -11.67° C on December 11, 1962. 
The number of days with temperatures 32.22° ¢C (90° 
F) and above average 45 per year and days with 
temperatures below 0° C average 35 days. The mean 
number of frost free days is 253, extending from 
mid-March to mid-November. 


Precipitation in the study area (1951-1980) 
averaged 130 cm per year (NOAA, 1985). Annual 
precipitation patterns are highly variable due to 
the episodic occurrence of frontal passages, 
tropical storms and hurricanes. Storm size and 
frequency for a given season are quite variable. 
Winter and spring are drier (20 and 21 percent of 
annual precipitation) averaging 4.8 and 4 storms 
per month, respectively. Fall is wetter (24 
percent of annual precipitation) with 4 storms per 
month. Summer is the wettest season (35 percent of 
annual precipitation) with the greatest variation 
in storm frequency and size being due to the 
frequent occurrence of tropical storms and 
hurricanes. On an annual basis, over 30 percent of 
the storms result in less than 1 cm precipitation, 
but comprise less than 5 percent of the total 
annual rainfall. Because of the high probability 
of tropical storms, the return frequency for 
storms >10 cm in 24 hours is only 2 years (Purvis 
and McNab, 1985). The average (1935-1986) water 
budget for Georgetown indicates a slight deficit 
(1 to 17 mm) occuring from April to August. 


North Inlet is the only tidal pass between the 
marsh-estuarine system and the ocean. Four smaller 
connections (South Jones, No Man’s Friend, Sign, 
and Haulover Greeks) at the southern end of the 
system support tidal exchange with the Winyah Bay 
Estuary (Fig. 1). Despite its relatively small 
surface area (30 square miles) Winyah Bay has the 
largest estuarine drainage area in the 
southeastern US (9,511 square miles; Dowgiallo et 
al, 1987). The primary water exchange occurs 
between Town Creek at North Inlet and the ocean 
(79 percent), but inputs from Winyah Bay occur 
during southwesterly winds and/or during high 
river discharge. Twenty-one percent of the water 
exchange in the North Inlet Estuary occurs through 
South Jones and No Man’s Friend (Kjerfve et al., 
1982). Of that volume, 80 percent of the exchange 
occurs through South Jones (Kjerfve, 1978). Sign 
Creek flows into the lower portion of South Jones 
Creek. Haulover Creek is narrow and shallow, 
therefore, restricting flow to within a few hours 
of high tide (Michener and Allen, 1982). Exchange 
between the North Inlet and Winyah Bay estuaries 
is restricted and opposing tidal flows meet at 
relatively fixed nodal points under most 
conditions (Schwing and Kjerfve, 1980). 


Major factors affecting southeastern coastal 
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environments include the recurrence of large 
rainfall deficits (droughts) and rainfall excesses 
(tropical storms and hurricanes). Precipitation 
conditions affect the water budget and primary 
production, as well as vertebrate and invertebrate 
communities. Objectives of this paper include 
examination of spatial variation in rainfall, 
temporal variation in wind speed and 
directionality, hurricane frequency, and effects 
of moisture conditions on the water budget and 
primary production. We will characterize the 1986 
drought and examine its impact on oyster 
settlement and fish community structure. In 
addition, we will examine the conditions which 
lead to intrusion of Winyah Bay waters into the 
North Inlet Estuary and the related changes in 
nutrient availability. 


METHODS AND MATERTALS 


Meteorological Data 


Real-time monitoring and archival of 
meteorological conditions at North Inlet have been 
implemented with a meteorologic system designed by 
the Climatronics Corporation and software 
developed in-house. The Climatronics system is 
designed to measure, transmit, and store wind 
speed, wind direction, air temperature, barometric 
pressure, solar radiation, precipitation, water 
level, water temperature and water conductivity 
data. The system is comprised of a 10 meter 
aluminum tower, individual sensors, mainframe 
with power supply and modem, a multiplexing 
recorder, and a microprocessor-based cassette data 
acquisition system. The software scans, edits, 
reformats and summarizes these data by hour for 
subsequent archival. The meteorological station 
was deployed June 1982 at Oyster Landing dock in 
Crab Haul Creek (Fig. 1). 


Belfort Universal Series gauges (model #5-780) 
were used to record rainfall amounts (inches) at 
various locations on Hobcaw Barony (Fig. 1). All 
gauges were situated with at least a 45 opening in 
the forest canopy above the gauge. Six gauges (Rl, 
R2, R5, R6, R7, R8) recorded data from 1978 
through 1982. 


Water Budget Calculations 


Water surpluses and deficits, determined from a 
monthly water balance, are fundamental to the 
comprehension of important hydrologic variables 
such as flood duration, runoff rates, soil 
moisture changes, and freshwater inflows to bays 
and estuaries. The water balance components 
developed by Thornthwaite and Mather (1955) 
provide quantitative estimates of these variables 
in a form that can be easily compared with 
biological or chemical processes occurring in a 
watershed. There are seven components in the 
Thornthwaite and Mather water budget: (1) 
potential evapotranspiration, PE; (2) soil 
moisture storage, ST; (3) actual 
evapotranspiration, AE; (4) precipitation, P; (5) 
deficit, D; (6) surplus, S; and (7) runoff, RO. PE 
is defined as the maximum amount of 
evapotranspiration that would be possible if the 





soil surface was covered by vegetation and there 
were no soil moisture shortages within the root 
zone. PE estimates are calculated from the 
empirical relationships between air temperatures, 
daylight hours, and the heat index (Thornthwaite 
and Mather, 1957). 


Calculations of surpluses and deficits were 
based on the daily temperature and precipitation 
data for the Georgetown weather station supplied 
by John Purvis at the Office of State Climatology. 
The North Inlet temperature data were taken from 
the Climatronics meteorological system at Oyster 
Landing. Precipitation at Oyster Landing was 
measured daily at 10:00 am with a standard Belfort 
5-700 Series rain gauge. The sandy soils along the 
coastal regions are generally well drained and 
were assumed to have very high soil moisture 
storages (16 inches). 


Nutrient Data and Winyah Bay Intrusion 


Water samples are collected daily at 10:00 am 
(eastern standard time) from one foot below the 
surface at Glambank Creek and at Oyster Landing 
adjacent to the meteorological station. Samples 
are transported immediately to the laboratory. 
Approximately 200 ml of each sample is filtered 
through precombusted Whatman GF/F glass fiber 
filters. Ammonia and nitrate/nitrite are preserved 
with phenol and mercuric chloride, respectively, 
stored below 4° C and analyzed within two weeks to 
minimize loss due to sample storage. Aliquots for 
total nitrogen and phosphorus were frozen in 
digestion tubes and stored separately. Ammonia, 
nitrate/nitrite, and total N and P were measured 
by automated colorimetric tests with a Technicon 
Autoanalyzer II or Orion Scientific Autoanalyzer 
System. Ammonia was determined by the phenate 
method (Technicon Industrial Method No. 154-71W) 
and nitrate/nitrite by cadmium reduction 
(Technicon Industrial Method No. 158-71W). Total N 
and P were determined using Technicon Industrial 
Method No. 329-74W following persulfate digestion 
(Glibert and Loder, 1977). 


Nutrient concentration data were log 
transformed to achieve normality and induce 
homogeneity of variances (Sokal and Rohlf, 1981). 
Regression analysis (SAS, 1985a) was then 
performed on the transformed data to remove any 
seasonal bias. Harmonic regression terms 
(Chatfield, 1984) which adjusted for seasonal 
(365.25), lunar (29.0), and tidal (14.5) 
periodicities were included in the initial model 
(R? = 0.05, p=0.0001). Exploratory analyses 
included principal component analysis and stepwise 
regression. Graphical analysis of residual error 
terms were utilized to eliminate problems 
associated with collinearity of predictor 
variables and to arrive at the reduced model used 
in the final analysis (Draper and Smith, 1981; 
SAS, 1985a). Analysis of variance was performed 
on the reduced model to examine the influence of 
individual parameters and their interaction on 
residual nutrient concentrations during 1983 (SAS, 
1985a). Frequency analysis of salinity data 
collected at Clambank Creek (1981-1987) indicated 
that 25 percent of the recorded salinities were 





less than or equal to 31 ppt (SAS, 1985b). For 
this analysis, we therefore defined low salinity 
at Clambank Creek to be at or below 31 ppt. 
Further, a low salinity event was defined as any 
period from one to multiple consecutive days 
during which low salinities were recorded. 


Oyster Settlement 


American oyster (Crassostrea virginica) spat 
settlement was monitored at three stations (Town 
Creek, Oyster Landing, and Old Man Creek; see Fig. 
1) from 1983 to 1986 (Kenny et al., 1989). 
Salinity typically remains high at Town Creek and 
Old Man Creek. Current velocities are generally 
lower at Oyster Landing and the site receives 
periodic freshwater runoff from the nearby forest 
and salinities may be depressed for extended 
periods. 


At each station two to four vertically oriented 
rope harnesses were used to maintain 225 em”, 5 mm 
thick asbestos plates in a horizontal plane at 
three different levels: 30, 70 cm above mean low 
and 30 cm below mean low tide. The harnesses were 
suspended from PVC pipes which were located 
approximately 3 m from the bank and concrete 
blocks held the arrays in place. Plates at each 
level were approximately 30 cm apart. One side of 
each plate was slightly textured and was placed on 
the top side throughout the study. The plates were 
replaced every 2 weeks (14 +/- 2 days) in the 
settlement season and every four weeks in winter. 
They were transported to the lab in upright 
slotted trays to prevent abrasion. If necessary 
the plates were gently rinsed with water to remove 
excess sediment. Counts of spat were made on the 
entire top and bottom sides of all plates with a 
dissecting microscope (5X). Plates were wire 
brushed under running water prior to reuse. 


Analysis of variance (SAS, 1985a) was performed 
on log transformed count data (logl0 (total 
harness count+l)) utilizing a complete block 
design for the statistical model (Sokal and Rohlf, 
1981). Tukey’s Studentized Range Test was used for 
all multiple comparison tests (Mize and Schultz, 
1985). Post hoc reduction of the number of 
multiple comparisons examined was not performed to 
ensure that all statistical interpretations were 
conservative. 


Fish Community Structure 


Biweekly seine collections of fishes, shrimps, 
and crabs from a tidal pool were initiated April 
15, 1983 and continued through the present. The 
sampling site (Fig. 1; adjacent to OL) is located 
near the meteorological station at Crabhaul Creek 
and consists of a depression in the creek bed, 
which is isolated during low tide. The pool is 
approximately 10 m in diameter with a maximium 
depth of 1 m. The bottom substrate is primarily 
muddy with scattered oyster shells. The pool is 
surrounded by Spartina alterniflora along the 
steepest banks and live oyster reefs occur at the 
pool’s inflow and outflow. Salinity ranges from 1 
ppt during heavy rains to 33 ppt during droughts, 
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extreme high tides, or strong NE, E, or SE winds. 
Water temperatures range from 3° C during the 
winter to 38° CG during summer months. The tidal 
creek extends to the edge of a loblolly pine 
forest and merges with an intermittent stream 
which contributes freshwater following major 
rainfalls. 


Salinity and air and water temperatures were 
recorded prior to the seine collections. Two 
sequential hauls with a 0.25 in. bag seine that 
stretches across the width of the pool comprise 
the biweekly sampling effort. All fishes, shrimps, 
and crabs are identified, enumerated, weighed, and 
up to 100 individuals of each taxon are measured 
from each seine haul. 


RESULTS AND DISCUSSION 
Wind 


Fall (September, October, November) and winter 
months (December, January, Febuary) are dominated 
by northeasterly winds (mean = 58 and 37 percent, 
respectively; Fig. 2). Spring months (March, 
April, May) exhibit the greatest variation in 
average wind directions (NE = 22 percent, NW = 18 
percent, SE = 22 percent, SW = 37 percent). 
Southwesterly winds dominate the summer months 
(June, July, August; 49 percent). On an annual 
basis, northeasterly and southwesterly winds are 
most common (each accounting for 34 percent). 
Northwesterly and southeasterly winds are less 
frequent (19 and 13 percent, respectively). 
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Figure 2. Average percentage of winds originating 
from four quadrants on a monthly basis. 


Highest seasonal average wind velocities (daily 
mean) consistently originate from the northeast 
(fall, 3.8 m/s; spring, 4.4 m/s; summer, 4.1 m/s; 
winter, 3.7 m/s). Lowest seasonal average wind 
velocities typically originate from the 
southeastern quadrant (fall, 2.2 m/s; spring, 3.4 
m/s; winter, 2.6 m/s). During the summer, the 
lowest average wind velocities are associated with 
northwesterly winds. The greatest seasonal 
variation in wind velocity is associated with 
southeasterly winds during the summer and winter 
months (C.V. = 55 and 56, respectively). The 
passage of tropical storms and hurricanes 
partially explain this variability. The lowest 
variability (C.V. = 19) is associated with summer 
northwesterly winds. The highest daily average 
wind velocity (13.15 m/s; SE) was recorded July 
18, 1984 and was related to the passage of a gale. 
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Precipitation Variation 


Rainfall events resulting in 2.54 ecm 
precipitation or less are typically associated 
with low spatial variability (Fig. 3). Larger 
rainfall events resulting in 2.54-5.08 cm of 
precipitation exhibit higher spatial variability. 
These events occurred most frequently during 
spring months and were usually associated with 
high and variable winds. Spatial variability 
steadily decreased with increased precipitation 
(5.08-15.24 cm). These events were marked by the 
passage of large frontal systems. 
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Figure 3. Mean daily rainfall versus daily 
coefficient of variation (CV). 


Highest spatial variability, when analyzed on a 
monthly basis, was associated with late spring 
through early fall events (Fig. 4). During these 
months, localized thunderstorms were frequently 
observed with variability in wind velocity, 
direction, and duration influencing the spatial 
variability of precipitation during these storms. 
From October through April lower spatial 
variability is associated with the passage of 
large fronts. These fronts typically cover a 
greater area and are of longer duration than the 
isolated thunderstorms of summer. 
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Figure 4. Mean monthly rainfall versus monthly 
coefficient of variation (CV). 


Hurricanes and Tropical Storms 


Tropical storms or hurricanes affect the South 
Carolina coast approximately once every 2.6 years 
(Gentry, 1971). From 1901-1985, 22 tropical storms 
or hurricanes have made landfall on the South 
Carolina coast (Purvis et al., 1986). Only eight 
of these were class 2 to class 4 intensity. No 
category 5 hurricanes have hit the South Carolina 
coast during this period. Two storms (September 
17, 1945; September 29, 1959, Gracie) were 
classified category 3 and Hazel on October 15, 
1954 was a category 4 storm. Hazel made landfall 


near Little River, South Carolina with 106 mph 
winds and tides up to 16.9 feet. Gracie made 
landfall on St. Helena Island and maintained 
hurricane strength more than 100 miles inland. 
Damage occurred from Beaufort to Charleston with 
heavy rains. Prior to 1900 three hurricanes with 
direct impact to the North Inlet system occurred 
in 1822, 1854, and 1893 (Purvis et al., 1986). The 
1854 storm caused severe flooding with the 
Waccamaw Neck inundated with salt water from 
Waverly to Pee Dee. The great storm of 1893 
occurred August 28 with winds approximately 120 
mph. The storm struck at high tide and submerged 
whole islands. Two other hurricanes that affected 
the Georgetown area occurred on October 3 and 13 
in 1893. Heavy rains were associated with both 
storms. 


Hurricanes or tropical storms are an annual 
threat to Georgetown County. For any particular 
year there is a 46 percent chance of a tropical 
cyclone (34 or more knots wind speed) and a 13 
percent chance of 2 such storms occurring in a 
given year (Purvis and McNab, 1985). The 
calculated return frequency of a 14 foot storm 
surge is approximately 100 years (Myers, 1975). 
The hurricane season along coastal South Carolina 
is June to October with the maximum storm 
occurrence frequency from September 24-30 (Purvis 
and McNab, 1985). The predominant storm direction 
is from the southwest. 


During the period 1935-1986, 31 tropical storms 
or hurricanes passing within 75 miles of 
Georgetown resulted in precipitation at the 
Georgetown National Weather Service monitoring 
station. The percentage of Georgetown’s total 
annual precipitation due to tropical storms or 
hurricanes ranges from 0 to 24 percent (Fig. 5a). 
During years when these storms occur, they are 
responsible for an average of 12.3 percent of the 
total annual rainfall. Total annual precipitation 
averages 128.8 cm during years affected by 
tropical storms or hurricanes versus 120.0 cm for 
those years not affected (Fig. 5b). 


Water Budget Results 


The water balance is not consistent from year 
to year. Although North Inlet usually has high 
winter surpluses and high summer deficits, there 
can be periods when precipitation exceeds PE 
during summer (1985; Fig. 6) and when PE exceeds 
precipitation during winter (1984-85, Fig. 6). 


Drought conditions at North Inlet are not the 
norm. When comparing the average water budget at 
North Inlet with that of Georgetown (8 km away), 
it was found that the winter surpluses were 
approximately the same; however, the deficits of 
summer were far more extensive at North Inlet. 
Differences between sites were a function of the 
time period that was used in our comparisons 
(1937-1986 for Georgetown, 1979-1986 for North 
Inlet) and not differences due to location. 
Regression analysis of precipitation data revealed 
a highly significant correspondence between 
Georgetown and North Inlet (R*=0.883). 


25 


© 20 | 
o | | 
= 18 i nl ! 
f 1 | 
o oO: i i ' : 
ao PF HW 
| i 5 i | 
| 
4 
0 
200 
& 
2 
c 
= 150 
7 
a 
a 
& 
® 100 
Dees 
o. 


50 
1935 1940 1945 1950 1955 1960 196519701975 1980 1985 


Year 


Figure 5. (A) Percentage of annual precipitation 
due to hurricanes and tropical storms (first 
storm=black, second storm=open). (B) Total 
annual precipitation for Georgetown, SC. 
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Figure 6. Water budget for Hobcaw Barony and 
North Inlet, SC. 


Primary Production Measurements 


North Inlet marsh grass, Spartina alterniflora, 
production and biomass have been measured by Dame 
and Kenny (1986) using the harvest technique and 
by Morris and Haskin (1988) using the census 
method. The two techniques produce similar 
estimates of standing biomass that range from a 
low of 50 g/m? in midwinter to a high of 1000 g/m? 
in early fall. Spatial and temporal differences of 
live standing biomass can however, be very 
significant. Dame and Kenny (1986) found that 
maximum creekside biomass ranged from 500 g/m? to 
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1150 g/m? during October to November, while 
maximum highmarsh biomass (i.e. areas not 
frequently inundated by tides) ranged from 280 
g/m? to 350 g/m? during September to October. 
Spatial and temporal differences in biomass can be 
accounted for by differences in hydrodynamic 
forcing and the energy subsidies associated with 
tides and flooding frequencies (Odum et al. 1979, 
Conner et al. 1981). 


Relationships between variation in marsh 
biomass and climatic events such as El Nino, are 
implied in the work of Dame and Kenny (1986). 
Climatic variations can produce annual differences 
in marsh salinities, hyrodology, and nutrient 
distributions. According to Dame and Kenny (1986), 
maximum biomass in the North Inlet creekside marsh 
occurred in 1983 (1150 g/m?) and 1984 (1100 g/m?) 
when winter surpluses of rainwater were extremely 
high and when mean winter salinities were 
extremely low (15 ppt). What must also be included 
in this climatic equation however, is the power of 
the wind to move river water past the tidal node 
thereby altering the spatial and temporal 
distribution of nutrients, sediments, and 
productivity in North Inlet. This is discussed 
later in more detail. 


The standing stock of marsh grass is only 
indicative of marsh productivity. A more direct 
measurement of marsh production in North Inlet was 
done by Morris and Haskin (1988). They followed 
the growth and decay of every single blade of 
grass in six 100 cm? plots. They found a maximum 
production rate of 270+/-87 g/m* in the midmarsh 
habitats of North Inlet in July of 1985 (Fig. 7). 
They also found that production rates peak one or 
two months before peak standing stocks, 
productivity has been decreasing since 1985, and 
annual trends in production at different sites are 
parallel. Although no climatological interactions 
are analyzed, Morris and Haskin conclude: 
"Parallel growth trends between sites, both annual 
and seasonal, suggest that productivity is 
responding to at least a regional scale 
phenomenon." We see good evidence for their 
conclusion. For example, salinity data at North 
Inlet indicated bay water intrusions during 
periods of maximum production, water quality data 
in the bay indicated high nitrogen concentrations 
(55 ug At./1), and water budget data indicated 
water surpluses during maximum production with 
increasing drought conditions after 1985. Other 
factors such as high phosphorus concentrations in 
the creek waters or variations in sea level and 
sediment salinity (Morris 1988) during 1985 could 
explain the high productivities, but these 
parameters are also indirectly controlled by 
climatic variables. 


Drought 


Droughts are a normal part of the climate in 
the southeast and have occurred 17 times in the 
past 100 years (Guttman and Plantico, 1987). 
Drought conditions occurred during 20 to 25 
percent of the months from 1895 to 1986. Droughts 
in South Carolina last from a year and a half to 
seven years. Major droughts occurred in the 
1920's, 1930’s, 1950’s, and late 1980's. The most 
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Figure 7. Spartina alterniflora monthly primary 
production estimates + 1 S.E. 
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persistent drought extended from 1952 to 1957. The 
end of the 1950’s drought was marked by a rapid 
shift to wet conditions which persisted until the 
late 1970's. The 1980’s brought a return to drier 
conditions. 


The 1986 drought affected the entire southeast 
coast. Although not unusually long (one year), 
this drought was the most severe on record since 
precipitation records began in 1876 (Dowgiallo et 
al., 1987) and especially affected the Carolinas 
and Georgia. During April (the month most 
deficient in rainfall for the southeast), 
Charleston, SC had a 73 percent increase in the 
number of clear days and the number of clear days 
for the year was greater than the long-term 
average (Dowgiallo et al., 1987). Air temperatures 
during 1986 averaged 3° F above normal in 
Charleston, SC (Dowgiallo et al., 1987). In 
addition, Charleston had 83 days with maximum 
temperatures above 32.22° C for 1986 compared to 
the long-term average of 49. 


Monthly precipitation was consistently below 
average from January through July and the maximum 
deviation was recorded in July (107 mm; Fig. 8). A 
slight positive deviation from normal was observed 
in August and was related to the passage of a 
single storm. Drier conditions returned in 
September and October. Air temperatures were above 
average during 10 months in 1986 (Fig. 9). High 
positive temperature deviations were recorded in 
Febuary and November (3.3 and 3.2 C, 
respectively). Annual flow of the Pee Dee River 
was 50 percent of normal (Dowgiallo et al., 1987). 


Spat settlement 


Spat settlement started in late April or mid 
May and ended in late October or mid November 
during each of the four years. Settlement was 
initiated during a narrow range of water 
temperatures (21.6 - 23.2 C) and was continuous 
with two or more peaks during most summers. The 
first peak typically occured in early June and the 
second occured in late July or early August. The 
second peak was usually more intense and generally 
accounted for 30 percent of annual set. Spat 
settlement was significantly higher in 1984 
(p<0.05) at both Town Creek and Old Man Creek 
(Fig. 10). Oyster Landing also exhibited highest 
settlement in 1984 although total settlement was 
not significantly higher than that observed in 
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Figure 8. (A) Deviation in 1986 monthly 
precipitation values from 30 year average. 

(B) Monthly precipitation during 1986 (shaded) 
and 30 year average (open). 
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Figure 9. (A) Deviation in 1986 monthly 
temperature values from 30 year average. (B) 
Monthly temperature during 1986 (shaded) and 30 
year average (open). 


1983. Settlement was lowest during 1986 at all 
three sites and this pattern could be related to 
the absence of a late summer pulse in late July 
and early August. The presence of high densities 
of late stage larvae in the water column in late 
July indicated reproductive success, but a 
corresponding settlement peak did not follow. 


July and early August water temperatures were 
substantially higher in 1986 than values recorded 
for the previous three years (Fig. 11). Also, 
salinities were high and unusually stable during 
this period. The combined effect of high 
temperatures and high stable salinities was likely 


responsible for the absence of a late summer 
settlement peak. The interaction of high 
temperature and high salinity has been 
demonstrated to cause high mortality in Ostrea 
edulis larvae (Robert, et al., 1988). In 
addition, high intertidal air temperatures and 
dessication stress may have adversely affected 
spat survivorship in the intertidal zone. 
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Figure 10. Mean total annual spat settlement per 
harness (Oyster Landing, OL; Old Man Creek, OM; 
Town Creek, TC). 
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Figure 11. Mean daily temperatures and salinities 
observed at Glambank Creek, Town Creek, and Oyster 
Landing during summer 1986. Bars represent 
monthly means + 2 S.D. for the period 1983-1985. 


The area from Charleston SC to Beaufort SC 
experienced heavy oyster mortalities in 1986 as a 
result of Dermo disease (Perkinsus marinus). The 
dieoffs followed drought and high temperatures and 
may have been related to reduced nutrients, 
temperature and salinity stress enhancing oyster 
susceptibility to disease (Dowgiallo et al., 
1987). 


Fish Community Structure 


Biweekly collections of the more than 40 
species of fishes, penaeid shrimps, and crabs 
which utilize high marsh creeks as "nursery" 
habitat have indicated that community composition 
and population dynamics are influenced by 
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meteorological processes on both short and long 
time scales. Distinct seasonal patterns of 
occurrence for the dominant species have been 
documented and the relative abundance of the 
numerically dominant species is predictable from 
year to year. With the exception of a few species 
(e. g. mummichog), most species recruited to the 
shallow creeks from remote spawning locations, 
including deep ocean areas. The timing of first 
occurrence for many fishes and shrimps was usually 
within two or three weeks of the five year average 
for that species and the succession of species 
from January through August was predictable. 
However, the size (length) distribution at the 
time of first occurrence for spot, white mullet, 
and other dominant species was variable among 
years. Apparent growth rates and the maximum size 
of individuals collected in the high marsh creek 
were generally similar among years. Delays in the 
arrival time of young of the year and/or later 
classes of spot and white mullet appeared to 
correspond with interannual variations in both 
water temperature and salinity. 


Patterns of abundance within each of four years 
indicated that spot were the most abundant fish 
each year and that spot abundance varied 
considerably within and among years (Fig. 12). 
Spot abundance was highest in 1987 following an 
unusually long period of low salinity, a condition 
which is known to enhance larval spot recruitment 
(Allen and Barker, 1988) and was lowest in 1986 
when unusually high temperatures and salinities 
occurred (Fig. 12). High mummichog and silver 
perch abundance also reflect the drought 
conditions during the summer of 1986. During 1984 
and 1985, large mid-summer fluctuations in both 
salinity (one increase and one decrease) and water 
temperature apparently resulted in major 
reductions in spot abundance. Short-term changes 
in weather patterns have also been observed to 
influence the catch composition during the summer. 
For instance, significant rainfall during the days 
preceding the biweekly sampling reduced salinity 
in the tidal creek pool and lower abundance or 
absence of some species was usually observed. 


More analyses of the abundance, length 
distribution, weight, and species composition data 
sets will be conducted to quantify the climatic 
and short-term meteorological influences, but 
preliminary interpretations suggest that much of 
the long and short-term variability can be 
explained by precipitation and temperature 
patterns. Fluctuations in climate which affect the 
patterns of "nursery" habitat utilization may have 
significant implications for the success 
(abundance and size) of many economically 
important fisheries in southeastern estuaries. 


Bay Intrusion 


Nitrogen concentrations in North Inlet average 
3 to 10 percent of the concentration present in 
Winyah Bay. In the North Inlet estuary, the major 
nitrogen form is dissolved organic nitrogen (60 
percent, 20.25 ug At./1l) with particulate nitrogen 
(34 percent, 11.52 ug At./1) also making up a 
substantial fraction. It is estimated that 
approximately 10 percent of the particulate 
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Figure 12. Oyster Landing seine composition for 
the period 1984-1987. 


nitrogen is living biomass (Erkenbreker and 
Stevenson, 1978). Ammonia comprises 76 percent 
(1.73 ug At./1) of the inorganic nitrogen, but 
only 5 percent of the total nitrogen present. 
Nitrate nitrite is less than 1 percent (0.55 ug 
At./1) of the overall nitrogen available in North 
Inlet. In Winyah Bay, inorganic nitrogen (nitrate, 
nitrite and ammonia) comprises 28 percent of the 
total nitrogen (106.42 ug At./1l) present (DHEC, 
1988). Nitrate/nitrite concentrations (16.57 ug 
At./1) are slightly higher than ammonia (14.07 ug 
At./l). The remaining 72 percent of the nitrogen 
is dissolved organic. No studies have focused on 
particulate or gaseous forms of nitrogen in Winyah 
Bay. Seasonal nitrogen patterns in North Inlet are 
strongly related to the temperature cycle while 
variation in Winyah Bay is moderated primarily by 
freshwater discharge. 


Nitrate-nitrite nitrogen concentrations in 
North Inlet are inversely related to salinity 
(R*= 0.45) suggesting external sources for this 








nutrient. Two sampling stations affected by 
freshwater runoff are Oyster Landing (runoff from 
undeveloped forest system) and Clambank Creek 
(intrusion from Winyah Bay). Examination of 
nitrate-nitrite nitrogen concentrations measured 
daily at Clambank Creek and Oyster Landing during 
the period September 1978 through May 1986 
indicated that average concentrations are higher 
at Clambank Creek than Oyster Landing (0.757 ug 
At./1 and 0.559 ug At./1, respectively; N=1577). 
Higher nitrate-nitrite nitrogen concentrations in 
Winyah Bay relative to North Inlet (Fig. 13) and 
intensive studies of nutrient exchange at the 
creeks connecting the two systems (No Man’s Friend 
Creek and South Jones Creek) suggested that Winyah 
Bay may be an important source of nutrients for 
North Inlet (Michener and Allen, 1982). Analysis 
of variance (Table 1) indicated that wind 
velocity, wind direction, maximum daily river 
height and the interaction among the three 
components were related to nitrate-nitrite 
nitrogen concentrations at Clambank Creek. 
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Figure 13. Nitrate-nitrite mean monthly 
concentrations at Winyah Bay and North Inlet. 


Table 1. Analysis of variance for Log-transformed 
nitrate-nitrite nitrogen concentrations during 
1983 at Clambank Creek. 


Source of Variation at SS F Pr>F 


Wind velocity (WV) 1 0.244 5.44 0.020 
Wind direction (WD) 3 0.314 2.38 0.071 
Rain A. 0.114 2.59 0.109 
River height (RH) 1 0.837 19.02 0.0001 
WV*WD 3 0.042 0.31 0.815 
WV*WD*RH 4 0.433 2.46 0.046 
Error 203 8.927 

Total 216 11.233 


Full Model: Pr>F=0.0001, R%=0.205 


Analysis of the 15 low salinity events (112 
days total) that occurred during 1983 revealed 
that southwesterly winds were most frequently 
associated with the initiation of the intrusion 
event (7 of 15). Intrusion duration, minimum 
salinity, and maximum nutrient concentrations were 
related to predominant wind direction during 
events, wind velocity, and maximum river height. 
For example, one 2-day event (Jan 24) was 
associated with northwesterly winds and below 


average maximum daily river height on day 1. The 
lowest salinity (28 ppt) was recorded on the first 
day and increased on day 2 (30 ppt). Maximum 
observed nitrate-nitrite nitrogen concentration 
was 0.40 ug At./l. In contrast, during one 14-day 
event (March 3), salinities dropped to 6 ppt and 
nitrate-nitrite nitrogen concentrations reached 
5.01 ug At./l at Clambank Creek. Maximum river 
height was above average on day 1. Southwesterly 
winds accompanied the initiation of the event and 
were dominant throughout the event (12 of 14 
days). During this event, maximum nitrate-nitrite 
nitrogen concentrations reached 1.87 ug At./1l and 
salinities dropped to 2.0 ppt at Oyster Landing. 
The sharp decline in salinity at Oyster Landing 
was related to a single 1.8 cm rain storm and 
freshwater runoff from the surrounding watershed. 
Higher nitrate-nitrite nitrogen concentrations at 
Clambank Creek are likely related to the proximity 
of this site to Winyah Bay. Nutrient inputs 
occurring from intrusion events undergo mixing and 
are diluted considerably prior to reaching more 
distant portions of the North Inlet Estuary (e.g. 
Oyster Landing). 


CONCLUSIONS 


Spatial variability in rainfall and seasonal 
changes in wind velocity and directionality are 
characteristic features of the climate at the 
North Inlet Estuary. Impacts of hurricanes, 
tropical storms, and droughts are related to the 
recurrence frequency, intensity, and aerial extent 
of the event. Precipitation resulting from 
hurricanes or tropical storms may be considered an 
integral component of the regional water balance. 
Direct landfall of a Glass 4 or 5 hurricane would 
likely result in noticeable changes to the North 
Inlet estuarine system. Tidal circulation patterns 
and sedimentation rates are predicted to be 
affected appreciably. Droughts affect the North 
Inlet Estuary less frequently than tropical storms 
or hurricanes (recurrence frequencies of 5.9 and 
2.6 years, respectively). Droughts occurring for 
extended periods over broad spatial scales can 
have significant impacts on productivity within 
the North Inlet Estuary. Low recurrence rates of 
severe droughts coupled with variable salinity and 
temperature tolerances of flora and fauna may have 
implications for community dynamics as well as 
reproductive success and recruitment. Intrusion of 
Winyah Bay waters into the North Inlet Estuary 
occurs frequently (i.e. 30 percent of observations 
during 1983) and may be a dominant factor 
affecting seasonal productivity patterns. Enhanced 
sedimentation, closure of connecting creeks, or 
opening of new connections between the North Inlet 
and Winyah Bay Estuaries would likely result in 
pronounced alteration of the available nutrient 
pool and subsequent primary and secondary 
production in the North Inlet Estuary. 


Prior to examination of the potential effects 
of global climatic change (i.e. enhanced CO, 
levels, sea level rise, etc.) on regional and 
local biomes, it is necessary to understand 
ecosystem responses to naturally occurring chronic 
and acute climatic events. In this paper, we have 
examined the natural variability associated with 
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several meteorological parameters and, where 
possible, related ecosystem responses to 
individual climatic events. Further analysis of 
the database and additional experimental studies 
will provide us with a better understanding of the 
relationship between climate and ecosystem 
processes and enable us to more accurately predict 
ecosystem responses to global and regional 
climatic change. Also, with an increased research 
effort directed towards the analysis of the 
interactions between climate, material flows, and 
ecological processes, the North Inlet LTER 
research will bring us closer to an understanding 
of how regional phenomenon alter the long-term 
structure of a landscape. 
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LAKES AS INDICATORS OF AND RESPONDERS TO CLIMATE CHANGEL/ 


Dale M. Robertson2/ 








Abstract--Lakes integrate short-term fluctuations in the 
atmospheric climate. Interannual changes in the response of 
thermal features, namely, epilimnion and hypolimnion 
temperatures, thermocline depth, duration of ice cover and 


freeze and break up dates, 


can therefore be used to detect 


and quantify past changes in the climate which are often 
obscured by short-term variability and secular changes in 


meteorological measurements. 


Given future climatic scenarios 


and the quantitative relationship with air temperatures, 
future lake conditions can be estimated. Extensive ice cover 
and water temperature data for Lake Mendota, Wisconsin, 
beginning in 1856, are used to demonstrate the utility of 
historical lake records to detect past changes in mean air 
temperature. Future lake conditions are forecast in response 
to projected "greenhouse" warming scenarios. 


Keywords: Air temperature, water temperature, historical 
lake records, “greenhouse” warming 





INTRODUCTION 


"Climate" is defined as the collective state 
of the earth's atmosphere for a given place for a 
specified interval of time (Landsberg, 1960). 
Therefore, when one describes climatic changes, 
the time frame and meteorological parameters must 
be specified. Here, climate and climatic change 
will primarily refer to air temperature and 
changes in mean air temperature, respectively. 
Climatic changes have been occurring since the 
earth formed and will continue to occur in the 
future. Time scales of climatic change range 
from changes over millions of years to changes 
over days. To detect historical and future 
climatic changes, the climate or weather 
conditions must have been recorded in a manner 
allowing a rigorous comparison of climatic 
conditions through time. Many indicators have 
been used to detect past climatic changes; 
however, most only detect changes over specific 
time scales (Heicht and others, 1979). 
Indicators used to detect changes over time 
scales from 10's to 100's years include 
meteorological (Budyko, 1977; Wahl and Lawson, 
1970), phenological (Novitzkie and Devaul, 1978), 
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tree ring (Fritts and others, 1971; Meko and 
others, 1985), and lake sediment data (Webb and 
Bryson, 1972; Bernabo, 1981; Swain, 1978). In 
addition to these, ice cover and water 
temperature records of lakes also serve as 
climatic indicators (McFadden, 1965; Tramoni, 
1985; Palecki and Barry, 1986; Strub and others, 
1985; Robertson, 1989). 


Recently, major questions dealing with 
climate change have focused on anthropogenic 
effects, such as: Have human disturbances to the 
atmosphere and land affected historical climates 
and will they affect future climates? Will 
increasing concentrations of atmospheric CO> and 
other “greenhouse" gases cause surface 
temperatures of the earth to increase? If air 
temperature increases, how will this affect the 
various ecosystems on the earth, specifically 
aquatic ecosystems? In this paper, I demonstrate 
how lake ice cover and water temperature records 
for a specific lake, Lake Mendota, WI, can be 
used to detect and quantify past local changes in 
mean air temperature that otherwise be 
undetectable using past meteorological records. 
Future mean conditions of the lake are forecast 
using projected climatic warming scenarios and 
the calibrated response of the lake parameters. 


Changes in mean air temperature over time 
scales ranging from 10's to 100's of years are 
often difficult to discern using biological 
indicators because the effects of air temperature 
are usually integrated with effects of other 
meteorological conditions, primarily 
precipitation and humidity (Swain, 1978). 


Biological indicators also have variable response 
times dependent on the life span or growth rate 
of the organism. Air temperature records 
themselves often have secular changes or 
systematic biases incorporated into the records 
which may be larger than subtle climatic changes. 
These systematic biases can be caused by changes 
in the equipment, station location, station 
design, and time of observation (Schaal and Dale, 
1977; Karl and Williams, 1987). Some systematic 
biases can be removed using coinciding data from 
nearby stations (Karl and Williams, 1987) or 
using theoretically based corrections to adjust 
for some known changes in observational 
techniques. However, changes in early 
observational techniques, which are usually not 
documented, can be almost impossible to remove 
resulting in early air temperature records being 
of questionable quality and often very 
misleading. 


Ice cover (freeze and break up dates, and 
total days of ice cover) and water temperature 
records for a lake represent an integration of 
local weather conditions. These parameters often 
integrate daily fluctuations in weather 
conditions, primarily air temperature, in a 
manner less subject to secular biases than other 
climatic indicators including direct 
meteorological measurements. Ice cover and water 
temperature data are not internally recorded; 
each has to be observed and documented by an 
observer. Therefore, lake records, as well as 
any parameter documented by observers such as 
meteorological data, may have systematic biases 
incorporated into the records caused by changes 
in observational techniques or criteria. 

However, for some lakes, such as Lake Mendota, 
the transition to and from complete ice cover are 
very rapid and can often be resolved to one or 
two days. For the ice cover records of such 
lakes, variability introduced by observer 
subjectivity is small compared to the variability 
of the system. 


Annual differences in ice cover and water 
temperature for small to medium sized lakes are 
primarily driven by changes in air temperature 
for specific seasons (McFadden, 1965; Palecki and 
Barry, 1986; Robertson, 1989). Freeze dates 
represent an integration of fall air 
temperatures; break up dates represent the period 
between freezing and break up (winter and early 
spring) and total ice duration represents the 
entire fall, winter, and early spring period. 

The length of the integration period represented 
by the ice records depend on lake depth, with 
deeper lakes integrating a longer time period 
(McFadden, 1965). Water temperatures represent 
an integration of weather conditions from the 
spring, summer, and fall periods. 


Lake records have had limited use as climatic 
indicators because long-term continuous ice cover 
and water temperature measurements are extremely 
rare. Few lakes have water temperature records 
that extend more than 20 years or exist prior to 
1920. Possibly the most extensive water 
temperature records available exist for Lake 


Mendota, WI, which begin in 1994 and continue 
intermittently to present. Long-term ice records 
are a little more common, and a few exist from as 
early as 1840 in North America (Tramoni, 1985). 
The Longest continuous records available for any 
lake in North America begin in 1855, also for 
Lake Mendota (Tramoni, 1985). 


LAKE MENDOTA 


Lake Mendota is a dimictic, eutrophic lake 
located in Madison in south-central Wisconsin 
(43°4'N, 89°24'W). The lake has a surface area 
of 3,940 ha, a maximum depth of 26 m, 4 mean 
depth of 12.4 m, and a maximum fetch of 9.8 km 
(Wisconsin Department of Natural Resources, 
1961). 


DETECTING CLIMATIC CHANGES 


Detecting prior climatic changes from lake 
records or other climatic indicators is generally 
approached by examining a time series for 
statistically significant changes in mean 
conditions. Ice records of Lake Mendota are used 
to demonstrate how lake indicators can be used to 
detect climatic change. In Figure i, mean ice 
duration (as demonstrated by the seven year 
moving average) of Lake Mendota has significantly 
changed since 1856. In fact, most of the 
climatic change appears to have been rather 
abrupt around 1890. Since 1980, there appears to 
possibly be another transition to shorter ice 
cover. 


Climatic changes are generally thought to 
occur over several years. To determine whether a 
significant change has occurred in the indicator, 
the data are tested to determine whether there is 
a significant slope or change through time. 

Often the indicator time series has significant 


LAKE MENDOTA ICE DURATION 
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Figure 1. Total days of winter ice cover (ice 
duration) for Lake Mendota. The seven year 
moving average, plotted in year four, is 
superimposed on the observed data as an estimate 
of mean ice durations. 
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autocorrelation; therefore, ARIMA analyses are 
used to remove the autocorrelation, prior to 
testing for significant changes. There is a 
significant downward slope or change to shorter 
ice durations for Lake Mendota since 1856 (a< 
0.05). However, when examining the time series 
of the indicator, often the change appears to 
occur rather abruptly rather than gradually. One 
technique which has been used to test whether a 
rather abrupt change better describes the change 
in the indicator is ARIMA intervention analyses 
(Box and Tiao, 1975). This technique tests for 
significant slopes or changes before and after a 
possible abrupt change. This technique can be 
used with or without significant autocorrelation 
in the dataset. The timing of the change is 
defined by the user. The “best" model is often 
selected by comparing t-ratios of consecutive 
yearly interventions. The “best" model 
describing the change in the mean ice duration of 
Lake Mendota had a constant mean of 118 days from 
1856 to 1888 (Figure 2A). In 1888 a rather 
abrupt change occurred, modeled as a step change 
followed by a stationary mean of 103 days. The 
decrease in ice duration of 15 days is indicative 
of higher mean air temperatures after 1888. 
Whether the change occurred precisely in 1888 or 
over several years only slightly modifies 
estimated changes in mean conditions. 


Often when examining a time series, a period 
of unusual years appears present. For example, 
the ice duration from 1980 to 1988 appear to be 
shorter than that prior to 1980 (Figure 1). The 
seven year moving average ice duration, from 1982 
to 1988, is the shortest on record, even if the 
record short ice duration of 1983 was replaced by 
the mean ice duration from 1889 to 1979. One 
technique to test whether the mean of a specific 
period is significantly different from that of 
another is a Student t test (Burr, 1974). The 
period from 1980 to 1988 had significantly 
shorter ice duration, by 12 days, than the prior 
period (a < 0.05). Therefore, the “best" 
description of the mean ice duration of Lake 
Mendota is demonstrated in Figure 2A. 


Similar procedures to describe mean freeze 
and breakup dates were used and the results are 
demonstrated in Figure 2B, C, respectively. A 
change in mean freeze date occurred in about 
1888, but no significant change was observed in 
recent years. From the freeze date data, it was 
unable to be determined whether the first, 
cooler, period ended abruptly or whether it was a 
gradual transition period. For the purpose of 
further analyses, it was assumed the first period 
ended abruptly in 1888. This gave a change in 
mean freeze date from Julian day 348 to 356, a 
change of 8 days based on mean freeze dates. 
Three climatic periods, similar to those found 
with ice duration records, were found from break 
up dates. The first period lasted from the 
beginning of record to 1888, when mean break up 
date changed rather abruptly from Julian day 101 
to 94, a 7 day change. The second period lasted 
from 1889 to approximately 1979. The later mean 
break up date of the first period indicated a 
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cooler climate prior to 1888. No significant 
trends occurred in the first or second periods. 
Since 1980, mean break up date appears to be in 
transition to earlier dates. Presently, the 
moving seven year average, from 1982 to 1988, has 
obtained the earliest on record at Julian day 85. 
The mean break up date from 1980 to 1988 was 
Julian day 86, an 8 day earlier mean break up 
date than that of the second period, indicative 
of a warmer period. 
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Figure 2. Observed annual ice cover records for 
Lake Mendota with the "best" estimate of how the 
mean for each parameter has changed. The 
idealized means were determined from several 
ARIMA intervention models fit to the observed 
data. (A) Ice Duration. (B) Freeze Date. (C) 
Break Up Date. 











QUANTIFYING CLIMATIC CHANGES 


To quantify these changes in terms of 
specific meteorological parameters, we must 
determine what that change means in terms of air 
temperatures. This is done by determining how 
the indicator has responded to changes in that 
meteorological parameter, i.e., calibrate the 
indicator. The calibration can often be very 
complex because of the integrating nature of the 
indicators. 


Two approaches have been used to determine 
how climatic indicators respond to changes ina 
specific meteorological parameter: the dynamic, 
or mechanistic approach, and the empirical or 
statistical approach (Wright, 1980). Since the 
calibration techniques are much different for ice 
cover and water temperature, both are described 
below. 


In the dynamical approach, mathematical 
models are developed to simulate the response of 
specific indicators to weather or climatic 
conditions, based on equations simulating the 
physical processes occurring in the system. All 
important parameters must be included in the 
model to properly simulate the system. The 
models are calibrated using concurrent historical 
records for the climatic indicator and measured 
meteorological data. The response of the 
indicator to a specific meteorological parameter, 
such as air temperature, is then determined by 
running the model for various simulations with 
only the specified meteorological parameter 
altered. The change in response of the indicator 
in these simulations is then used to determine 
the relationship between the indicator and the 
specified meteorological parameter. This 
approach is more commonly used than the empirical 
approach because only a few years of measured 
data are required. 


In the empirical approach, we exploit the 
fact that nature has been performing experiments 
with the climatic system for a long time; the 
results of those experiments are our 
observational records (Wright, 1980). These 
records are used to develop statistical 
relationships between the important 
meteorological parameters, including air 
temperature, and the climatic indicator. If 
statistically significant relationships involving 
air temperature are found, the equation(s) are 
then used to determine the pattern and response 
of the indicator to changes in air temperature. 


Whenever possible, dynamic modeis should be 
used in concert with empirical models to 
determine if the results agree. A similar 
quantitative response using two independent 
approaches adds more confidence to the results 
than results found only using one technique. 


Ice Cover 
Freeze and break up dates were analyzed 


separately to quantitatively determine the 
magnitude and seasonal timing of the climatic 


changes rather than calibrating ice duration as a 
total winter indicator. Two techniques have been 
used to develop lake ice cover records into 
temperature indices: fixed period regression 
analyses (Rannie, 1983; Tramoni and others, 1985; 
Palecki and Barry, 1986; Robertson, 1989) and 
variable length air temperature integration 
(McFadden, 1965; Robertson, 1989). 


Regression analyses have been used to 
calibrate ice records of specific lakes by 
regressing specific air temperature summary 
parameters, such as a monthly or bi-monthly 
average air temperature, with observed freeze or 
break up dates. To calibrate the freeze date for 
Lake Mendota, annual freeze dates from 1884 to 
1988 were regressed against several combinations 
of monthly mean air temperatures. The earlier 
years in the record were not included because of 
secular biases in meteorological records prior to 
1884. Average November to December air 
temperatures were most strongly correlated with 
observed freeze dates, R2 = 0.61 (Figure 3A). 
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Figure 3. Calibration of mean freeze date into 
an index of fall and early winter air 
temperatures. (A) Fixed period analysis. 
Observed freeze dates from 1884 to 1988 versus 
the deviation from the long-term average November 
to December air temperature (1884 to 1988). (B) 
Moving integrated air temperature analysis 
("Freeze™ model). 
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The slope of the regression suggests mean freeze 
date changes 4.32 days (+ 0.6 days, 95% C.I.) per 
1.0°C change in average November to December air 
temperature. Deviation from the average November 
to December air temperature, from 1864 to 1988, 
was used as the x-axis to simplify comparison of 
model results. Interpolation of the regression 
slope suggests the 8 day earlier freeze date 
before 1888 corresponds to average November to 
December air temperatures which were 1.7°C lower 
than those after 1888. 


The second technique to calibrate ice records 
involves using a moving index which accumulates 
through time, such as a moving mean air 
temperature (McFadden, 1965) or integrated air 
temperature (Robertson, 1989). When the moving 
index reaches a specified threshold value, the 
specific event such as freezing or break up is 
estimated to occur. When calculating simple 
moving average air temperatures, days early in 
the period of the moving average are weighted 
similar to days just prior to freezing or break 
up; whereas, integrated air temperatures can 
weight the more recent air temperatures more 
heavily. Therefore, integrated air temperatures 
were used. Integrated moving air temperatures 
were used to predict daily water temperatures and 
ultimately when the surface temperature became 
0°c and the lake froze. This was done by 
weighting the prior day's predicted water 
temperature by a specific value and adding to 
this the change estimated for that day, which is 
dependent on the difference in air and water 
temperatures (Eq. 1). Therefore, the most recent 
days are more influential to the present daily 
water temperature value. Eq. 1 can be adjusted 
to add the effects of daily synoptic conditions 
such as strong winds which could cause the lake 
to freeze later or break up earlier than expected 
without modifying the estimated day's water 
temperature. 


Wr a Wry + a + K * (Any —s Wp_1) (1) 
or 

Wp = a+ b* Ap_y +c * Wpiy (2) 
where: 


Wp = Water temperature on Day T 

Wp_y = Water temperature on Day T - 1 
Ap.1 = Air temperature (°C) on Day T - 1 
a, b, c, and K are Empirical Constants: 


a= 0.056; b = 0.946; c = 0.326 


Initially, Wr-1 = 28.70 on 1 October 


Integrated moving air temperatures were 
developed into a "Freeze" model (Eq. 2) by 
calibrating the coefficients such that the 
estimated surface water temperature (Wr), 
calculated starting on 1 October of each year, 
became < 0°C very near to the day the lake was 
observed to freeze. The “best" coefficients for 
estimating annual freeze dates for Lake Mendota 
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were determined by iteratively calculating the 
coefficients which minimized the sum of square 
errors for a 40 year reference period, 1948 to 
1987. This minimization process was done using 
the Nelder-Mead simplex procedure (O'Neill, 
1971). Additional terms were added to Eq. 2 in 
an attempt to adjust for synoptic conditions. 
These terms had the effect of raising or lowering 
the threshold criteria but did not affect future 
predictions of Wp, i.e., these terms did not 
alter the following Wp-1- Therefore, with a 
strong wind lower water temperatures may be 
required for the lake to freeze. Additional 
terms, however, did not significantly improve the 
model. 


For annual models, such as the "Freeze" model 
or a similar "Break up" model, the estimated 
effects of changes in mean air temperature are 
made by comparing average estimated freeze or 
break up dates of a reference period, a 30 year 
period from 1948 to 1977 was used here, with 
average predicted dates obtained from the models 
after the daily air temperatures for those years 
are raised or lowered by the simulated change in 
air temperature. Daily air temperatures for all 
30 years were lowered 1, 2, and 3°C and raised 1, 
2, 3, 4, and 5°C. The average freeze dates 
estimated by the "Freeze" model were then 
regressed on the simulated change in air 
temperature (Figure 3B). The slope of the 
regression suggests a change in freeze date of 
4.83 days per 1.0°C change in air temperature, 
very similar to the 4.32 days per 1°C change 
estimated using the fixed period analysis. 
Interpolation of the regression slope suggests 
the 8 day earlier freeze date before 1888 
corresponds to average November to December air 
temperatures which were 1.5°C lower than those 
after 1888. This approach to estimate the 
effects of changes in air temperature assumed 
that the daily variability in weather conditions 
for the simulated changes in air temperature was 
and would be similar to that of the reference 
period. 


A similar analysis for break up dates 
suggests a change of approximately 6 days per 1°C 
change in mid winter to early spring, January to 
April, air temperatures. Therefore, the changes 
in mean break up dates (Figure 2C) suggest the 
period prior to 1889 was 1.1°C lower than the 
second period and the period after 1979 was 1.2°C 
higher than the second period. 


Water Temperatures 


The DYnamical REservior Simulation Model 
(DYRESM) (Imberger and Patterson, 1981) was used 
to simulate the physical processes and ultimately 
water temperatures in Lake Mendota responding to 
the forcing meteorological conditions. DYRESM is 
a one-dimensional numerical model used to predict 
water temperature and salinity profiles of small 
to medium sized lakes and reservoirs given 
initial temperature and salinity profiles, lake 
morphometry, meteorological conditions, and 
inflow and outflow information. For each 
simulation DYRESM was started on the day after 
ice break up when the lake has been shown to be 


very near 3°C (Neese and Bunge, 1956; Robertson, 
1989). DYRESM was calibrated for Lake Mendota 
using two years of daily weather data and 
approximately biweekly measured water temperature 
profiles from shortly after ice break up until a 
few weeks prior to freeze up for 1958 and 1966 
(Robertson, 1989). These years were selected 
because the thermal structures varied 
considerably between years. The “best" values 
for the empirical coefficients used in DYRESM 
were determined by iteratively calculating the 
values which minimize the overall average of the 
mean daily absolute errors for both years (Eq. 
3). The minimization process was done using the 
Nelder-Mead simplex procedure (O'Neill, 1971). 
The calibrated model simulation for 1958 is 
compared with observed data in Figure 4. DYRESM 
simulated the daily temperature progression very 
well, with the possible exception during late 


fall. 
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MDAE = Mean Daily Absolute Error 
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= DYRESM simulated water 
temperature 


Tmodel 


610 1418 22 26 





HONG G09 Goonooog dG 


10 
15 
20 
o 
‘i 
254— 


DYRESM model verification for Lake 


Figure 4. 
Mendota in 1958. Computed vertical temperature 
profiles (solid connected dots) and observed 


(open squares) temperatures. Profiles were 
measured and compared approximately biweekly from 


To determine how changes in air temperature 
influence the thermal structure of the lake, 
DYRESM was run for 11 years for each of 6 air 
temperature regimes. Eleven years of daily 
meteorological data (1958-63, 1966-67, 1971-73) 
were used as the base or reference air 
temperature regime. Five other air temperature 
regimes were obtained by adding -3, ~1, +1, +3, 
and +5°C to each daily average air temperature 
for all ll years. Since annual ice break up 
dates or the date at which the simulation began 
should change with air temperature, forecasted 
ice out dates were determined for each year based 
on the "Break Up" model, using daily adjusted air 
temperatures as input for the model. 


For each annual simulation, morphometrically 
weighted monthly mean epilimnion or shallow layer 
temperature (0 to 7m), hypolimnion or deep layer 
temperature (15 to 18m), and average monthly 
thermocline depth were computed. Average monthly 
values for these three parameters were computed 
for each air temperature regime and subtracted 
from that for the the reference regime. 
Differences from the means of reference period 
represent the average effect caused solely by the 
change in air temperature. 


Changes in average monthly epilimnion 
temperatures were related directly to changes in 
air temperature, with water temperatures 
increasing with air temperature increases and 
decreasing with air temperature decreases (Figure 
5A). However, the changes were smaller than the 
air temperature change and varied seasonally, 
with the largest response occurring in spring. 
The response of epilimnion temperatures varied 
from approximately 65% of the change in air 
temperature in May to approximately 40% for most 
summer months. 
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Julian day number 129 (8 May) to 325 (20 
November) and are in order of day number from 


left to right. Each date is offset by 6°C. 
date has the same relative x scale, indicated by 
the scales on the top of the figure for the Ist 

and 8th profiles. 


Each 
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Changes in average monthly hypolimnion 
temperatures were also related directly to 
changes in air temperature, with mean water 
temperatures increasing in response to air 
temperature increases and decreasing in response 
to air temperature decreases (Figure 5B ). The 
response was again smaller than the change in air 
temperature and varied seasonally from 
approximately 45% of the change in air 
temperature in May to 25 to 30% for most of the 
summer. 


Simulated thermocline depths displayed little 
response to changes in air temperature (Figure 
5C). There was a slight indication of shallower 
thermoclines from August to October with 
increasing air temperature; however, these 
changes were not statistically significant. 


PROJECTED RESPONSE TO FUTURE CLIMATIC CHANGE 


The future response of the thermal structure 
and ice phenology of Lake Mendota can be made 
based on climatic warming scenarios estimated 
from various General Circulation Models (Hansen 
and others, 1984; Manabe and Wetherald, 1986; 
Schlesinger and Zhao, 1988) in response to 
increased atmospheric CO» and other "greenhouse" 
gases in concert with the quantitative 
relationships found above. If all other non-air 
temperature parameters remain unchanged, a 
uniform monthly increase in air temperature of 
1°C should cause the responses listed in Table 1 
for Lake Mendota. With further increases in air 
temperature these responses should be 
approximately additive, i.e., a 2°C increase 
should cause the responses to double, etc. 
However, with a 4 and 5°C increase in air 
temperatures, results suggest the lake will often 
freeze and break up more than one time per winter 
which in turn will further shorten total ice 
duration (Figure 6). The frequency of refreezing 
should increase with continued warming. With a 
5°C increase, expected to occur around 2050 to 


Table 1. Seasonal response of Lake Mendota to a 1°C increase in air 
temperature, all other meteorological parameters were assumed 
to remain unchanged. 





Season Response 
Spring - Entire water column temperatures higher by 0.65°C. 
Early Summer _ - Higher epilimnion temperatures by 0.5 to 0.69C. 


- Earlier stratification. 
- Higher hypolimnion temperatures by 0.3°C. 


Mid - Late - Higher epilimnion temperatures by 0.4 to 0.59C. 
Summer - Little change in thermocline depths. 
- Higher hypolimnion temperatures by 0.3°C. 
Fall - Higher epilimnion temperatures by 0.4 to 0.5°C, 


~ Longer stratification. 

- Higher hypolimnion temperatures by 0.4 to 0.5°C. 
- Mean freeze date delayed by 5 days. 

- Mean break up date earlier by 6 days. 

- Mean total ice duration shorter by 11 days. 


Early Winter 
Late Winter 
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Figure 5. Estimated changes in the thermal 
structure of Lake Mendota in response to changes 
in air temperature. The simulated change in air 
temperature is labelled on each curve. All 
changes are differences from the mean monthly or 
seasonal values. All monthly and seasonal 


changes have the 95% confidence limits. (A) 
Change in mean epilimnion temperture. (B) Change 
in mean hypolimnion temperature. (C) Change in 


mean thermocline depth. Mean epilimnion and 
hypolimnion temperatures are defined as the 
morphometrically weighted mean water temperature 
from 0 to 7m and 15 to 18m, respectively. 
Thermocline depth is defined as the depth of the 
maximum temperature gradient. 
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Figure 6. Expected future response of mean ice 
cover to climatic warming as projected by the 
Transient General Circulation Model from the 
Goddard Institute of Space Studies (Hansen and 
others, 1984). Projections were made using the 
"Freeze" and "Break up" models. Historical 
annual records and the best description of 
historical means are included (similar to Figure 
2). 


2060, mean ice duration should be reduced by 65 
days and results of the "Freeze" model suggest 
the lake should remain ice free all winter with a 
frequency of 1 out of 30 years. With increasing 
air temperatures the possibility of the lake not 
freezing in a particular winter will increase 
further. 


CONCLUSIONS 


Changes in the response of the thermal 
features of lakes, namely, epilimnion and 
hypolimnion temperatures, thermocline depth, 
freeze and break up dates, and total ice 
duration, are good quantitative indicators of 
past changes in air temperature. These 
parameters integrate the daily fluctuations in 
weather conditions, primarily air temperature, in 
a manner less subject to secular biases than 
other climatic indicators including 
meteorological measurements. Ice cover records 
are especially good climatic indicators because 
each ice parameter has a large response to small 
changes in air temperature. 


Other studies have estimated future climatic 
warming in response to increased levels of 
atmospheric CO and other “greenhouse"™ gases 
using General Circulation Models. Higher air 
temperatures should cause the changes in water 
temperature and ice cover described for a 1.0°C 
increase to be multiplied by a factor equal to 
the increase in air temperature. With an 
increase of 4 to 5°C, winters in which the lake 
freezes and breaks up more than once should 
become more common and further shorten ice 
duration. With an increase of 5°C, mean ice 
duration should be shortened by 65 days and the 
lake should remain ice free with a frequency of 
about 1 in 30 years. 


Because lake ice cover has a large 
response to changes in air temperature and 
because lake ice cover integrates air 
temperature over the seasons when climatic 
warming is projected to be largest (winter and 
spring), a lake's ice cover may be one of the 
best early indicators of climatic warming 
induced by the "Greenhouse Effect." 
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AVAILABILITY AND QUALITY OF EARLY WEATHER DATA, AND 
IMPLICATIONS FOR CLIMATE CHANGE ASSERTTIONS- 
THE ILLINOIS EXAMPLE 


Wayne M. Wendlandt/ 
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Abstract.--Early weather observations from Illinois 


are examined for quality. 


Discontinuous records 


are common, aS are changes in time of observation, 
exposure and method of calculating means. Quality 
control of current data continues to be a problem. 


Keywords: quality climate data, early weather 
observations, Illinois 
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INTRODUCTION 


The first formal and systematic meteorol- 
ogical observations (as opposed to occasional 
records in diaries etc.) in the United States 
were made along the eastern seaboard in the late 
1700s. Only several decades later did such 
observations begin in the central United States. 


The first systematic observations were 
made by committed individuals, like Thomas 
Jefferson, Benjamin Franklin on the east coast, 
and Increase A. Lapham in Wisconsin (Witnah 
1961). The latter’s observations begin in 1827. 
Lapham was particularly interested in beginning a 
network of observing sites, and played a promin- 
ent role in the establishment of a national 
weather service. 


Networks of systematic observations began 
operation under the direction of several or- 
ganizations. The U.S. Army began a_ network 
during the War of 1812. The network was aug- 
mented by private individuals, and by state- 
supported observers in Illinois, Indiana, Iowa, 
Kansas, Michigan, Missouri, Nebraska, Nevada, New 
Jersey and Ohio (Witnah 1961). 





1/state Climatologist, Illinois State Water 
Survey, Champaign IL 61820, and Prof. of 
Geography, University of I1linois-Urbana. 
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Between the 1870s and 1900, funding to the 
Signal Service (which directed the meteorological 
effort) increased ten-fold, and there were more 
than 3,900 reporting weather stations in the U. 
S. by 1904. Table 1 shows the number of current 
station records continuous to the 1840s in 
Illinois by decade. Individual station opera- 
tions may have begun earlier, however no record 


persists. 


The Smithsonian Institution began to 
distribute raingages to individuals by 1853 
(Smithsonian Institution 1855). Meteorological 
data were published by several, and sometimes, 
multiple sources during the mid-1800s, including 
the Smithsonian, and/or some state agency (e.g., 
the Illinois State Weather Service, Chicago 
around the turn of the century), and/or privately 
published sources, e.g., Blodget (1857). Multi- 
ple publishing was common. 


Table 1.-—-Number of [Illinois weather station 
records beginning by decade. 


1840s 1850s 1860s 1870s 1880s 1890s 


1 7 7 9 93 57 
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EARLY OBSERVERS AND OBSERVATIONS IN ILLINOIS 


The accounts of Andrew Ellicott, of Cairo 
IL, refer to the severity of the 1796-1797 winter 
claiming a minimum temperature of -17.2°C on 22 
Dec 1796. 


Table 2.--Early weather observers in Illinois 
after Smithsonian Institution (1855). 
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LOCATION BEGINNING DATE OBSERVER 
Alton May 1849 Johnson 
Athens Jan 1843 Prof. J. Hall 


Dr. S. B. Mead 
Surgeon General S. 
Meacham, S. Brooks, 
I. I. Langguth 
Assistant Surgeon 
Dr. Ryhiner, A. G. 
Bandelier 

T. Dudley & Coffin 
Dr. M. K. Brownson 


Augusta Aug 1833 
Chicago®/ Jul 1832 


Elgin?/ Jan 1824 
Highland Jan 1841 


Jacksonville®/ Apr 1849 
Joliet®/ oct: 1843 


Monroe 1849 Main 

Upper Alton®/ 1849 James 
Warsaw May 1840 Ben Whitaker 
Waukegan®/ 1849 Joslyn 


8/observations continue through the time of this 
writing (though, perhaps discontinuous) 


Many of the earliest observations cont- 
inued only for 1 or 2 years, probably ceasing 
because of lack of continued interest on the part 
of the observer, or because the settlement 
disbanded. Six of the cities in Illinois which 
began observations in the 1870s are no longer 
cities (see Table 2). Clear Lake (Putnam Co.) 
became McNabb in 1900; Prairieville (Lee Co.) is 
now a rural delivery route of Dixon; Irishtown 
(Clinton Co.) remains only as a township name 
today; Brush Hill (Tazewell Co.) is now only the 
name of a rural school near Pekin IL; Jordon’s 
Grove (Randolph Co.) disbanded during the 1880s; 
and Three Mile (Perry Co.), whose population in 
1880 was 27, now only remains as the name of a 
prairie occupying large portions of 2 townships. 


The earliest observations in Illinois are 
from Albion (1856), Alton (1843), Marengo (1855), 
Ottawa (1852), Pekin (1855), White Hall (1854), 
and Woodstock (1857). Two of these, Pekin and 
Woodstock, no longer function as active weather 
stations. 


Published mean and extreme data from 
Illinois first appeared in the 1850s. Blodget 
(1857) published mean temperatures for 7 Illinois 
sites, Chicago (Ft. Dearborn), data beginning in 
1832, Ft. Armstrong (Moline), Ottawa, Augusta, 
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Athens, Ft. Des Moines, and Ft. Dodge in 1854. 
The means were composed of 2 to 12 year’s daily 
data, as opposed to the current use of 30-year 
means. Incidently, there is substantial evidence 
(e.g., see Court 1967-68, Lamb & Changnon 1981) 
that means of shorter periods are more appropri- 
ate when the means are to be specifically used 
for estimating the likely parameter value for the 
next year. Because of these discrepancies, it is 
suggested that the term "normal" not be inter- 
changeably used with mean. 


The Smithsonian Institution was first to 
promote and sponsor systematic meteorological 
observations in the U. S, interestingly, some 
decades before the efforts of the Mannheim 
Meteorological Society in Europe during the 
1890s. The latter made great strides distribut- 
ing instruments and encouraging routine observa- 
tions at sites around the world. In 1873, the 
Smithsonian published (Smithsonian Institution 
1873) mean monthly precipitation values for 40 
Illinois sites, based on 2 to 20 year’s data. In 
1876, mean temperature values for 93 Illinois 
sites, the data covering the same general time as 
those for precipitation (Smithsonian Institution 
1876) appeared. In 1885, mean monthly precipi- 
tation data from 78 sites in Illinois (Smith- 
sonian Institution 1885) were published. 


Fifty-six Illinois weather stations began 
a continuous period of observation prior to 1900 
and continue to the present (listed in Table 3). 
The earliest of these daily records from Illinois 
is from Ottawa. 


Note that the earliest sites (1860s) tend 
to be found in northeastern Illinois, and in the 
southeast along the Wabash River, with another 
score or so beginning operation in the southwest 
part of the state in the 1880s. 


ESTABLISHMENT OF A NATIONAL OBSERVING NETWORK 


The occasional observations made by 
individuals in Illinois and elsewhere prior to 
about 1850 evolved into routine, systematic 
observations first made under the auspices of the 
U.S. Department of War, when Calhoun was Secreta- 
ry, the program being under control of the 
Surgeon General (Maury 1857). The directing 
agency then became the U.S. Department of Agricu- 
iture, and finally the U.S. Department of Com- 
merce. 


One of the first to tout systematic 
observations taken at a common time was Lt. 
Matthew Fontain Maury (USN). In a 1855 letter to 
the Editor of the American Farmer, Maury urged 
farmers to join in systematic weather observing, 
as he had successfully urged sailors earlier 
(Maury 1857). His comments are most perceptive, 
namely, "We want not only corresponding observa- 
tions as to the time, but we want them made with 
instruments that are alike, or that can be com 
pared; with them we may expect to find out 


Table 3.--Illinois weather observing sites for which observations are 
available before 1900 and continue to the present, including the beginning 


year of continuous data. 





STATION NAME BEGINNING DATE 
OF CONTINUOUS DATA 


OTTAWA 1852 
MARENGO 1855 
FLORA 1869 
AURORA 1870 
HAVANA 1870 
ROCKFORD 1872 
SYCAMORE 1881 
GRIGGSVILLE 1882 
MASCOUTAH 1882 
MC LEANSBORO 1882 
PALESTINE 1882 
WHITEHALL 1884 
FATRFIELD 1885 
BELLEVILLE 1887 
CASEY 1887 
DU QUOIN 1887 
GREENVILLE 1887 
HENNEPIN 1887 
HOOPESTON 1887 
JACKSONVILLE 1887 
KANKAKEE 1887 
MT. CARMEL 1887 
OLNEY 1887 
VIRDEN 1887 
SPARTA 1887 
WHEATON 1887 
URBANA 1888 
RUSHVILLE 1889 
CARLINVILLE 1991 
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something certain and valuable concerning the 
movements of this grand and beautiful machine, 
called the atmosphere." (op. cit., p. 488). Not 
only did he recognize the importance of standard 
exposures, but also the importance of a common 
time of observation. 


CHANGES IN OBSERVING TACTICS WITH TIME 


Prior to the mid 1800s, observations from 
official sites were usually made at two or more 
of the following times: 0500, 0600, 0700, 1400, 
1700, 1800, and 2100 local time. Mean and 
extreme daily temperatures could therefore be 
calculated by a variety of methods. Standard 
methods for calculating the daily mean included: 


[T (morning) + T (evening)] / 2 
(Ty Ty 4 oe T57 ] / 3, OC 
ia + T14 + +97, + T21] / 4 


Because of the numerous ways used to 
calculate daily means, the homogeneity of temper- 
ature records reaching back into the 1800s are 
suspect until verified. Even if daily means were 





STATION NAME BEGINNING DATE 
OF CONTINUOUS DATA 
CARBONDALE 1892 
PARIS 1892 
WALNUT 1892 
BLOOMINGTON 1893 
GALVA 1893 
JOLIET 1893 
ANNA 1894 
DECATUR 1894 
MONMOUTH 1894 
TUSCOLA 1894 
WATSEKA 1894 
HILLSBORO 1895 
KASKASKTA 1895 
LAHARPE 1895 
MINONK 1895 
MORRISON 1895 
MT. CARROLL 1895 
MT. VERNON 1895 
PEORIA 1895 
RANTOUL 1895 
TISKTILWA 1895 
CHARLESTON 1896 
DIXON 1896 
HARRISBURG 1898 
PANA 1898 
ALEDO 1899 
VANDALIA 1899 


obtained from maximum/minimum thermometers 
throughout the time of record, changes in the 
time of observation must be known to accommodate 
the time of observation bias (Karl and others 
1986) since these biases can be as much as 2-3°C. 


Aside from station relocations and changes 
in the time of observation, mean monthly tempera- 
ture data from prior to about the turn of the 
century to after are likely not homogeneous. 
Therefore the use of means to determine change or 
trend over time may be improper. In addition, 
only after the mid-1880s in the United States was 
instrument exposure emphasized. At this time 
Hazen and Aitken strongly argued for proper ther- 
mometer exposures, including a radiation screen 
with ventilation. Prior to that time, ther- 
mometers were generally located in an unheated, 
north-facing room, or hung outside in a (hope- 
fully) shaded location. 


OBSERVATIONS AT URBANA 

Beginning in 1873, maximm and minimum 
temperature, 24-hour precipitation, snowfall, 
wind direction, cloud amount, and "intensity" of 
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storms were recorded by Prof. Burrill2/, Director 
of the University of Illinois Experiment Station. 
Perusal of these observations shows that January 
and February of 1875 were very cold indeed, e.g., 
January*‘s mean temperature was 14.0°F (current 
mean: 24.7°F), a value only exceeded 3 times 
since 1888 (1918, 1977 and 1979). Similarly the 
February 1875 mean temperature was 15.0°F (cur- 
rent mean: 29.3°F), coldest of any February of 
record. 


Regular, systematic observations on the 
University of Illinois campus began primarily 
because of 2 significant occurrences, (1) passage 
of the Hatch Act in March 1887 (which authorized 
federal support to establish agricultural experi- 
ment stations, and (2) the foresight of the 
University Trustees, who 6 days after law pas- 
sage, decided that the University should es- 
tablish a weather station. 


Changnon and Boyd (1963) report that the 
Urbana site was first located near the southeast 
corner of the current Lincoln Hall. In June 
1897, the observing site was moved about 200m to 
the south southeast, adjacent to Morrow Plots. 
The site was moved twice more at the Plots, but 
never more than 20m from the second location (in 
April, 1904 and again in June 1948). 


The initial raingage was essentially the 
same as one still used today, i.e., a 20.3 cm (8 
inch) diameter can. The max/min thermometers 
were mounted in a louvered wooden shelter. Soil 
temperature was measured by thermometers placed 
in glass tubes which penetrated the sod to the 


appropriate depth. 


Station operation was initially under the 
direction of the Agricultural Experiment Station 
at the University of Illinois-Urbana. With time, 
Agronomy Department personnel became interested 
in the operation of the site, as well as other 
faculty in the University, notably from the 
Depts. of Soils, and Geography. In August 1902 
the Urbana station became a U.S. Weather Bureau 
Cooperative Station, a status which continues to 
today. During the 1930s, interest in the Morrow 
Plots weather station relaxed somewhat, in part 
because a the University South Farms were growing 
in use. Indeed, a second observing site was in- 
stalled at the South Farms under the direction of 
the Agronomy Dept., the operation of which also 
continues until today. 


The Illinois State Water Survey has 
supervised the operation of the (now) National 
Weather Service Cooperative station since May 
1948. The site remained at the Morrow Plots 
until the summer of 1985, when it was moved to 
the Water Survey Research Center, located about 5 
km south southwest of the Morrow Plots site. 
Recording instruments continue to be maintained 
at the Morrow Plots, and temperature comparisons 
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2/purrill m.s., University of Illinois Library. 
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between the 2 sites will be made after several 
years of mutual data have been collected. 


MEMORABLE 19TH CENTURY WEATHER EVENTS IN ILLINOIS 


Records of early tornadoes in Illinois are 
sparse. After the Morman community left Nauvoo, 
a French communal group (Icarians) settled to 
restore the Morman temple. The temple was 
totally destroyed by a tornado 27 May 1850 
(Ludlum 1970). 


On 13 June 1857, Pana experienced a 
tornado which destroyed 25 houses, and essential- 
ly ruined the remaining buildings in the town 
from the effects of wind and hail. In another 
occurrence, nineteen persons were killed, and 50 
injured in the Ellison (Warren Co., population: 
500) tornado, 31 May 1858. 


Iudlum (1968) reports several severe 
winters and winter storms in Illinois. At times, 
it is difficult to determine whether the storms 
were more severe or more frequent than today, or 
whether the early settlers were simply inex- 
perienced with the severity of Midwestern winte- 
rs. Regardless, the winter of 1830-31 was 
intense and is now generally known as the Winter 
of the Deep Snow in Illinois. Records from Ft. 
Armstrong (Moline) and Ft. Dearborn (Chicago), 
and those of Dr. Samuel Mead of Augusta IL, a 
quality observer serving under the direction of 
the Smithsonian, all indicated cold temperatures 
for the winter as a whole, and heavy snowfalls of 
2 to 3 feet occurring at the end of December. 
This severe winter with a resulting poor corn 
harvest in northern Illinois, compelled northern 
Tllinoisans to seek grain in southern Illinois, 
who had largely escaped the extreme conditions of 
the north. The furnishing of grain by the south 
was reminiscent of the Biblical account, leading 
to the use of the term "Little Egypt" still given 
to southern Illinois. 


The winter of 1836-37 is known as the 
Winter of Starving Time in Illinois (Ludlun, 
1968). In December a strong cold front passed 
through northern and central Illinois dropping 
temperatures by more than 22°C in less than 24 
hrs. A temperature change of that intensity in 
Novembers, today, has a recurrence frequency of 
about once each 4 years (Wendland 1987). A man 
riding a horse from Chatham to Springfield to 
procure a marriage license began the trip in 
rain, and arrived in Springfield in snow with 
temperatures falling rapidly through 0°c. Upon 
arrival, his great coat was frozen solid, and he 
was so firmly frozen to the saddle, that he and 
the saddle were carried into a house and thawed 
next to the fire. Ft. Dearborn reported a mean 
December temperature (method of calculation 
unknown) of -6.7°C (4.2 below the 1951-80 mean), 
that of January was -5.1% (+0.8), that of 
February was -5.6°C (-2.2), and that of March was 
-2.0% (4.1). 





The winters of 1855-56 and 1856-57 were 
also memorable in Illinois. During the former, 
temperatures at Dubuque IA were less than 0°C for 
all observations from 21 December through 15 
February, 57 days. From 22 December through 11 
January, the maximum temperature at Sandwich IL 
(60 mi west of Chicago) never exceeded -17.2°C. 
The winter of 1856-57 was very cold in the Upper 
Midwest, with Ft. Ripley MN reporting -45.5°C on 
10 February 1857. 


THE STATE OF CURRENT CLIMATOLOGICAL DATA IN 
TLLINOIS 


The Morrow Plots site continues to be 
monitored for temperature and precipitation, 
however a full component of parameters is now 
continuously measured at the Water Survey Resea- 
rch Center. 


Somewhat more than 110 NWS Cooperative 
stations daily record temperature, precipitation, 
snowfall and snow on the ground today. These 
data are quality-controlled and published by the 
National Climatic Data Center (NCDC), Asheville 
NC. Copies of the original forms are on file at 
the ISWS, as is the published information from 


the NCDC. 


A comment is appropriate concerning 
changes in quality-control exercised on climate 
data through the years. Funding declines require 
declining services, and often, such cuts are made 
in "hidden" components, as in quality-control. 
Witnah (1961) gives several examples where 
Weather Service observations and services were 
cut during the mid-1800s for this reason. 


In addition to the Coop stations referred 
to above, NWS supports 5 first order stations in 
the state: O’Hare, Rockford, Moline, Peoria and 
Springfield. These sites operate 24 hours a day, 
monitoring virtually all components of the near 
surface atmosphere; strongly focused on flight 
operations. Because (1) these site are all 
located in the northern half of the state, (2) 
because of the increasing interest in alternative 
energy sources, and (3) because of the Water 
Survey’s commitment to monitor the whole water 
enviroment of the state, about 5 years ago the 
ISWS, supported by the Illinois Dept. of Energy 
and Natural Resources, installed 17 automatic, 
recording weather stations in the state, each 
monitoring temperature, precipitation, wind, 
solar, soil temperature and soil moisture. One 
need only recall the Water Survey’s timely 
assessments of the 1988 drought to realize only 
one of the benefits of such a network. 


The Water Survey was a pioneer in es- 
tablishing a means for the dissemination of near 
real-time climate data to the user public by es- 
tablishing the Climate Assistance Service (CLASS) 
several years ago (Changnon and others 1987). 
This computer-based system of data and informa- 
tion has been widely used by private citizens and 
agencies to assess the current state of climate 


and water resources. CLASS provided a foundation 
for a regional data and information system 
currently operated for the 9-state region by the 
Midwest Climate Center. 


The Water Survey represents the largest 
community of meteorologists and climatologists in 
one location in Illinois, indeed, in most states. 
It’s data, services and research have added much 
to the understanding of the atmosphere. We are 
committed to maintain that posture. 
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LONG~TERM CLIMATE TRENDS AND AGRICULTURAL 


PRODUCTIVITY IN SOUTHWEST MICHIGAN! 


James R. Crum2 
G. Philip Robertson 
Fred Nurnberger 


Long-term meteorological data (1888-1988) were obtained 
for the Kalamazoo weather station to describe the 
precipitation and temperature patterns for the last 100 
years in the area surrounding the W.K. Kellogg Biological 
Station's Long Term Ecological Research site. Also, corn 
(Zea mays) yield data were obtained from the Michigan 
Department of Agriculture to examine relationships between 
the amount and timing of precipitation and harvested corn 
yield for the period 1945-1985 for Kalamazoo County, 
Michigan. 


Long-term temperature trends for Kalamazoo County, 
Michigan show a 0.5-1.0 °C increase in mean annual 
temperature since 1930. Nonurban records of less duration 
from a nearby rural site do not show this increase in 
temperature; possibly indicating this increase may be due 
to urban heat load from the City of Kalamazoo. Annual 
precipitation patterns show little change over the 1888- 
1988 period other than a greater frequency of above-normal 
amounts for years post 1930, and the strongest correlation 
between precipitation and yield occurred for the June-July 
precipitation period. 


Keywords: climate, corn, soils, variability 
agricultural productivity, Michigan 


INTRODUCTION dependent to one degree or another on the amount 
and timing of rainfall to satisfy crop demands 
Agricultural ecosystems differ from most for water, even where supplemental irrigation is 
natural systems in their disturbance regimes, in used. In Southwest Michigan, annual field crops 
the restricted diversity of both plant species (mainly corn and soybeans) are planted in late 
and plant life history patterns, in the levels April and May and harvested in October or 
of inputs (e.g. water, fertilizers, pesticides) November. Rainfall before and within this 
added to the systems, and in their degree of period is critical to crop success. Too much 
human management. Agricultural ecosystems are precipitation in the spring delays planting, 
managed to produce a crop that will return an decreasing the length of the growing season and 
economic profit. In all systems production is grain yield. Too little precipitation delays 


crop emergence and makes the maturing crop 
susceptible to late season drought or frost. 
Midseason rainfall, during pollination and 
grain-fill periods, is equally critical: too 


Ipresented at the 1988 LTER Climate Variability little precipitation from late July to early 
Workshop, University of Colorado, Mountain September drastically reduces yield under rain~ 
Research Station, Boulder, Colorado, August fed conditions and too much precipitation delays 
1988. harvest and losses occur from stock lodging and 
grain spoilage. 
2associate Professor of Soil Science, Michigan Between rainfall events ail ecosystems 
State University, East Lansing, MI; Assistant depend on water storage in soil. Soils vary 
Professor of Soil Ecology, Kellogg Biological greatly in their ability to store water and 
Station, Michigan State University, Hickory buffer the effects of prolonged dry 
Corners, MI; and State Climatologist, Michigan conditions. Thus the timing of rainfall events 
State University, East Lansing, MI. is especially critical in well-drained soils 
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such as those typical of SW Michigan. 

In this paper we 1) describe precipitation 
and temperature patterns for the last 100 years 
in the area surrounding the W.K. Kellogg 
Biological Station's Long Term Ecological 
Research (LTER) site in agricultural ecology, 
and 2) examine relationships between the amount 
and timing of precipitation and harvested corn 
yield since 1945 for the county in which the 
LTER site is located. 


SITE DESCRIPTION AND METHODS 


Soils and Agriculture of Southwest Michigan 


Fourteen thousand years ago Michigan's 
Lower Peninsula was covered by glacial ice 
advancing from three basins: (1) the Lake Erie 
basin to the southeast; (2) the Saginaw-Huron 
basin to the due east; and (3) the Lake Michigan 
basin to the west. Upon retreat of the glacial 
lobes, some 12,000 years BP, glacial outwash 
(mainly sand and gravel) was deposited over much 
of Southwest Michigan, including the area around 
the Kellogg Biological Station (KBS). A variety 
of soils have formed on this sandy-textured 
parent material: Mollisols in areas of flat 
topography that supported prairie vegetation and 
Alfisols in areas with more rolling terrain. At 
KBS, principal soils are Alfisols of the 
Kalamazoo series (Fine-loamy, mixed, mesic Typic 
Hapludalfs) and Oshtemo series (Coarse-loamy, 
mixed, mesic Typic Hapludalfs). These two soils 
are quite similar, differing mainly in the 
thickness of the clay loam textured upper 
argillic horizone The Kalamazoo has a thicker 
upper Bt horizon than the Oshtemo and therefore 
has a greater water holding capacity. 

Agriculture in southwest Michigan is as 
varied as the soils of the region. Major 
agricultural commodities include grain and 
forage crops (corn, soybeans, wheat, oats, hay) 
fruit crops (apples, cherries, peaches, grapes, 
pears, plums, blueberries) vegetables 
(asparagus, carrots, cauliflower, celery, 
cucumbers, tomatoes, lettuce, onions) flowers 
(bedding plants, potted plants, cut flowers) and 
animals (beef and dairy cattle, hogs, 
poultry). Michigan ranks among the top five 
states for many of the above crops (particularly 
fruits and vegetables) and among the top 10 for 


diversity are 1) the moderating climatic effect 
of Lake Michigan, 2) inherent soil fertility, 
and 3) proximity to several large population 
centers. 


General Climate of SW Michigan 


KBS is about 65 km north of the Indiana 
border and 95 km east of Lake Michigan. It is 
east of Michigan's fruit belt and the snow belt 
that extends north-south along Lake Michigan's 
eastern shore. The climate is primarily 
continental in character but the marine 
influence of Lake Michigan is sometimes 
expressed (Strommen 1971). The most noticeable 
lake effect, the result of prevailing westerly 
winds, is increased cloudiness and milder 
temperatures during the fall and winter months. 


Long-term Meteorological Data 


The long-term meteorological data (1888- 
1988) on which much of our results are based 
were obtained from the Kalamazoo, Michigan 
weather station located approximately 15 km 
southwest of the KBS LTER site. The station was 
opened in June of 1867 on the grounds of the 
State Regional Psychiatric Hospital (85 °, 36 
w; 42 °, 17 'N 3; elevation 288 m). It has 
been in continuous operation since opening. The 
Hospital is located on one of the highest 
topographic positions in Kalamazoo (population 
85,000); effects of urban heat load on the 
temperature patterns are unknown. 


y 


RESULTS AND DISCUSSION 


Long-term Trends 


Because day-to-day weather in SW Michigan 
is controlled largely by the movement of 
pressure systems across the continent, KBS 
seldom experiences prolonged periods of either 
hot, humid weather in the summer or extreme cold 
during the winter (Table 1). The 30-year (1951- 
1980) mean annual temperature is 9.2 °C and the 
mean annual precipitation is 86.2 cm (Table 
1). Figures 1 and 2 show the mean annual (MA) 
temperatures and total precipitation for the 
years 1888-1988. The deviation from the 
calculated 30 year normal (1951-1980) is quite 


| 
| 
| 


most. Major reasons for this agricultural large for both MA temperature and precipitation, 


Table 1. Temperature and Precipitation Normals for the Kellogg Biological 
Station, 1951-1980. 


Seb es tet it oe ton cm ete see Sh Se A SN ns an st SS ei DH mt mM Me ts eh le ot mE HD MO MD A te A 


ee ae ce a ns th Me ME ‘th et el St Nn a le lt mise tenth i el es ma a Wm Sn a Heh Nuh Nh Ch el OS A sh tute hyo tem in Nt Yn: mech 


Temp(C) 
max “0.9 1.0 6.6 14.9 21.6 26.6 28.8 27.9 24.0 17.4 8.7 1.8 14.9 
min ~9.6 -8.9 -4.0 2.3 8.1 13.4 15.7 14.9 11.1 5.4 -0.2 -6.1 3.5 
mean ~5.2 -3.9 1.3 8.6 14.8 20.0 22.3 21.4 17.6 11.4 4.3 -2.2 9.2 
Precip(mm) 
mean 45 36 51 89 80 106 86 89 75 73 68 57 855 
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Figure l. Mean annual temperatures (a) and 
5-year running means (b) mean annual term for 
the longest-running Kalamazoo County weather 
station over the period 1888-1988. 
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Figure 2. Mean annual precipitation (a) and 


5-year running means (b) of mean annual 
precipitation for Kalamazoo County, Michigan. 


suggesting high year to year variability. Five- 
year running means suggest that both temperature 
(Fig. lb) and precipitation (Fig. 2b) may be 
gradually increasing: recent temperatures 
appear to be .5-1.0 °C higher than in the 1880- 
1930 period, and post~1930 precipitation means 
are more frequently above the 100 y norm. The 
effect of urban heat load on this trend is 
unknown but potentially significant; available 
data (1930-present) from the Gull Lake weather 
station (located on the Kellogg Biological 
Station grounds) shows no consistant increase in 
temperature or precipitation since 1930. 
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Soil Water Balance 


The primary soil on the KBS LTER site is in 
the Kalamazoo series. Available water-holding-— 
capacity (WHC) of this soil is approximately 150 
mm to 1.5m depth. The soil-water balance 
(Fig-3) indicates that precipitation is greater 
than potential evapotranspiration (PET) from 
about the middle of September through late 
April, and that stored soil water together with 
precipitation is used through the summer 
months. Typically a soil-water deficit occurs 
in August, with a duration dependent on the 
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timing of summer rainfall. For example in 1988 
KBS received below normal precipitation from May 
through mid-July and then above normal for the 
remainder of the year. A soil water deficit 
occurred early and strongly affected crop 
yields, though soil-water recharge occurred the 
remainder of the year. Despite the growing- 
season drought, soils were at or near field- 
capacity at the beginning of the 1989 crop year. 
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Figure 3. Monthly precipitation means and 


annual soil water balance for a Kalamazoo Soil, 
150mm water-holding-capacity. Precipitation and 
PET based on 30-year (1951-1980) means. 


Yield vs. precipitation. 


We used Michigan Agriculture Statistics for 
Field Crops (Borum et al., 1950, 1960; Fedewa et 
al., 1983, 1986; Hines et al., 1970, 1976, 1980) 
for the years 1945 through 1985 for Kalamazoo 
County corn yields together with Kalamazoo 
weather data to test the hypothesis that annual 
amounts of precipitation are not as important as 
precipitation that occurs during the crop 
growing season, when PET is greater than 
precipitation (Fig. 3). 

Corn yields have increased remarkably in 
Kalamazoo County since 1945 (Fig. 4). This 
increase is due mainly to increased fertilizer 
use, improved crop varieties, and the use of 
herbicides and pesticides. In order to relate 
yields to precipitation patterns, we attempted 
to remove such changes in management as a factor 
by expressing Kalamazoo yields as a ratio of 
Kalamazoo yields to yields from Lenawee 
County. Lenawee County, located in SE Michigan 
near Lake Erie, is the most productive (greatest 
yield per area) county for corn in Michigan 
because the soils generally have a higher water- 
holding capacity which better buffers against 
short-term drought. Our hypothesis predicts 
that changes in the Kalamazoo:Lenawee production 
ratio are due strictly to climatic factors. 
Management factors are presumed to be equivalent 
for the two counties. Monthly weather data for 
Lenawee were accumulated from the Adrian 
Station, located in the geographic center of the 


56 


etd (bu//on) 


Yielet Ratic 


140 ecg aa 
Larkiwes 
1Z¢ 
Kalomecs 
100 
80 
66 
49 
20 
g = ~ C ew 
1940 1956 1986 ly7d 1980 bos 
aor 
Figure 4. Corn yields 1945-1985 in Kalamazoo 


and Lenawee counties. 


county. To remove differences due to total 
annual precipitation we have also expressed 
precipitation as a ratio of Kalamazoo vs. 
Lenawee counties. 

We found very weak correlations between 
rainfall and corn yields for both annual (Fig. 
5) and key-month (Fig. 6a-d) periods. 
Nevertheless, correlation was stronger (r=0.36, 
n=41, p<.025) for the July-August period (July- 
August precipitation ratios vs. annual yield 
ratios) than for either annual precipitation 
(r=.22, n=41, p>.5) or precipitation in any 
other 2~month period (r<.06). This suggests 
that annual yields are indeed most sensitive to 
precipitation deficits in the July-August 
period, as predicted above based on our 
knowledge of the soil water balance (Fig. 3). 
Earlier in the growing season precipitation is 
not as important a predictor because stored 
soil-water is generally available and the 
plant's water demand is not as great. 





Precipitation Ratic 


Figure 5. Ratios of corn yield (Kalamazoo: 
Lenawee counties) vs. annual precipitation 
(Kalamazoo:Lenawee) for the period 1945-1985. 
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Figure 6. Ratios of corn yield (Kalamazoo: 
Lenawee counties) vs. precipitation for the 
periods of a) April-May, b) May-June, c) June- 
July, and d) July-August, 1945-1985. 
CONCLUSIONS 
That our correlations are not stronger 1. Long-term temperature trends (1888-1988) for 
could be due to a number of factors. First, Kalamazoo County, Michigan show a 0.5-1.0 6 
precipitation data are from a single county increase in mean annual temperatures since 
station while yield estimates are from the 1930. Some portion of this increase is 
entire county. Summer precipitation is mainly likely due simply to urban heat load at the 
in the form of afternoon showers and weather station site. Nonurban records of 
thunderstorms (Strommen 1971) with relatively shorter duration (Gull Lake, 1930-1988) do 
small impact areas. A single weather station not show consistantly increasing temperature 
may thus only roughly represent county-wide for the post-1930 period. 
means. A second source of error is county-wide 
differences in water holding capacities among 2 Annual precipitation patterns show little 
soils under corn cultivation. Small differences change over the 1888-1988 period other than 
in WHC will significantly affect a crop's a greater frequency of above-normal 
ability to weather short-term, midseason precipitation years post-1930. 
drought, and thus have a significant effect on 
the relationship of midseason precipitation to 3. County-wide corn yields and single-station 


county-wide yields. Finally, bimonthly (e.g. 
June-July) means may not sufficiently resolve 
precipitation patterns that affect corn 
growth. Corn is most susceptible to drought 
stress during tasseling, silking and early ear 
growth. 


precipitation patterns were weakly 
correlated. The strongest correlation 
between precipitation and yield occurred for 
the June-July precipitation period. The 
soil water deficit is also highest for this 
period. 
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GLIMATE VARIABILITY AT NIWOT RIDGE IN THE TWENTIETH CENTURY 


David Greenland 4 





Abstract.--The climate of Niwot Ridge displays abnormally low 
temperatures in the early 1980's and an overall downward trend 
in temperature. Some cyclicity appears in the data series of 
annual values of both temperature and precipitation. Silver 
Lake annual precipitation data show only a slight downward trend 
when the extremely high value for the year of 1921 is excluded. 
Implications of these changes for the alpine tundra ecosystem 


are briefly discussed. 


Keywords: Climate Change, Climate Analysis, LTER, Ecosystem 


Response. 


INTRODUCTION 


The climate of Niwot Ridge, Colorado is espec- 
ially important because of the large amount of re- 
search in physical geography and ecology which is 
performed there. It has been the site of research 
in the International Biological Programme, the 
Man and the Biosphere Program, and, more recently, 
the Long-Term Ecological Research Program (LTER). 
It is over 16 years since Barry (1972, 1973) pro- 
vided a much quoted analysis of the climatic ob- 
servations at the site. Since that time meteoro- 
logical observations have continued to be made and 
it is fitting that a new analysis of the record 
be undertaken, not only because of the increasing 
frequency of LTER-related and other studies at 
the site but also because of the current general 
interest in climatic variability and global 
change. 


This paper focuses on monthly mean and annual 
values of temperature and precipitation for the 
Cl and Dl stations on Niwot Ridge and on other 
nearby stations. The establishment of the Niwot 
Ridge stations is described by Greenland (1987). 
Data associated with the Niwot Ridge stations are 
presented for 1952-1987 and additional precipita- 
tion data from Silver Lake, 1 km. southwest of 
Cl, are analyzed from 1915. The Ridge, along 
with the nearby Green Lakes valley to the south, 
is the locus of research operations in the LTER 





1/ Presented at the LTER Workshop on Climate 
Variability, Niwot Ridge Green Lakes Valley LTER 
Site, CO, August 21-23, 1988. 


2/ Professor of Geography, Department of Geog- 
raphy, Campus Box 260, University of Colorado, 
Boulder, CO 80309. 





program. The site and its synoptic climate is 
described briefly. Most attention in the climate 
analysis is given to the years 1951-80, the time 
selected by the LTER Intersite Climate Committee 
for cross site comparison. Considerable effort 
was given to estimating early and missing values 
in the Niwot Ridge data using values from the 
nearby Allenspark observation site. The resulting 
data series are analyzed to give a climatic de- 
scription of the site including information on 
the long-term average water balance and temporal 
variation. Station data are: 





Station Latitude Longitude Elevation 
Degrees N. Degrees W. Meters 

Como, 

Cl 40 02' 105.32! 3048 

Niwot Ridge, 

DL 40 03' £05.37" 3749 

Boulder 40 02' 105-6" 1638 

Allenspark A122 * TOS 32 2591 

Silver Lake 40 02' Lo 3a 3109 


Readers who are interested in climatic variables 
and values other than monthly mean and annual 
temperature and precipitation should refer to 
other sources listed in the Front Range alpine 
and subalpine bibliography compiled by Halfpenny 
et al. (1986). A more comprehensive version of 
this paper has been presented elsewhere 
(Greenland, 1989). 


THE SITE AND ITS SYNOPTIC CLIMATOLOGY 
The Niwot Ridge site is in an alpine tundra 
ecosystem. The Ridge stretches eastwards from 


the Continental Divide and the area of most 
studies is found between the Cl station at 3048 m 
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and the Dl station at 3750 m with the treeline 
being found at approximately between 3350 and 
3500 m. On the ridge and in the Green Lakes 
valley to the south, there are many distinct 
topoclimates associated with such factors as 
saddles and knolls, moraines and other glacial 
and periglacial features, semipermanent snow 
banks, and permanent ice and lakes. 


The synoptic scale climate of the area is 
controlled by the mid-latitude, continental lo- 
cation and by the elevational and topographical 
situation. The high elevation causes low air 
temperatures at all times of the year. Air 
temperatures are effectively further depressed 
by the high velocities of the wind which gives 
vise to low wind chill values. The mid-conti- 
nental location leads to a large temperature 
range between summer and winter but this large 
range is more marked at the lower elevations. 
Precipitation-carrying storms are brought over 
the site in the winter and spring by the upper 
westerly flow. In these seasons, snow is 
brought from the west at higher elevations. 
is also brought from the east, at the lower 
elevations, by cyclonic easterly, upslope flow 
developing to the east of the Continental Divide. 


Snow 


Table 1 
SUMMARY STATISTICS O1 


TEMPERATURE 
eg. C. 
Jari Feb Mar April May 
Mon Mean 15.15 ~le.?F% 11.23 ~#.01 -0. 66 
An Mean ~3.71 St Dev 6.06 
Mean Mx T  -10.09 ~9.61 ~?.60 ~3.32 2.68 
Mean Mi T -16.21 ~15.94 -14.66 ~-10.70 ~4.40 
Mean Temp Warmest Month 8.23 St Dev 0.89 
Mean Temp Coldest Month “13.15 St Oev 1.35 
Anrwal Range of Monthly Mean Temps 21.36 
Num months with mean temp >0 4 
Num manths with mean temp >15 0 
Highest monthly mean in record per 10.48 
Lowest monthly mean in recerd per ~1?.58 
PRECIPITATION 
mat 
Jan Feb Mar April May 
Mor mean 102.5 80.1 127.68 100.1 87.5 
Mean annual total 930.0 
Hettest year in period l4er 
Oriest year in pericad 541 
Monthly totals during wettest year in period Year 
Jar Feb Mar April May 
209 121 255 220 45 
Monthly totals during driest year in period Year 
Jan Fab Mar April May 
ee 20 54 52 el 
Total precip in months with temp of eel 
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100 


June 
16 


These storms are responsible for the spring 
maximum of precipitation. In the summer, rain- 
fall is produced from localized convectional 
storms. Fall is the driest season. 


METHODS OF DATA PROCESSING AND ANALYSIS 


Missing and Biased Data 





There are two principal problems in creating 
monthly and annual summaries for the Cl and DL 
stations. First, there are occasional data gaps 
when data were simply not collected. Second, 
until November 1964 there was no wind shield 
around the precipitation gauge at the Dl site. 
This gave rise to significant under-catch in pre- 
cipitation, particularly in the winter. 


Estimates of the missing values caused by the 
above problems were given by standard climato- 
logical techniques of using regression analysis 
or by ratios with observed values from a nearby 
station (Brooks and Carruthers, 1953). The most 
suitable station for use in this analysis is that 
of Allenspark located 28 km almost due north of 
the Cl station. It is the nearest station that 
has a long enough climate record and is in the 


July Aug Sept Det Now Dec 
8.235 7.06 3.3559 ~2.13 “6.94 -11.77 
12.43 11.00 7.35 1.34 -5. 99 -B.46 
4,03 a.1i1 ~0. 57 ~5.59 -12.50 ~-15.07 
July Aug Sept Get Now Dec 
50.5 Be. 1 49.1 40.3 bo. 3 30.1 
ard 
July Rug Sept Bet Nov Dec 
36 12 495 it4 73 15? 
13954 
July Rug Sept Bet Nov Bec 
61 45 64 35 55 43 








same climatic regime as the Niwot stations with 
its high latitude, east of the Continental Divide 
location. The Allenspark data and their treat- 
ment are described by Greenland (1987). During 
the analysis of the Niwot data it was determined 
that any slight discrepancies in values quoted by 
Marr (1961) and Barry (1972, 1973) are due to 
variations in the ways in which missing values 

for individual days are treated. Estimated 

values of monthly precipitation for the Dl station 
for the years 1951-1964 are used in the data pre- 
sented in Table 1. A new study comparing precipi- 
tation values from unshielded and shielded rain 
gauges for a limited time during the 1951-1964 
period suggests that the synthetic values are 
quite reasonable (Clark and Halfpenny, pers. comm. 
1989). 


Water Budget Computations 





Water budget computations were made based on 
the Thornthwaite technique (Thornthwaite and 
Mather, 1957) and using an adaptation of a program 
by Willmott and Rowe (1985), which was based on 
and is more fully described by Willmott (1977). 
The technique is not fully applicable for this 
area owing to its neglect of such factors as 
windspeed, humidity conditions, and a direct net 
radiation input. It is useful, however, to give 
a relative comparison between the Dl and Cl sites 
and also between other LTER sites. The latest 
version of the program, which incorporates snow- 
pack water equivalent and snowmelt, gives a 
vivid impression of the importance of snowmelt in 
the Niwot Ridge area. 


Time Series Analysis 


The data were preprocessed as follows. Month- 
ly data values in Lotus worksheet files were con- 
verted to ASCII files. The ASCII files were up- 
loaded onto the University of Colorado VAX 8550 
cluster. SAS data sets were created for time 
series modeling and ASCII data sets of annual 
data values were created for spectral analysis. 


The prepared data sets were plotted as time 
series and examined. Trends were investigated by 
fitting quadratic regressions. Exploratory 
modeling of the annual data values took the form 
of identifying the best fitting Autoregressive 
(AR) and Autoregressive Integrated Moving 
Average (ARIMA) models using an iterative process. 
Autocorrelation functions (ACF) and Partial 
Autocorrelation functions (PACF) were used to 
guide initial estimates of parameters. 


In order to perform the spectral analysis first 
the time series were smoothed and detrended. 
Next the ACF of the detrended series was checked 
for Markov persistence using the first three 
values of the ACF. Confidence intervals for which 
Markov persistence did not exist were dismissed. 
The spectral peaks in frequency cycles per year 
were translated into periods. The Nyquist fre- 
quency prevents the detection of periods of less 
than 2 years. Most of the present analyses show 
a sharp rise at the two year period perhaps due 
to the l-year autoregressive component evident in 
most of the time series. 


Annual precipitation values for Silver Lake 
were prepared for time series analysis in the 
following manner. Unpublished microfiche data 
for 1910 to 1985 and published NWS data for 1910 
to 1950 were combined to form what was believed 
to be the most accurate data set. Data from 
Allenspark for 1945 to 1985 and from Cl from 1951 
to 1987 were assembled for regression analyses to 
help simulate missing data values. There were 
too many missing values between 1910 and 1913 so 
the final modified data set which was analyzed 
consisted of the 74 data points between 1914 and 
1987. The analysis was then performed in the 


same manner as for the Cl and Dl data as previous- 
ly described. 


DISCUSSION OF RESULTS 


Climate Description 


Summary statistics for Dl and Cl for the years 
1951-80 are presented in Tables 1 and 2. These 
clearly show the severity of the climate at Dl 
on the tundra. Important characteristics include 
(1) the existence of only four months with temper- 
ature above freezing, (2) the large annual range 
of monthly mean temperatures, (3) the relatively 
large variability in intra-annual precipitation, 
(4) the winter and spring precipitation maximum, 
and (5) the small amount of precipitation arriving 
in a "growing season" defined as one having 
monthly mean temperatures above O°C. Cl station, 
below the treeline, shows (1) an extended grow- 
ing season, (2) an almost equally large range of 
monthly mean temperatures, and (3) a more clearly 
defined spring precipitation maximum. The rela- 
tively small amount of precipitation arriving in 
the "growing season" is misleading in this case 
because of the importance of snowpack. Whereas 
much of the snowfall arriving at Dl is blown away 
and never becomes available to the plants, at Cl 
winter snowfall largely resides as a snowpack on 
the ground between the trees and becomes avail- 
able to the vegetation during the spring melt. 


Water Budget 


Water budget analysis suggests that about half 
of the precipitation at Cl is lost by evapo- 
transpiration while the rest is lost by surface 
runoff of percolation into the soil. At Dl, in 
contrast, only about a third of the precipitation 
is lost by evapotranspiration, while the rest of 
the greater amount of precipitation is lost as 
runoff or, more realistically, percolation into 
the soil. The importance of the spring thaw is 
seen at both sites. 


Temporal Variation 


Some of the most interesting findings of the 
present study are associated with the time series 
analysis of the data. Although analyses were 
performed for both the Cl and the Dl stations only 
diagrams for the latter are described here in de- 
tail. Pertinent points relating to the Cl ana- 
lyses are summarized but in general the same over- 
all characteristics appearing in the Dl data are 
also seen in the Cl data. 
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Table 2 
SUMMARY STATISTICS Cl 


TEMPERATURE 
Deg. C. 
: 5 c Nov Dec 
Jan Feb Mar April May June July Rug Sept _Get v 
Mon Meari ~7.43 -6.70 -5.51 -1.24 3.75 8.63 11.91 10.95 7.40 3.16 -3.26 -6.02 
Ari Mean 1.30 St Dev 6.60 
5 2 : 5 ? 4.13 9.26 1.78 ~1.41 
Mean Mx T -3.09 -1.49 ~-0.21 4.06 3.59 15.47 19.158 17.91 14.1 ; ' A 
oon Mi T 3-11.90 -11.91 -10.80 -6.55 -1.89 1.78 4.64 3.99 0.67 -2.96 ~B.30 10.62 
Mean Temp Warmest Month 11.91 St Dev 0.96 
Mean Temp Coldest Manth ~7.49 St Dev 1.58 
Arirwal Range of Monthly Mean Temps 19.40 
Num months with mean temp >U 6 
Num months with mean temp >15 0 
Highest monthly mean in record per 13.75 
Lowest monthly mean in record per -11.34 
PRECIPITATION 
mim 
z U ‘ Bec 
Jan Feb Mar April May June July Rug Sept _Uct Nov 2 
Mon mean 50.3 45.2 67.2 85.4 85.3 54.8 65.3 59.3 50.2 37.0 46.9 45.0 
Mean annual total 6592.0 
Wettest year in period 949 
Uriest year in period 3935 
Manthly totals during wettest year in period Year 195? - : 
es Jari ‘ Feb Mar April May June July Aug Sept Det Nov Dec 
65 46 S2 199 185 b4 gg 76 43 be 29 1 
Monthly totals during driest year in period Year 1954 
4 Jan e Feb Mar April May June July Aung Sept Oct Now Dec 
13 10 33 20 47 3 78 42 @l 28 Je 18 
Tatal precip in months with temp >U 352 


Annual temperatures at D1 (Figure 1) are note- 
worthy for their year to year variability from 
1951 to 1980 with a standard deviation of 0.64 
°C around a mean of -3.71 °C and for the relative- 
ly low temperatures in the early 1980s. A quad- 
ratic regression curve fit to the data series 
indicates a general decreasing trend with the 
steepness of the descent being greater in the 
second part of the series. There is a linear 
trend over the data period of -0.0013°C per year 
for the first part of the period from 1951 to 
1980. Although Barry (1973) did not find a signi- 
ficant trend for the Dl data, he did find a down- 
ward trend of -0.07°C per year for the Cl data. 
The linear trend for the whole data set including 
the low temperature values of the early 1980s 
yields a downward rate of - 0.05 °C per year. 
Five year running means of this data series 
emphasize low temperatures of the early 1980s. 
There also is a suggestion of a decreasing trend 
from 1953 to the mid-1960s and a slight increase 
from here to the beginning of the 1980s. 


It is important to consider whether the period 
of low temperatures in the early 1980s, which 
might be called the Niwot Minimum, is real or not. 
Since the Minimum appears at both Dl and Cl, it 
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is unlikely that the feature is due to some charac- 
teristic of one of the sites. The observing 
personnel did change in 1980, but calibration pro- 
cedures and instrument did not and the new ob- 
server is not aware of any instrumental or pro- 
cedural reason for the appearance of the Minimum. 
The likelihood of observational or instrumental 
error would be decreased if similar changes were 
seen in nearby, independent systems. Unfortunate- 
ly there are many missing data for the Allenspark 
station between 1979 and 1983 and in 1984 its 
location was changed from a ridge site to a 

valley site 46 m lower and 2.4 km to the north. 
Thus, temperatures at Allenspark for 1984 and 1985 
that are significantly lower than the 1951-1980 
mean, although being consistent with the Niwot 
Minimum, may reflect the site change. Annual 

mean temperatures at Nederland (8.5 km south 
southwest of Cl) for 1982-85 are about 1.5°C 

lower than in 1980 and 1981, but are higher than 
in the previous decade. Annual mean temperatures 
at Boulder during the Niwot Minimum are lower 

than in the previous three decades. Other pieces 
of information that support the reality of the 
Minimum exist. First, there were positive mass 
balances of ice in many of the glaciers near 

Niwot Ridge in the years of the early 1980s 
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Figure 1. Annual Mean Temperatures at Dl Station 1950-1987 


(Caine, pers. comm. 1986). Second, in the Saddle 
(on Niwot Ridge about halfway between Cl and D1), 
snowpack remained unusually late in 1982 and 1984. 
In 1984, the summer snowpack did not melt out 
until the winter snowfalls came. Thus, although 
it will have to remain a judgement call, the event 
is more likely to be real than not. 


One of the important implications of the Niwot 
Minimum, if it is real, is that any search for 
warming due to global greenhouse gas enhancement 
will be unlikely to be fruitful at the Niwot site 
because of the influence of the temperature values 
of the early 1980s on the time series analysis. 


A time series model forecast provides two pieces 


of information. First is the forecast itself and 
second is time series diagnostics that go to pro- 
duce the best forecast model. Autoregressive 
Integrated Moving Average Forecasting Models 
(ARIMA) are characterized by the statistics p, d, 
and q where p is the number of periods in an 
analysis where a partial autocorrelation function 
is used, d is the number of difference times em- 
ployed to render the time series stationary from 
an otherwise potentially nonstationary one, and q 
is the number of periods (years) in the weighted 
moving average. Thus, in a general model such as 
ARIMA (p,d,q) the values of p and q will be 


somewhat indicative of the potential cyclicity in 
the data (in this case in years). 


The best forecast model for Dl annual temper- 
ature values simulates trends quite well although, 
typically for these kinds of models, under- 
estimating peak and low values. There is also a 
one year lag of forecast to observed values 
(with observed values occurring one year before 
they are forecast by the model). It is not 
possible to estimate what confidence to attribute 
to the forecast. With this caveat in mind, it is 
worth noting that the model projects higher 
temperatures for the turn of the decade followed 
by somewhat lower values in the first half of the 
next decade and higher values in the second half. 
Model parameters indicate some importance of 2 and 
5 years value recurrence at Dl and 2, 13, and 14 
year recurrence at Cl. 


Spectral analysis suggests significant sta-~ 
tistical periodicities in Dl temperature values at 
9.00 and 7.20 years and in Cl temperature values 
at 12.00 and 9.00 years. High-frequency cycles 
near two years are believed to be a function of 
the methodology rather than reality. 
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Annual precipitation values at Dl demonstrate 
high inter-annual variability. The 1983 value, 
which occurred during the first five years of LTER 
studies at the site, was the highest in the 37- 
year record. The summer (JJA) precipitation for 
1984 was also the highest of all comparable values 
in this record. There is a small upward trend in 
the record when fitted with a quadratic regression 
curve. Although when scaled to emphasize trend 
the 5-year running mean does suggest an upward 
trend with dry years in the mid 1950s, around 
1970 and in the late 1970s. The best ARIMA fore- 
case model does well in simulating high and low 
values, but shows a lag of one year between fore- 
cast and observed values. It does not project any 
large changes into the next century. Spectral 
analysis of the Dl annual precipitation record 


manifests significant peaks at 7.19 and 9.00 years. 


The Silver Lake annual values of precipitation 
again show high year to year variability but are 
important in that they are able to place some of 
the Niwot Ridge values into a better temporal per- 
spective. For example, 1983 is shown here to have 
only the second highest precipitation value in the 
years back to 1951. More importantly an idea of 
the largest possible precipitation for the area is 
shown by high values in 1921 about which more will 
be said later. A linear regression of these data 
shows only a small downward trend but the 5-year 
running means suggest that this may be due to the 
1921 value and that the upward trend previously 
seen in the Dl data might be continuous from the 
dry period of the mid-1920s. The best ARIMA model 
forecast does less well in simulating amplitude 
and suggests only a continuation of the small up- 
ward trend through to the end of the century. No 
highly significant low frequency periodicities 
are shown by the spectral analysis although in- 
teresting high spectral values significant at the 
80% confidence level are found at 14.71 and 12.35 
years. The last of these parallel the 12.05 year 
periodicity found in Dl precipitation at the same 
confidence level. 


The year 1921 is identified as an extraordinary 
high precipitation year at Silver Lake and pre- 
sumably in the Niwot Ridge area. Monthly values 
(Table 3) indicate winter, spring, and the month 
of June had high values. Potential snow equiva- 
lents at 1/10 and 1/15 ratios for the winter 
months are also provided. If all precipitation 
fell as snow in April and none melted or blew 
away and the higher conversion ratio is used the 
15.8 m (51.84 ft) falling on the previous snow- 
pack would certainly have buried much of the 
vegetation. 


Implications for the Ecosystem 


It is difficult to say exactly how the climate 
changes described above would affect the eco- 
system, if at all. 


Some preliminary points are in order. First, 


it is not necessarily annual values of temperature 
and precipitation that have the most marked effect 
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on the alpine ecosystem. In many cases, the snow- 
pack development in the spring, and the growing 
season conditions are the most important. To the 
extent that annual values are indicative of these, 
then the annual values might have greater utility. 
With respect to this it has been shown that the 
summer temperature values play a large role in 
affecting the annual temperatures for Cl and D1 
(Greenland, 1987). Second, it is hard to deter- 
mine the effect of runs or patterns of years 

with similar values as opposed to individual 

years of high or low values. Greenland et al. 
(1985a) showed that the 1983 and 1984 years show- 
ed only smail differences in above -— ground net 
primary productivity despite quite large differ- 
ences in summer soil moisture. The fact that the 
majority of the alpine tundra biomass is below 
ground will create a buffer to inter-annual and 
shorter climate variations. Greenland’ (1984) 
further demonstrated that values of the surface 
summer energy budget in the alpine tundra could 
vary more in one year between the moist and dry 
surfaces at the Saddle than between a dry year 

and a wet year over the same surface. It is quite 
likely that it is the early growing season snow- 
melt and rainfall which is important for plant 
growth and survival. Snowpack distribution and 
ablation rates are critical in the phenology of 
the vegetation communities at the Saddle. 


The above points go some way to explaining the 
fact that there is at present no evidence for 
dramatic changes, such as species disappearance, 
in the vegetation of Niwot Ridge over the last 
40 years or so. Neither the long term trends in 
the annual temperature and precipitation data, 
nor the cyclicities, nor the periods of several 
or individual years with high or low values of 
these variables appear to be related to any ob- 
vious changes in the vegetation. This testifies 
to the resilience of the flora and how it is fair- 
ly well adjusted to the normal year to year 
variability of the present day climate. This is 
not to say, however, that there have been no 
changes in the vegetation associated with some of 
the climate changes described here. The vegeta- 
tional changes, if they exist, may be quite 
subtle and concerned with variations in the re- 
lative proportions of different species existing 
at any one time at a locality on the Ridge. Some 
scientists who have been working on the Ridge in 
the last decade believe they have noticed such 
changes (Halfpenny, pers. comm., 1988). Current- 
ly, a very detailed Geographic Information 
System is being constructed which should have the 
ability to detect such changes if they occur. 


There are certainly many questions to answer. 
Are the correct climatic variables being measured? 
The answer to this lies in what ecological 
questions require resolution. Annual values of 


Gren tand, D.E. 1984. Energy budgets over a 
moisture gradient in alpine tundra. Paper 
delivered at Association of American Geographers 
Meeting. Washington, D.C. 











Table 3. Precipitation 
Silver Lake Elevation 319097 m (10,200 feet?) 
Annual Average precip 1910 - 1930 847 mm 
1921 anmual precip I725 mm (48.20 ins) 
L921 Monthiy Precio Values 
Inches 
d F M a Pl ul d 
7,51 6.65 G.31 15.16 3.909 12.36 4,15 
nm 
L186 Loe aid Bes 7 S14 1oS 
Snow equivalent in meters at 1/10 ratia 
Leg bet wal a? 
Snow equivalent in meters at 1/1 ratio 
2.8 ae Sue eS 


temperature and precipitation address broad 
questions such as what is the range of values of 
the climatic environment to which the vegetation 
and fauna is adapted. Day to day, smaller, but 
in many ways more vital, questions may have their 
answer in surface and near surface energy budget 
and related microclimatic measurements. On an 
intermediate scale, there are questions such as 
do wet or dry years and/or months effect the 
release of Dissolved Organic Carbon into the soil 
water and fluvial systems? On the century to 
millenial scale there will be some time scale 
over which the vegetation zonation with respect 
to altitude (Greenland et al, 1985b) will change. 
It would be very interesting to determine the 
approximate size of this scale. In the inter- 
action between the atmospheric and biospheric 
scientist, the latter should take the lead in 

the dialog by first specifying, with the aid of 
the former, what questions require answers. 


CONCLUSIONS 


In addition to the specification of the climate 
by means of the summary statistics and water 
budget data, the most important features of the 
climate of Niwot Ridge that arise from the fore- 
going analysis are as follows: 

1) The unusually low temperatures of the early 
1980s, designated here as the Niwot Minimum, were 
probably real but recovered to the long term 
average in the late 1980s. 

2) The downward trend in temperature noted by 
Barry (1973) at Cl is now apparent at D1 and has 
continued. 

3) Potentially useful Autoregressive Moving 
Average forecasting models were produced to pro- 
ject annual temperature and precipitation values 
to the end of the century and to help diagnose 
the existing record. The diagnosive function did 
not turn out to be very useful. Whether the pro- 
jected values are useful or not will only become 
apparent in future vears. 

4) Spectral analysis indicated statistically 
significant periodicities in annual temperatures 
at 12.0, 9.00, and 7.20 years and in annual pre- 
cipitation values at 7.19 and 9.00 years. 


VYalues at Silver Lake in 1921. 
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5) Precipitation at Dl and Cl has shown an up- 


to 1987. This conclusion 
recalling that the 1951-64 D1 


ward trend from 1951 
must be qualified by 
data are estimated. 
6) Annual precipitation values at Silver Lake 
from 1913 to 1987 show a slight decreasing trend 
but. when the extreme value of the 1921 year is 
removed, this is actually a slight increasing 
trend. 

7) 1921 has been identified as an extremely 
unusual year at Silver Lake with a precipitation 
value of 1735 mm (68.30 ins) which is over twice 
the average of 755 mm (29.75 ins) for the 1914- 
1987 period. 


None of these changes have had obvious effects 
on the flora of the locality. Assuming there have 
been some changes, this would indicate a number of 
possibilities. First, the climatic variables 
used here are not sensitive enough to be related 
to vegetation changes. Second, the vegetation 
is not being monitored in enough detail or by 
means of the most sensitive variables that would 
indicate change. Overall, however, the lack of 
dramatic change in the vegetation testifies to 
its resilience within the range of values of the 
climatic variables reported here. 
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CLIMATE VARIABILITY IN THE SHORTGRASS STEPPE’ 


Timothy G.F. Kittel? 





Abstract.— Climatic variability at the Central Plains Experimen- 
tal Range LTER in northeastern Colorado was evaluated from 
annual, decadal, and century perspectives. Spectral analysis 
separated climatic variation into 2-7 y and >10 y scales. The 
record is dominated by cold-wet and warm-dry years. El Nifio 
years had on average higher annual precipitation than years with 
cold water anomalies in the equatorial Pacific. At the decadal 
level, variations in temperature and precipitation parallel those 
for the Central U.S. and broader regions. The temperature 
record has a significant 75-y trend of +2°C. Such variations in 
climate have impacts on the structure and function of the short- 


grass ecosystem. 


Keywords: Colorado, Central Plains Experimental Range, Long- 
Term Ecological Research, scale, El Nifio/Southern Oscillation, 
ecological response to climate variability, Vegetation Index. 


INTRODUCTION 


Ecosystems are responsive to climate variability at 
a range of temporal scales. Temperature, precipitation, 
and other climate variables (e.g., humidity, cloud cover) 
vary year-to-year and decade-to-decade in ways that are 
likely to influence ecosystem composition and function. 


I present an analysis of variability in temperature 
and precipitation at a range of temporal scales (annual, 
decadal, century) for the shortgrass steppe region in 
the vicinity of the Central Plains Experimental Range 
(CPER) in northeastern Colorado. Questions I address 
are whether there are characteristic scales of variation 
in temperature and precipitation, whether these scales 
are the same for the two variables (e.g., are dry 
periods characteristically warm or cold?), and whether 
there are long-term trends in the record. These are 
important questions because they seek to characterize 
periodic and unidirectional changes in the external 
forcing on ecosystems. Such understanding provides a 
basis for the study of between-year and long-term 
changes in ecosystem dynamics. 


‘Presented at the Long-Term Ecological Research 
(LTER) Workshop on Climate Variability and Ecosys- 
tem Response, Nederland, CO, August 21-23, 1988. 


"Research Associate, Natural Resource Ecology 
Laboratory and Cooperative Institute for Research in 
the Atmosphere, Colorado State University, Ft. Collins, 
CO 80523 


In the context of interannual variability, I examine 
possible long-distance links (teleconnections) between 
the climate of the steppe region and a source of global 
climate variability, El Nifio and the Southern Oscilla- 
tion (ENSO). ENSO consists of quasi-periodic (general- 
ly 3-7 y) oscillations in ocean-climate dynamics of the 
tropical Pacific Basin (Philander and Rasmusson 1985). 
As part of these oscillations, El Nifio is a strong 
warming of surface ocean waters in the central and 
eastern equatorial Pacific. ENSO teleconnections play 
a role in the climate dynamics of the extratropics 
(Yarnal 1985); impacts on ecological dynamics have 
been suggested for sites in the United States (e.g., 
Strub et al. 1985, Gosz 1988). 


Finally, I discuss observational and modeling 
studies on the impacts of climate variability on the 
steppe ecosystem. 


METHODS 


Site Description 


The CPER (40°49’N, 104°46’W, 1650 m above sea 
level) is operated by the USDA Agricultural Research 
Service (ARS) and is an NSF Long-Term Ecological 
Research (LTER) site. The vegetation at the CPER is 
shortgrass steppe dominated by perennial grasses 
Bouteloua gracilis and Buchloé dactyloides, with 
important succulent (Opuntia polyacantha) and sub- 
shrub components (Moir and Trlica 1976). Much of the 
surrounding region is grazed native grassland and 
cropland. 
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The climate of the shortgrass steppe is characteris- 
tic of mid-latitude, continental semi-arid regions; it is 
classified as BSk (cold, steppe climate) in the Képpen- 
Geiger system (Miller 1982). The mean climate and 
water budget is described by Parton and Greenland 
(1987) based on the period 1951-1980. Monthly mean 
temperature peaks in July (21.6°C) and is lowest in 
January (—3.1°). Mean annual precipitation is 309 mm, 
with maximum monthly precipitation usually occurring 
in May and June. Winter precipitation is associated 
with mid-latitude cyclonic storms and summer rainfall 
with cumulus convective storms that often initiate 
along the Rocky Mountain Front Range to the west 
and intensify as they progress eastward. Mean month- 
ly precipitation does not exceed monthly actual evapo- 
transpiration year-round (Parton and Greenland 1987). 


Climate Data 
The CPER 


I used temperature and precipitation data from two 
climate stations on the CPER, one operated by ARS 
and the other by LTER. The ARS station is at an 
elevation of 1645 m above sea level (asl); the LTER 
station is at 1643 m asl. The stations are ~5 km 
apart. ARS data are for the period 1940-1973 (tem- 
perature data start in 1948, precipitation data are 
continuous from April 1947). LTER data used here are 
from July 1970 to May 1988. The records of these 
stations were concatenated, with the LTER data replac- 
ing ARS data in October 1970; this was considered 
appropriate because the stations are in close proximity 
and have strongly correlated traces for the period of 
overlap (monthly temperature r’=0.99, monthly precipi- 
tation r’=0.74; Parton and Greenland 1987). Monthly 
mean temperatures were calculated as an average of 
monthly mean minimum and maximum temperatures. 


Grover 


To analyze long-term behavior in regional climate, I 
used data from a National Weather Service cooperative 
station approximately 16 km west of Grover. This 
station, officially named Groverl0W, has a 58-y record 
(1911-1969). It is roughly 25 km east-northeast from 
and 100 m asl below the LTER station. The station 
was moved several times during its history. For most 
of its operation (1929-1962), it was at 40°52'N, 
104°25’ W, and 1547 m asl. A significant move 
occurred in 1962 and is reflected in the temperature 
record (Fig. 1a). 


Neither temperature nor precipitation monthly 
records for Grover are without gaps. To create a 
continuous record, gaps of usually one month (occasion- 
ally up to three months) were filled with corresponding 
long-term monthly means for 13 years in the tempera- 
ture record and two years in the precipitation record. 
An advantage of this substitution is to retain years 
with nearly complete monthly data, although it reduces 
interannual variability and trends and blurs spectral 
signals. 


Analysis 


For analysis of interannual variability, blocking the 
data into 12-month periods is often arbitrary (e.g., 
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calendar years). In this study, I used water year 
(October to following September) to associate fall and 
winter precipitation with the following growing season. 
Consequently, this annual period includes the approxi- 
mate integrating period over which an important 
climate-ecosystem interaction, soil water recharge, 
operates. Water years are identified by the year the 
period ends: e.g., water year 1982 is October 1981- 
September 1982. Five-year running means of these 
data were constructed to investigate decadal and 
longer-term behavior. 


Spectral analysis was used to identify characteristic 
scales of variation in the yearly data. This analysis 
characterizes periodic behavior in a time series and 
identifies frequencies at which there is high covariance 
between two time series (i.e., high coherence). This 
analysis was accomplished using time series analysis 
software (BMDP/PIT; Thrall and Engelman 1985). 
Statistical software (SPSS/PC+; SPSS and Noruéis 
1986) was used to test for linear trends in the data. 
The time series analyses were run on detrended data 
and used a wide band width. A wider band width 
strengthens the statistical stability of the analyses, 
although it reduces the spectral resolution. 


Spatial Variability 


Spatial variability may give rise to problems in 
interpreting these station data in a regional and 
ecological context. Point measurements of precipitation 
may not represent the 6280 ha CPER because winter- 
time winds redistribute snow across the landscape and 
summer convective precipitation events are often patchy 
in their occurrence. For example, growing season 
precipitation varied across the site by as much as 96 
mm in 1988, an amount greater than the standard 
deviation of annual precipitation (Hazlett 1990). These 
effects can be ecologically important, influencing local 
soil water recharge (Burke et al. 1989) and contribut- 
ing to spatial variability in biomass accumulation 
(Lauenroth and Milchunas, unpublished data) and 
community structure. 


Southern Oscillation Index 


To evaluate possible teleconnections with El Nifio 
episodes and the Southern Oscillation (ENSO), I 
compared CPER data to a Southern Oscillation index 
(SSTso) based on central and eastern equatorial Pacific 
sea surface temperature anomalies developed by Wright 
(1984). The sign of this index is opposite to that of SO 
indices based on sea level pressure: SST xo is positive 
during El Nifio episodes and negative during cold water 
("anti-El Nifio" or "La Nifia") episodes. Wright’s (1984) 
record for this index is for 1950-1983. Water year 
averages of SSTs,. were included in correlation and 
spectral analyses. 


RESULTS AND DISCUSSION 
Water Year Time Series 
Station Comparison 
During the period of record overlap (1948-1969), 


there is a high correlation between CPER and Grover 
water year precipitation totals (r=+0.67, p<0.001; 








Table 1, Fig. 1b). However, precipitation during the 
overlap period differs between the two sites (p<0.05, t- 
test, Table 1). Greater precipitation at Grover likely 
reflects a regional gradient in precipitation (Cowie and 
McKee 1986): northeast of CPER, Grover is less in the 
rain shadow of the Front Range. 


There is also a difference in temperature means 
between sites during the overlap period (p<0.01, t-test; 
Table 1, Fig. la). This is primarily due to a diver- 
gence in temperature series after 1962 that reflects the 
change in Grover’s location. For the period of overlap 
up to this station change, annual temperature means 
are not significantly different (p>0.5, t-test) and are 
highly correlated (r=+0.71, p<0.005; Table 1, Fig. 1a). 
The difference in temperature means calculated for 
each station’s record (AT opge_gr=t+0.5°C, Table 1) may be 
due to climatic trends and is discussed later. 


Power Spectra 


Both high and low frequency variability is evident 
in the temperature and precipitation time series (Fig. 
1). The temperature power spectrum for Grover (Fig. 
2a) shows peaks in variance for oscillations at periods 
around 214-4 y and greater than ~20 y. The spectrum 
for precipitation (Fig. 2b) is broader with less well 
defined peaks at ~2 and 2%-5 y, as well as at >10 y. 
This reveals a natural break in temporal dynamics 
between annual and decadal scales. Dye (1983) found 
comparable breaks in the precipitation spectra for 
stations throughout Colorado. 
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Figure 1.— (a) Water year mean temperature, °C, and 
(b) precipitation, mm/y, for CPER (- - - -) and 
Groverl0W ( ). Horizontal lines are record means 
for each station. 
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Table 1— Summary of CPER and Groverl0W (GR) 
climate data for record and for period of overlap (1949- 
69 for T.,,, 1948-69 for ppt,,,). Numerical subscripts 
indicate other periods (e.g., 12-62 = 1912-62). Years 
are water years. Trends over station record (and also 
over 1912-62 for Grover) and corresponding significance 
levels (two-tailed t-test) are given. Station comparison 
over the full overlap period and over the overlap until 
1962 are shown based on correlation analysis (r) and 
paired two-tailed t-test for difference of means (t). 
Level of significance for trends, r, and t are: * = 
p<0.05 and ** = p<0.01. T,,,, = water year mean 
temperature (°C); ppt,,, = water year precipitation 
(mm/y). Elevations (m asl) are approximate because of 
station location changes. sd = standard deviation. n/a 
= not applicable. 


CPER GR A(CPER-GR) Ir t 
Elevation 1643 1547 +96 
Record: 
dg 
dates 1949-87 1912-69 
mean 9.0 8.5 +0.5 
sd 1.3 1.0 +0.3 
trend +2,4** +2.1** 
trend,s.¢0 n/a +1.7** 
PPtaa 
dates 1948-87 1913-69 
mean 317 347 —30 
sd 88 95 —7 
trend +52 —15 
trend. ¢2 n/a +1 
Overlap: 
Tee, 
mean 8.2 9.0 -—0.8 +0.14 3.2** 
sd 0.9 0.7 +0.2 
MEAN Ag. 69 8.6 8.6 0.0 +0.71** 0.5 
PPtinn 
mean 315 356 —41 +0.67** 2.4* 
sd 102 93 +9 
MEAN 4s 69 300 373 —-73  +0.77** 4,5** 





Coherence (h?) is analogous to the squared correla- 
tion between two variables within a given spectral 
band. A high coherence identifies a frequency at which 
variables strongly co-vary. There are peaks in the 
coherence between Grover temperature and precipita- 
tion series (Fig. 3a) in the short-period region (214-7'% 
y, h’>0.3) and in the long-period region (>20 y, h?>0.4). 
Analysis for CPER reveals a similar pattern with 
strong coherence at 3-3% y (h’=0.5) and at 20 y 
(h’?=0.4). The phase of these spectral correlations (Fig. 
3b) is such that temperature leads precipitation by ~% 
cycle across all frequencies; i.e., annual mean tempera- 
ture and annual precipitation fluctuations are out of 
phase. Significance of the coherence and phase is 
discussed below in the context of annual and decadal 
variability. 
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Figure 2. Power spectra for Grover water year 
(a) temperature and (b) precipitation. Band width = 
0.1167 y’. 


Interannual Variability 


Comparison of temperature and precipitation time 
series (Figs. 1a and b) suggests that peaks in annual 
temperature often coincide with low precipitation years 
and vice versa (Table 2). Indeed, annual temperature 
and precipitation are negatively correlated (r=—0.41, 
p<0.02). This is in accordance with the phase shift of 
¥% cycle found in the coherence results (Fig. 3b). The 
alternation of warm dry years with cold wet years 
implies control by oscillations in position and intensity 
of the mid-latitude polar jet and subtropical high 
pressure centers: e.g., cold wet years occur when the 
jet is farther south than average, warm dry years 


Table 2.— Contingency table of cold vs. warm and wet 
vs. dry years for Grover. Years are blocked based on 
water year mean temperature and precipitation greater 
vs. less than corresponding local (5-y running) mean to 
adjust for decadal trends. (7=5.175, df=1, p<0.025). 
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Figure 3.— Cross spectral (a) coherence, h’, and (b) 
phase for Grover water year temperature and precipita- 
tion. 


when farther north. While the predominance of cold- 
wet and warm-dry years is significant (p<0.025, 

Table 2), the alternative combinations of temperature 
and precipitation anomalies are present in the record 
(Table 2). Because controls over summer versus winter 
precipitation-generating mechanisms are relatively 
independent, control over variation in annual statistics 
is probably complex. Analysis of seasonal statistics 
may more readily reveal continental-scale controls over 
regional climate variability. 


Complexity in climate variability at the regional 
level is illustrated by decoupled variation in climate at 
CPER and Niwot Ridge, ~100 km to the southwest 
(Greenland 1990, this volume). The difference in 
behavior is likely because the climate of Niwot Ridge 
(3000-3750 m asl) is more strongly dominated by 
westerly flow across the continental divide. 


El Nifio Teleconnection? 


Spectral analysis revealed peaks in temperature and 
precipitation spectra and their coherence at periods of 
212-74 y (Figs. 1 and 2a). This range of periods is 
similar to that recognized for the Southern Oscillation 
(2-10 y, Rasmusson and Carpenter 1982) and so hints 
at a possible link to El Nifio episodes. 


There is a positive correlation between the annual 
Southern Oscillation SST Index (SST,,) and CPER 











annual precipitation (r=+0.30, p=0.11), such that 
precipitation tends to be higher during El] Nifo years. 
There is no significant correlation with annual tem- 
perature. However, at periods of 3-4 y there is 
spectral coherence between SST s, and CPER tempera- 
ture (h?=0.5), as well as between SST,, and CPER 
precipitation (h?=0.3). The ratio of annual precipitation 
during years with a strongly positive SST. (El Nifio 
years) to those with highly negative SST, is 1.13 
(Table 3), although only 5 of 9 El Nifio years between 
1951 and 1983 had greater than average precipitation. 
Sheaffer (personal communication) found that the 
corresponding El Nifio to “anti-E] Nifio" ratio for 
summer precipitation in northeastern Colorado is on 
the order of 1.33. By comparison, the Sevilleta LTER 
region (central New Mexico) is more strongly influenced 
by tropical Pacific climate dynamics: the El] Nifo to 
"anti-El Nifio" ratio is 1.59 for annual precipitation and 
2.45 for winter-spring precipitation (Molles and Dahm, 
in Gosz 1988). 


Observational and modeling studies of ENSO 
indicate that strongest teleconnections should be in 
winter, associated with the polar jet. For northern 
central North America inclusive of the shortgrass 
region, these teleconnections include increases in winter 
temperature and surface pressure (van Loon and 
Madden 1981, van Loon and Rogers 1981, Rasmusson 
and Wallace 1983). This is in accordance with greater 
annual temperatures observed at the CPER during El 
Nifio years (Table 3). Increases in winter surface 
pressure suggest that El Nino years should have drier 
winters (Yarnal 1985). Sheaffer (personal communica- 
tion) found this to be the case for NE Colorado: 
winter precipitation El Nifio to "anti-El Nifio" ratio is 
<0.75. However, depressed winter precipitation has 
only a small impact on total annual precipitation 
during El] Nifio years since most precipitation at the 
CPER comes during summer. 


Table 3.— Comparison of CPER water year mean 
temperature (T,,,,, °C) and precipitation (ppt,,,, mm/y) 
for years with high vs. low eastern and central 
equatorial Pacific sea surface temperatures (SST). 
High SST years correspond to E] Nifio episodes and 
were selected as those years with annual SST, 
exceeding the mean + % its standard deviation (s.d.). 
SST, falls below the mean — % s.d. in low years and 
within the region prescribed by the mean + % s.d. in 
‘Near Mean’ years. Analysis is for water years 1951- 
1983. n = number of years. n/a = not applicable. 


SST PPtann Toe n 
High 343 9.3 9 
Near Mean 313 8.7 14 
Low 303 8.7 10 
High:Low 1.13 n/a 

A(High-Low) +40 +0.7 

Meang.43 318 8.9 33 


While a clear ENSO teleconnection is weak for the 
CPER, difficulty in finding such a signal in the climate 
of the steppe is expected. Four factors contribute to 
the problem. (1) Teleconnections are partly a function 
of mid-latitude circulation patterns prior to ENSO 
episodes (Yarnal 1985). Consequently, the character of 
extratropical teleconnections shifts over decades with 
shifts in hemispheric circulation (Carleton 1987, 
Sheaffer and Reiter 1985). (2) Southern Oscillation 
indices can fail to catch the timing or intensity of an 
ENSO episode because of variations in the location of 
equatorial SST and/or atmospheric anomalies (Yarnal 
1985). This is because such indices are empirical and 
do not adequately reflect ocean-climate dynamics. (3) 
Large-scale studies show that ENSO-linked variation in 
climate for northern and central North America is less 
consistent in nature than for other parts of the 
continent (van Loon and Madden 1981, Blackmon et al. 
1979). (4) ENSO teleconnections account for only a 
small portion of extratropical climate variability 
(Wallace and Blackmon 1983). 


Decadal Variability and Long-Term Trends 
Annual Temperature 


A warming trend (AT=+3°) through the late 1930’s 
and subsequent cooling trend until 1950 (AT=—1.5°) are 
observed in 5-y running means of the Grover data (Fig. 
4a). A ~3.5° warming in the 1970’s is shown in the 
CPER data. These patterns agree with the general 
behavior of Northern Hemispheric temperatures 
reported by Hansen and Lebedeff (1987) and with that 
for eastern Colorado (Doesken et al. 1989) and the 
Midwest United States (Diaz and Karl 1988). The 
amplitude of the 1920’s to 1950 oscillation in the 
Grover record is magnified from that for the Northern 
Hemisphere (+0.4° warming, —0.15° cooling) and the 
Midwest (+1.5°, —1.0°), reflecting the steppe’s continen- 
tality. While a regional trend in the 1950’s and 1960’s 
is not clear due to the Grover station location change 
in 1962 (Fig. 4a), the trend in CPER data clearly 
follows that for other stations in eastern Colorado 
(Doesken et al. 1989) and the central United States 
(Diaz and Karl 1988): a slight warming in the early 
1950’s followed by cooling into the late 1960's. 


The long-term linear trend in the Grover tempera- 
ture record (until the 1962 station move, 1912-1962 
water years) is +1.7°C; that for the 1949-1987 CPER 
record is +2.4° (Table 1). The accumulative 75-y trend 
is 2.1°, not adjusted for station differences. The slopes 
of the CPER.and Grover trends are significantly 
different from zero (p<0.01, t-test), in spite of high 
year-to-year variability. This contrasts with Diaz and 
Karl’s (1988) conclusion that there has been no 
century-scale change in temperatures across the United 
States. Doesken et al. (1989) showed that trends in 
mean annual temperature in eastern Colorado are 
dominated by increases in winter temperatures. 


The result of a significant trend must be viewed 
cautiously before concluding that a regional change in 
climate has taken place. These results are only for 
two stations, each of which has a poor history with 
respect to location changes and data gaps. Other 
possible station changes such as in instruments, 
environs (e.g., buildings, vegetation), and time of day of 
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Figure 4.— As Figure 1, but for five-year running 
means. 


observation could potentially play a role in changing 
station means. 


Annual Precipitation 


The five-year running means of annual precipitation 
show significant decadal variation (Fig. 4b), such as the 
1930’s drought. Again there is considerable disagree- 
ment between Grover and CPER during the period of 
overlap. However, the CPER record agrees with the 
independent analysis of United States precipitation by 
Bradley et al. (1987) and shows the wide-spread 1950’s 
drought. 


Linear trends in Grover and CPER precipitation 
series are not significant (p>0.25, t-test; Table 1). This 
is in accordance with Bradley et al. (1987) who found 
no trend in annual precipitation for the United States 
since the 1850's. 


Covariance 


There is a general pattern of temperatures increas- 
ing from the 1910’s until the mid-1930’s, followed by a 
decrease into the late sixties and an increase since 
then (Fig. 4a). Superimposed on this apparent 50-y 
wave in temperature, recognized by Wigley and Raper 
(1987) for Northern Hemispheric temperatures, is a 
more rapid oscillation (15-20 y period) in precipitation 
(Fig. 4b). As a consequence, changes in the 5-y 
running means of temperature and precipitation follow 
the same pattern during some decades and opposite in 
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others. However, there is spectral coherence of 
temperature and precipitation at periods of 20 y or 
greater with a % cycle phase shift (h’>0.4, Fig. 3). 
This suggests that, as at the annual scale, multi- 
decadal variation in temperature and precipitation 
tends to be negatively correlated. This is supported by 
the pattern of droughts in the 1930’s and 1950’s being 
warmer than average and the wet 1910’s and 1960’s 
being colder. 


Some of the observed decadal or longer-term 
changes may be linked to continental-scale shifts in 
circulation patterns that are in response to hemispheric 
or global changes in atmosphere-ocean climate dynam- 
ics (Balling and Lawson 1982, Sheaffer and Reiter 
1985). Such shifts involve changes in the seasonal 
variation in intensity and position of the mid-latitude 
jet and the subtropical highs. Such changes signifi- 
cantly affect important abiotic controls over population, 
community, and ecosystem processes, such as the 
magnitude and timing of frosts and early growing 
season precipitation. 


ECOLOGICAL SIGNIFICANCE 


Identification of ecological responses to observed 
climatic variations is generally limited by a paucity of 
long-term biological data. However, some conclusions 
can be drawn based on available observations and 
modeling studies. 


Interannual Variability 


Net primary production of the shortgrass grasslands 
is strongly influenced by interannual fluctuations in 
precipitation (Lauenroth 1979). Lauenroth et al. (in 
preparation) compared field estimates of annual 
aboveground biomass production and annual precipita- 
tion at the CPER from 1941 to 1987. These data show 
that an extremely dry year reduces production for that 
and several years following. This suggests that dry 
episodes result in a reduction in the number of 
aboveground active meristems, constraining above- 
ground production in subsequent years, and cause 
belowground dieback, reducing the plants’ ability to 
acquire water also in ensuing years (Webb et al. 1978). 


Interannual variation in biomass development in the 
steppe is also influenced by the timing and amount of 
precipitation within the growing season. Seasonal 
dynamics of aboveground photosynthetically-active 
biomass can be inferred from satellite-based observa- 
tions, such as the normalized-difference vegetation 
index (NDVI) (Goward et al. 1985, Justice et al. 1985). 
NDVI is derived from data from the Advanced Very 
High Resolution Radiometer (AVHRR) on NOAA series 
polar-orbiting satellites. Comparison of seasonal 
variation in NDVI and precipitation (Fig. 5) suggests 
that in 1983 and 1984 high March-May precipitation 
resulted in a strong June peak in green biomass. In 
contrast, in 1982, high mid- through late growing | 
season precipitation corresponds to only a slight 
increase in late season biomass. These results suggest | 
that early growing season precipitation is more crucial 
for the development of aboveground biomass than late 
season rainfall. 











The frequency and timing of climatic events can 
also be critical for population and community level 
processes. For example, Lauenroth et al. (1987) 
showed that interannual variability in the occurrence of 
a combination of spring temperatures and precipitation 
events results in marked year-to-year variation in 
seedling establishment of blue grama (Bouteloua 
gracilis), with long-term consequences for community 
structure (Coffin and Lauenroth 1990). 


These results demonstrate that to assess the impact 
of interannual climate variability on shortgrass 
ecosystems, between-year differences in climate need to 
be evaluated at seasonal and event levels as well as 
the annual level. While the resolution of an analysis 
must be tailored to response times of ecological 
processes of interest, the sensitivity of a particular 
process may not be limited to climate variability at a 
single scale because of interactions with biological 
processes operating at finer and coarser scales (Allen 
and Starr 1982, O’Neill et al. 1988). Such interactions 
generate time lags in the response of ecosystems to 
abiotic forcing. 


Decadal and Long-Term Trends 


Decadal and long-term climatic trends are expected 
to impact ecosystem dynamics through climatic controls 


a VEGETATION INDEX 


NDVI 








on rates of net primary production (NPP) and decom- 
position. Using a grassland ecosystem model (CENTURY, 
Parton et al. 1987), Kittel et al. (1990) tested the 
sensitivity of a Colorado shortgrass ecosystem to 
directional climatic change. In the model experiment, 
annual temperature was increased by 4.9°C (roughly 
twice the 1949-1987 trend observed at CPER) and 
annual precipitation increased by 46 mm, based on a 
CO,-doubling climate change scenario (Hansen et al. 
1984). After 50 years, simulated aboveground NPP 
increased by 40 percent. This response was driven by 
increased precipitation and increased availability of 
nitrogen that was released as soil organic matter 
decreased. The net loss of soil organic matter was due 
primarily to increased decomposition driven by the 
increase in temperature. These results illustrate the 
potential importance of climatic shifts to the shortgrass 
steppe and the role of system dynamics in determining 
the response to such change (Schimel et al. 1990). 


Long-term climatic trends are also likely to influ- 
ence grassland community structure by changing 
species establishment rates and shifting interspecific 
competitive interactions. Such structural changes 
would influence the impact of climate change on 
ecosystem processes. 


b WEEKLY PRECIPITATION 
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Figure 5.— (a) Envelope of weekly-composited normalized difference vegetation index, NDVI, for a 50x50 
km region centered on the CPER in 1982, 1983, and 1984. (b) Weekly averages of CPER daily 
precipitation rate, mm/dy, for the same years. Dashed lines in (a) show a minimal threshold of 0.07 
below which NDVI values are assumed to indicate no green biomass. 
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CONCLUSIONS 


Climatic variability of the shortgrass steppe exhibits 
important characteristics that form a basis for evaluat- 
ing its impacts on the grassland’s ecological dynamics. 
These are: 


(1) There is a natural break in the temporal 
dynamics of the steppe climate between annual and 
decadal+ scales. At both scales, variation in tempera- 


ture and precipitation is generally 4% cycle out of phase. 


(2) At the annual scale, the phase relationship 
results in a predominance of cold-wet and warm-dry 
years, perhaps controlled by shifts in the position of 
the mid-latitude jet and subtropical highs. 


(3) There is a positive correlation between annual 
precipitation and ENSO episodes. However, there are 
limits to the empirical detection of an ENSO signal 
because teleconnections are dynamic and may be 
inconsistent in this region. 


(4) 20-50 y scale variations in the steppe tempera- 
tures match those in hemispheric and United States 
records and are of greater magnitude because of the 
region’s continentality. There is shorter period 
variability in precipitation at the decadal scale, such 
that in some decades temperature and precipitation 
change in the same direction and in other decades in 
an opposite fashion. 


(5) There is a statistically significant 75-y trend in 
temperature (~+2°C), but known and unknown 
problems in station histories warrant caution in 
attaching importance to this trend. No significant 
trend was found for annual precipitation. 


Ecological studies indicate that interannual and 
longer-term climate variability have significant conse- 
quences for ecosystem, community, and population 
dynamics in the shortgrass steppe. Climate-induced 
interannual variability in ecological processes can be 
complex, exhibiting significant year-to-year lags and 
sensitivity to climatic variation at within-season and 
event levels. 
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CLIMATE CHANGE AND ECOSYSTEM DYNAMICS 
AT THE VIRGINIA COAST RESERVE 
18,000 B.P. AND DURING THE LAST CENTURY! 


Bruce P. Hayden2 


Abstract. -- The forest cover of the Virginia Coast 
Reserve was studied using a Boreal FORET model 
parameterized using periglacial conditions derived 
independently of the pollen record of the site. 
Modern conditions of the landscape dynamics were 
studied using historical aerial photography and t he 
meteorological record of coastal storms. 
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The Virginia Coast Reserve (VCR) includes a ‘ ie Ny ease 
string of barrier islands on the coastal plain of ey 
the trailing margin of the North American Plate. 

It extends 100 km along the seaward margin of 

the Delmarva Peninsula. It encompasses 13 ay 
barrier islands (14,170 ha), broad intervening “ 
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inlets and extensive back barrier islands and 
shallow bays (Figure 1). Distinct communities fetes 
include sandy intertidal, open beach, grassland, * ‘ 
shrub thicket, pine forest, mud flat, salt marsh Cpe hy 

and estuarine lagoon. Terrestrial vegetation is d 
conspicuously shore parallel in general 

Structure, with sharp transitions or ecotones . 
betweeen communities. taf 
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Ecosystem dynamics of the VCR are closely tied 
to climate dynamics. Coastal storms are the 


dominant agents of landscape change in this 
coastal system. Storm waves and storm surges 


VIRGINIA BARRIER ISLANDS 
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the Eastern Shore of Virginia. 
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Figure 2. -- Cross-section of a barrier island 
showing the dominant morphological and 
vegetation zonation with the arrows indicating 
the extent of various marine processes a low 
dune case (top) and high dune case (bottom). 


this sandy coast and are responsible for the 
overall shore parallel zonation of the coastal 
landforms and communities on the barrier 
islands. The general form of the island zonation 
is shown in Figure 2. Woody species are 
restricted in their proximity to the sea by the 
presence of aeolian sea salts. Wind pruning of 
the vegetation by easterly sea breezes and 
northeasterly storm winds results in a direction 
oriented canopy morphology. The presence of 
shore parallel coastal dunes and wind roll 
vorticies contribute to the zonation of 
communities on the islands (Clarman 1977). 
The transition between the beach zone and the 
sandy grassland immediately inland is forced 
landward (eroded) by storms but aggressive 
growth by grasses encroach seaward during the 
growing season. Over the past four decades this 
ecotone over the entire VCR has been moving 
landward at about 5 m/yr (S.D. = 8 m/yr). 
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Had the VCR _ existed 18,000 years ago, its 
seaward margin would have been’ 100 
kilometers east of its current position and part 
of the boreal forest biome (Emery and others 
1967, Bonan and Hayden 1989). Peat deposits 
found at this distance offshore contain 
freshwater sphagnum peat and a _ decidedly 
boreal forest pollen assemblage (Emery and 
others 1967). With the post glacial sea level 
rise, the barrier islands "migrated" westward by 
the process known as roll-over. Beach sands 
from the shoreface are driven across the island 
during storm surges and deposited on the back 
side of the island. Thus erosion on the seaward 
side and accretion on the landward side of the 
island has the net effect of "moving" the island 
landward. 


At several places along the VCR tree stumps of 
forests once found on the back side of the 
barrier island now emerge from the beach face. 
In other locations lagoonside salt marsh peats 
outcrop on the beach. 18,000 years ago a boreal 
forest extended seaward 100 km. We have 
used a FORET model to simulate the structdure 
of this forest using climatic conditions at the 


a7 


end of the last ice age and have evaluated the 
output using fossil pollen assemblags from the 
area. Fossil pollen was not used to “tune” the 
FORET model as a totally independent test of the 
model was desired. By 5000 BP sea level rise 
had slowed and the islands were just a few 
kilometers seaward of their current position. 
Currently the rate of movement or trangression 
is about 5 m/yr. This rate is only exceeded on 
the sand starved barrier islands fronting the 
Mississippi Delta. If the charging of our 
atmosphere with CO2 gives rise to a global 
warming and an increased sea level rise, then 
the rate of transgression of the Virginia Barrier 
Islands will also increase. 


The model ISLAND is designed to connect 


climate, geomorphology, hydrology and 
vegetation cover in a _ stochastic-dynamic 
structure. The test site for ISLAND model 


development, parameterization and calibaration 
is Hog Island. Hog Island is the main research 
site for the Virginia Coast Reserve LTER. The 
southern end of this island is eroding rapidly 
and the northern end is accreting. At the 
northern end of the island, new landscapes are 
evolving while at the southern end old 
landscapes are being truncated. The history of 
shoreline erosion is shown in map form in 
Figure 3. 


In 1920 there was a town, Broadwater, Virginia, 
on the south end of Hog Island. There were 
about 100 buildings a coast guard station and a 
lighthouse. The former location of Broadwater 
is now hundreds of meters offshore. The rapid 
rate of landscape change on the Virginia Coast 
Reserve means that large changes can easily be 
monitored at the time scale of LTER research 
programs. The dominant driving forces for 
change on the Virginia Coast REserve are 
decidedly climatological. 


The landscape dynamics of the VCR is driven by 
astronomical and wind tides, winds, aeolian sea 
salts, storm waves and surges, fair weather 
swell, overwashed sands and saline water and 
the catchemnt of fresh rain water. Nutrients for 
the higher elevation fresh water elements of the 
landscape are by means of atmospheric 
deposition. 
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Figure 3. -- Shorelines of Hog Island from 1852 
to 1980. The site of the 19th century town of 
Broadwater at the southern end of the island 
was seaward of the current shoreline. 


THE VCR FORET MODEL: 18,000 B. P. 


The record of fresh water peats with arboreal 
pollen 80 km offshore dating to 18,000 B.P. 
offers an opportunity to test the climate driving 
functions of FORET type forest models against 
pollent data. The version of the FORET model 
we used in our studies was that of Bonan (1988) 
which is adapted from Shugart (1986). the 
basic components of Bonan's model are 
schematically illustrated in Figure 4. Bonan's 
gap model was designed to simulate the modern 
northern boreal forest and has been tested at 
four North American locations under a variety 
of local site conditions. The model was found to 
correctly simulate existing forest composition 
and dynamics. 
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Figure 4. -- Schematic of a stand element used 


in Bonan's FORET model of the boreal forest. 


30 
In order to apply the Bonan model to the | 


Virginia coast 18,000 years ago, climate and site 
conditions for the VCR _ were’ estimated 
indepnedent of the pollen record which was 
reserved for model validation. Climate and site 
parameters needed to run the VCR FORET model 
are listed in Table 1. Periglacial climatic 
conditions at the southern margin of the 
Laurentide ice sheet 18,000 years ago were 
taken from Moran (1972) and adjusted to take 
into account the distance between the edge of 
the ice cap and the VCR. 





———- MODERN 
--~- 18,000 YR BP. 


AIR TEMPERATURE (°C) 








Model Temperatures 
Mean monthly air temperatures from the 


modern record and those for 18,000 B.P. are 
shown in Figure 5. Monthly standard deviations 





of temperature for the modern period were Figure 5. -- Monthly mean temperatures for 
used for the period 18,000 B.P. Summer time the modern and 18,000 B.P. conditions at the 
temperatures were about 10°C cooler than Virginia Coast Reserve. 

modern and winter temperatures were about 

like modern. This latter condition arises 

because of the adiabatic warming of air moving 

southward and down over the ice sheet and Model Precipitation 

thus the preculding of extremely cold Precipitation estimates for the end of the last 
temperatures. One significant effect of this ice age were modern values except for the 
moderate glacial winter thermal climate is the hurricane season. Wendland (1977) reported a 
absence of permafrost on the VCR. very low frequency of Atlantic hurricnaces at 


the end of the last ice age. Modern summer/fall 
rainfall averages for the years with no 
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TABLE I 
AGEMAX Maximum age of species (years) 
DBHMAX Maximum diameter at breast height (cm) 
HTMAX Maximum height (m) 
G Growth parameter 
LITE Shade tolerance classification (1, tolerant; 2, intermediate; 3, intolerant) 
SMOIST The maximum percentage of the growing season that the species can tolerate soil 
moisture below the wilting point 

Sprouting 

N The tendency for stump sprouting 

DBHMIN The minimum diameter at breast height for sprouting (em) 

DBHMAX The maximum diameter at breast height for sprouting (cm) 
KTOL Fire tolerance (1, tolerant; 2, intermediate; 3, intolerant) 
NUTR Nutrient stress tolerance class (1, tolerant; 2, intermediate; 3, intolerant) 
IPFR Ability to grow on permafrost (1, good; 2, poor) 
IMO Ability to reproduce on moss~organic layer (1, tolerant; 2, intermediate; 3, intolerant) 
IBW Vulnerability to spruce budworm outbreaks (1, high; 2, low) 
ALC Light level at which reproduction is inhibited 
GDDMIN Minimum growing degree-days in the species’ range 
GDDMAX Maximum growing degree-days in the species’ range 
SWTCH Reproduction switches [SWTCH(I) is true if the species has serotinous cones. 


SWTCH{(2) is true if the species has copious, light, wind dispersed seeds. SWTCH(3) is 
true if the species can reproduce by layering]. 





hurricanes were used. In September, the month 
of highest hurricane frequency, more than 40% 
of the modern rainfall on the VCR is the result 
of passing hurricanes and tropical storms. The 
percentage is less in adjacent months. 
Aggregated over the year a 250 mm shortfall 
relative to modern rainfall totals results from 
the adjustment made. Summertime rainfall 
standard deviations were reduced to reflect the 
reduction in rainfall variance associated with a 
hurricane free climate. Monthly precipitation 
for the modern versus that of 18,000 B.P. are 
Shown in Figure 6 along with potential 
evapotranspiration (PET). 


Potential Evaporation (PET) 


Using calculated mean monthly radiation and 
temperature, elevation, and saturation vapor 
pressures. Mean monthly potential evapo- 
transpiration (PET) was estimated by the 
method of Jensen and Haise (1963). Figure 6 
contrasts the monthly mean water balance for 
the modern period and for 18,000 B.P. the PET 
deficit for 18,000 B.P. summer season is much 
lower than is typical for the modern conditons 
(note the scale difference in the two charts in 
Figure 6). 


Cloudiness_and_ Radiation 

Cloudiness estimates in tenths sky cover used in 
the model were arrived at by using modern 
values except in the summer. Summer values 
used were the average of the October through 
April values. The rational for this was that the 
termal contrst of winter and summer was little 
different across the seasons and the location of 
storm tracks in the middle latitudes was also 
seasonally stalbe. 
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Figure 6. -- Water balance for the Virginia 
Coast Reserve for the modern and 18,000 B.P. 
conditions. 


Monthly total, diffuse and direct solar radiation 
fwere calculated for the latitude of the VCR 
using the method of Liu and Jordan (1963) and 
adjusted for cloudiness. Modern Observed 
versus calcualted 18,000 B.P. mean monthly 
solar radiation for the VCR is given in Figure 7. 
Mean monthly radiation was not adjusted for 
the Milankovitch or orbital variations. 


MODEL RESULTS 

The results of the Bonan FORET model presented 
here are the average of 30 model runs. Biomass 
and % composition of the forest stand are 
calculated annually and the length of the model 
run is 500 years. Figure 8 shows biomass 
model output statistics for the average of 30 
plots (runs) using a fire cycle of 200 years. The 
model seems to achieve an_ equilibrium 
Structure after about 250 years. In this run the 
“boreal forst" simulated was a white pine 
dominatd forest with more typically boreal 
Species as co-equal subdominates in terms of 
biomass. More will be said about this white 
pine boreal forest later. Table II summarizes 
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Figure 7. -- Radiation at the Virginia Coast 
Reserve for the modern and 18,000 B.P. 
conditions. 


the model output statistics in % composition for 
the average forest stand simulated. Output 
statistics for fire cycles of 200 years and 50 
years are compared with pollen percentages 
from Emery and others (1967). It is clear from 
the table that at the general level the model 
simulation is rather similar to the observed 
pollen assemblage. Both model and pollen 
record s indicate a pine dominate forest with 
subdominates of spruce, birch and fir. The 
model output perhaps underestimates spruce 
percent composition. 
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Figure 8. -- Average biomass per species for 30 
model runs of 500 years each with a fire cycle 
of 200 years for climate conditions at 18,000 
BP. 








Fire cycie 
Mutneat rich Nutrient poor Dry site 

160 yr 200 yr 166 yr 260 yr 58 yr 100 yr 

Balsam fir 11.6% 10.8% 19.7% 17.9% 16.5% 12.3% 
Birch 18g 14.1 16.7 5.4 23.7 26.6 
Yellow birch 19.8 145 16.6 3.4 23.6 20.4 
White birch 6.6 66 61 6.8 67 6.2 
Red spruce #13 26.3 10.2 18.5 52 8.2 
Red cak ag 3 07 6.5 Rt ag 
Pine 56.5 544 58.6 56.8 57.5 579 
Fack pine 20.4 8.3 17.6 7S 35.3 26.5 
Red pine 6.2 O. 0.4 6.2 14 8.6 
White pine 35.¢ 45.8 48.6 48.4 20.8 8 





Previous investigatiors of Delmarva Peninsula 
pollen record have interpreted the record as 
evidence of a "classical" boreal forst and 
assumed that climatic conditions for the region 
were typical of conditions in the presenSt day 
boreal forest. It should be noted that we did 
not derive our climate data inut from inference 
from pollen records. if white pine. is not an 
available species in the model a forst of equal % 
compositon of pine, spruce, birch and fir results. 
This outcome is not consistent with the pollen 
record. 


We used a typical lapse rate for the atmosphere 
and adjusted our 18,000 B.P. VCR thermal 
climate to 1000 meters altitude so that we could 
simulate the forest of the Blue Ridge Moutains 
150 km west. Other climate data (modern 
values) were taken from the climate station Big 
Meadows. The simulation gave a typical spruce 
dominated boreal forest with both black and 
red spruce as the dominant’ spruces. 
Comparative pollen records for the crest of the 
Blue Ridge are not available. our model runs 
indicate that the boral forest of the mid-Atlantic 
region, from the Blue Ridge Mountians to the 
now-continental shelf area had a clinal shift 
from a spruce dominated boreal forest to a 
white pine dominated boreal forest. While pine 
pollen in published records do not go below the 
genera level, Harrison and others (1965) 
indicats that the pine pollen present at the end 
of the last ice age was large-grained and 
therefore probably white pine. In addition, 
white pine pollen is indicatd at Hack Pond, 
Virginia (Craig 1969). These records indicate 
perhaps a more general conclusion is in order: 
temperate latitude boreal forests in glacial 
times are white pine and spruce forests with 
spruce dominant in the higher altitudes and 
close to the glacial margin and white pine 
dominant at lower altitudes. This conclusion 
should stand as a working hypothesis awaiting 
further work. Two questions are important: at 
what elevation did the shift from a white pine 
to a spruce dominated forest occur; and at what 
latitude along the east coast did the white pine 
boreal forest replace the spruce dominated 
boreal forest. In subsequent runs of Bonan's 
Boreal FORET model for the VCR and using 3°C 
warmer growing season temperatures 
consistent with GCM model runs for the end of 
the last ice age, white pine was replaced by Jack 
Pine as the pine of the VCR boreal forest. 
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ISLAND 


ISLAND is a_ stochastically forced two 
dimensional dynamic model of a barrier island. 
The prototype version of the model was written 
at the University of Virginia (Rastetter 1989) 
and the ongoing island vegetion cover LTER 
monitoring is in part designed to provide critical 
parameterization for the model. In the 
prototype version of the model the following 
dynamics are included: shoreline’ erosion, 
movement of the beach/grassland ecotone, 
width of the active sand zone, coastal storm 
frequency, sea level change, dune construciton 
and erosion, island hydrology and_ island 
vegetation cover changes. 


The dynamics of the topographic dynamics of 
the barrier island are driven by aeolina and 
hydraulic sand transports. Of the two, hydraulic 
sand transport is by far the most important and 
this is driven by storm waves and storm surges. 
The maximum elevation that sand can _ be 
hydraulically delvated is around 4 meters 
above mean sea level. A small fraction of this 
sand can be elevated still higher by aeolian 
processees to form dunes. On the Virignia 
Barrier Islands dunes are rarely more than 2 
meters higher than the base of hydraulically 
transported sands. 


Horizontal changes in island topography are for 
the most part on the seaward side of the island. 
Storm waves and surges result in an offshore 
and across the island transport of snad. At 
present the net change in the isalnds is the 
transport of sand from the beach face across the 
islands. SAnd deposition on the subaerial 
portion of the island results in the general 
elevation of the island. This general elefvation 
of the island has to keep pace with the current 
rate of sea level rise (about 2 mm/yr). Sea 
level rise is an importna variable inour model 
ISLAND. 


Erosion 


Erosion of the seaward margin of the Virginia 
BArrier Island today averages nearly 5 m/year 
(Dolan and others 1979). We have records of 
this rate of change each 50 meters along the 
100 km VCR coast as well as the standard 
deviation of this rate of change (Figure 9) The 
model we are deveoping is a transect model and 
we plan to rund the model at each 5O meter 
location along the coast in order to arive at a 
three dimensional landscape model. 
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Figure 9. -- Historical erosion of shoreline of 


Hog Island over the last 5 decades. 


Topographic, hydrologic and vegetation transect 
data is being collected to parameterize the 
three-dimensional version of ISLAND. 


The rate of shoreline erosion is highly 
correlated with the frequency of storms 
offshore and our record of storminess off the 
VCR dates from 1885 (Figure 10). The 
correlation between shoreline eroision rate and 
storm frequency has been tabulated for each 
trasect for the period 1938-1989. While 
storminess off most of the U.S. Atlantic coast 
shows similar variation, shoreface dynamcis is 
most highly correlated with storms 100 to 200 
km offshore. Storm frequency recoreds for this 
location will be used to dynamically drive 
ISLAND landforms, hydrology and ecosystems. 


During storms, waves and _ surge penetrate 
inland as an overwash of oceanic wates and 
beach sands. The zone subjected to this 
overwash is scoured (cut) and then new sands 
are deposited (fill). In the growing season 
following the storm vegetation encroaches 
seaward and revegetates the overwash scar. An 
equilibrium between climate forces moving the 
active beach zone inland and the seaward 
encroachment of grasses seaward results. 
Because the shoreline is not constant in position, 
the boundary between the active sand zone and 
the edge of encroaching grasses also changes. 
The width of the active beach zone, the distance 
between the shoreline and the encroaching 
grassland, varies along the caost. Along the 
southern end of Hog Island the active beach 
zone is around 220 m wide. Along the central 
section of the island the active sand zone is only 
120 m wide. At the northern end there is a 
rapid shoreline accretion and new land is being 
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Figure 10. -- Frontal storm frequency off the 


Virginia Coast for the period 1885-1987. 


deposited. here the active beach zone is around 
520 m wide. 


ISLAND Dunes 


A dune erosion model has been constructed 
(DeKimpe 1987) which manages the demise of 
dunes as a function of the progress of the 
shoreline erosion and a model for the 
development of incipient dunes on_ the 
backshore is under construction. Woody shrub 
development on the ialnd is very much 
dependent on the protection against salt spray 
afforded by the presence of dunes. Studies on 
Assateague Island (Clarman 1977) indicate that 
the radial rate of growth and the elongation of 
shrub canopies are dependent on the direction 
of the sea breeze, onshore storm winds, and the 
presece of dunes. 


ISLAND _ Hydrology 


Barrier islands are typified by a _ shallow 
freshwater aquifer atop saline ground water 
(Bolyard, T. and others1979). The elevation of 
the water table surface of the fresh water 
aquifer is a function of the elevation of the sand 
surface, island width and precipitation input. 
Water table elevation in the model is driven by 
sand elevation and thus indrectly by sea level 
variations and by the frequency of coastal 
storms. Barrier Island vegetation is in 
generallly sensitive to ground water salinity 
(Schneider 1984). Vegetation along each 
transect is in part be governed by local model 
salinity. 


ISLAND Results 


Studies of barrier island vegetation dynamics 
(Clarman 1977, Elmer 1975, Schneider 1984, 
and Schroeder 1977) in addition to new 
vegetation transects are used to parameterize 
ISLAND. Detailed vegetation surveys in the 
early 1970s and an on going resurvey effort 
using photos and ground surveys are also used 
in the parameterization. Figure 11 shows 
schematically output from ISLAND for an 
eroding coastal unit. Island profiles, water table 
and vegetation cover are shown. Figure 11 
shows model output after 30 years following a 
simulated extinction of the grasses. The model 
is sensitive to grasses in that aeolian sand 
transport and deposition is involved in dune 
development which in turn alters the island's 
water table and thus the success of salinity 
sensitive vegetation. 


TOPOGRAPHY & 
SEA LEVEL 


GROUND WATER SALINITY 


1000 
ABUNDANCE = 500 


0 


BB Penenian KG annua [QQ mars crass 





Figure 11. -- ISLAND output for an 
experimental extinction of the grasses on the 
island. The results shown are abundances for 
perennial and annual terrestrial plants and for 
marsh grasses. 


CONCLUSIONS 


Research at the Virginia Coast Reserve LTER has 
three foci: investigation of environmental 
processes, monitoring landform and ecosystem 
change, and modeling ecosystem dynamics. 
Climate and climate dynamicas are importnat 
part of each of these efforts. Sea level change 
and storm climate change are important in 
landscape evolution and thuis ecosystem 
evolution while chagnes in weather elements 
(temperature, rainfall, cloudiness, etc.) over the 
time period of decades to centuiries have direct 
effects on the vegetation cover of the VCR. 
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Because the chnages on this sandy coast have 
been and continue to be very rapid, the VCR 
LTER is in an ideal position ot make substantial 
progress in the study of long term ecological 
processes within the contex of LTER funding 
cycles. 
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OVERVIEW OF CLIMATE VARIABILITY AND ECOSYSTEM RESPONSE? 


David Greenland and Lloyd W. Swift, Jr.? 


Abstract.-—— Unusual ecosystem responses are frequently 
driven by meteorological events. The frequency and 
Magnitude of these events and responses can be 
characterized through Long-Term Ecological Research. The 
LTER Climate Committee identifies four issues to be 
considered in future investigations: (1) the need to 
clarify terms and definitions used in discussing climate 
variability, (2) the importance of recognizing the various 
time and space scales of climate variability and ecosystem 
response, (3) the need to expand data beyond dependence on 
traditional summaries of temperature and precipitation, 
and (4) the value of insights gained from examining 
similarities and dissimilarities among climate episodes 
and ecosystem responses across LTER sites. 


Keywords: LTER, scale, climate change, air mass. 


INTRODUCTION 


An important contribution of long-term 
ecological research is the ability to place 
unusual ecological events in perspective. In a 
study of 380 ecology papers, Weatherhead (1986) 
concludes that "the danger of short-term studies 
may be that they experience too many unusual 
events. The reason for this unexpected conclusion 
may be that we tend to overestimate the importance 
of some unusual events when we lack the 
perspective provided by a longer study." He also 
notes that abiotic atmospheric factors, 
particularly precipitation and/or temperature, 
cause the great majority of unusual ecosystem 
events. In their report on long-term research, 
Strayer and others (1986) find that long-term 
studies are necessary to explore four major 
classes of ecological phenomena. They identify 
these phenomena as (1) slow processes, (2) rare 
events, (3) subtle changes in systems, and (4) 
complex processes requiring long-term multivariate 
studies to detect change. Note that the static or 
no change situation was not listed as an 


1 Developed from a discussion of the LTER 
Climate Committee held at the Workshop on Climate 
Variability and Ecosystem Response, August 21-24, 
1986 at Niwot Ridge-Green Lakes LTER Site, CO. 


2 Professor of Geography, University of 
Colorado, Boulder, CO, and Research 
Meteorologist, Coweeta Hydrologic Laboratory, 
Forest Service, U.S. Department of Agriculture, 
Otte, NC. 


ecological phenomenon. The first three of these 
classes of change may be closely correlated with 
climate data while climate data may be a 
significant variant in the fourth. Furthermore, 
according to Strayer and others, the measurement 
variables eventually selected in long-term 
research could be classified either as structural 
variables, such as species composition, or as 
functional variables such as primary productivity. 
Climate might be classified best as a set of 
functional variables, even though some of the 
functional relationships are not yet known. 


Following presentation of the papers in this 
volume, the LTER Climate Committee discussed the 
application of long-term studies to research on 
climate variability and ecosystem response. We 
focused on four main areas: (1) clarifying our 
terminology, (2) recognizing the importance of 
time and space scale in all aspects of such work, 
(3) developing and promoting climatic indices, 
other than standard expressions of temperature and 
precipitation, that may be useful in ecosystem 
studies, and (4) utilizing the similarities and 
dissimilarities between sites. 


TERMINOLOGY 


Even before the Workshop began, Committee 
members began debating the meaning of the term 
"climatic variability". The view was that the 
term, as used in the Workshop title, implies 
abnormality whereas climate variability (or any 
other kind of ecosystem variability) is the normal 
condition. Climate variability was described as 
consisting of a pattern of “episodes” and 
“events”. The time scale of the episode or event 
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is significant. The following points led the 
group into the concept of episodes versus events 
as it applies to Long-Term Ecological Research 
(LTER) . 


We defined weather as a real-time event, 
whereas climate is a synthesis or time integration 
of weather element values or weather systems. 
Climate has the component of expectation that a 
characteristic will occur. 


Clear terminology is necessary because weather 
events and climatic episodes have political 
ramifications when reported to the public by 
popular media. Reports must clearly state whether 
a particular heat wave or drought is a significant 
event or part of the variability that is an 
integral part of routine weather. Unfortunately, 
terms such as "climatic normal" and "normal 
climate” detract from the reality that variability 
is the normal characteristic of a climate system. 


Because of their long-term emphasis and 
extensive spread over diverse ecosystems, LTER 
sites are key elements in the national effort to 
detect changes in the climate system. Without 
data, the detection of change is strongly molded 
by human bias toward a time scale which 
corresponds to the human lifespan. Extreme 
weather events are ranked in severity against 
those events within the observer’s memory. The 
"Great Southeastern Drought of 1986-88" is most 
significant to those who cannot recall the drought 
of the mid-1920’s. A mature forest or natural 
grassland may be described as unchanging until 
someone obtains measurements of its 
characteristics over a period long enough to 
detect change. Investigators must be able to 
evaluate the state variables of an ecosystem 
before they can begin to relate change to climate 
variation. Consistently collected, long-term 
climatic data is a most valuable and necessary 
tool to categorize weather events and climate 
episodes. The LTER network provides opportunities 
to test for climate change and to relate it to 
ecosystem response without an anthropocentric time 
scale bias. This bias limits the ability to 
detect real changes between one episode and 
another. <A corcllary question might be, what 
other biases do we impose on the LTER ecosystems 
we study? For example, how do the sizes or 
sampling intensities of the LTER sites affect 
understanding of spatial scales? 


Both the cyclic nature of and rapid changes in 
climatic values have been noted in earlier 
chapters of this volume. Examples include data 
from Hubbard Brook (accumulated daily 
precipitation and temperature records), Northern 
Lakes (dates of lake freezing), and Niwot Ridge 
(recent annual temperature values). Step 
functions might be more common than smooth trends 
or cycles, even on an intra-annuel scale and for 
all elements. On a seasonel scale, phenological 
changes may force step~function shifts in 
microclimate values. Albedo, light transmission, 
and litter temperature, for example, change 
rapidly with leaf development or leaf fall. Other 
discontinuities in time series of climete 
variables are common at many LTER sites. These 
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discontinuities may be more important than regular 
trends because they "reset" the ecosystem. For 
example, every storm establishes a new state for 
soils, vegetation, and associated drainage 
systems. Rapid change is a prime characteristic 
of the interval (event) between episodes. 


For research on climatic cycles or step 
functions, members of the Committee cautioned 
against using methodology such as spectral 
analysis which looks for cyclic forms irrespective 
of the realities of the ecosystem. Such 
techniques should simply be used to search for an 
explanation of variance. Spectral analysis was 
applied to the Niwot-Green Lakes system and 
limited power was found to explain actual 
ecosystem operation. Other, more appropriate, 
methodologies may be available for investigating 
value discontinuities in ecosystem and abiotic 
variables. An example is the work of Walsh and 
Richman (1981) on the rotation of orthogonal 
principal components. By examining the sizes of 
anomaly fields, they were able to identify sister 
stations in both time and space and define 
discontinuities in the record. This technique 
could be very useful for extrapolating out from 
LTER sites. 


As a result of the above considerations, the 
Committee found that LTER scientists, and others 
working in the field of climatic variability, 
should be more specific than the term "climatic 
change" allows. For clarity, we should apply a 
distinction between "episodes" and "events". An 
"event" was defined as a single occurrence, such 
as a large rainstorm, often embedded in the 
functioning of the synoptic climatic scale. An 
"episode" was taken as a string of events and its 
duration is probably related to the time constant 
of the system. Some events and most episodes 
reset the time clock of the system. They result 
in a large change in the ecosystem at the time of 
eccurrence followed by a long tail of less obvious 
adjustments. 


The Committee recognized at least three 
versions of climate episode. First, climate 
episodes are defined by the data of the climatic 
time series bounded by their indications of 
changes of state. Second, the perceived climatic 
episode is often described by means of climatic 
data but is actually bounded by a time scale of 
human memory lasting between 40 and 80 years. 
Third, a climate episode may be best defined by 
responses of the components of the ecosystem. All 
are especially dependent on spatial scale and the 
latter is specifically apropos for Long-Term 
Ecological Research. 


SCALE 


In discussion, the Committee continually 
returned to questions about time and space scales 
in which episodes occur. Scale is an important 
consideration because it determines what kinds of 
questions can be asked about the operation of the 
ecosystem. Researchers must relate the scales on 
which climate systems operate to those scales on 
which the biotic parts of the ecosystems operate. 
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In making these time and space distinctions, we 
will be addressing the problems of complexity in 
the same sense as in the concepts of hierarchy 
theory. Those concepts were applied to ecology by 
such seminal works as Allen and Starr (1982) and 
O’Neill and others (1986). In organizing our 
ideas around specific time and space scales we 
will be dealing with an organized complexity 
instead of disorganized complexity. We will find 
that all parts of the system do not interact at 
the same time because of the very existence of 
different time and space scales. For instance, 
microbial respiration rates are more related to 
individual rain events than to gap/phase 
succession events in forests that have been 
subject to long-lasting droughts. This approach 
for simplifying organized complexity will enable 
us to structure our view of systems, but we may 
need to upgrade our key skills for sampling our 
systems. In all of these considerations, the 
functional ecosystem variables assume greater 
importance than the structural variables. 
Therefore, the climate veriables that relate to 
ecosystem function rather than to structure should 
be emphasized. 


Thus, we conclude that understanding climate 
variebility and ecosystem response demands thet we 
pey particular attention to space and time scales. 
We must bewere of arbitrarily imposed, human~— 
derived scales and concentrate on those sceles 
that emerge from the functioning of the ecosystem 
and climate systems. Research should specifically 
identify those functions and processes of the 
ecosystem that cannot keep up with potential rates 
of ebiotic change such as postulated global 
warming rates. 


INDICES FOR INTERSITE COMPARISON 


The LTER Climete Committee recognized a 
continuing need for consistency in obtaining and 
handling data across the Network. A set of time 
series analyses across all sites would be useful. 
Also useful would be new indices, not directly 
dependent on monthly and annual mean temperature 
and precipitation values, to extend the 
information base beyond our earlier work 
(Greenland 1987). One such index, the date of 
lake freezing, is ably demonstrated in this volume 
for the Northern Lakes site. However, this and 
another index related to storm surge data are 
specific for the LTER sites and ecosystems they 
represent. 


Other data sets exist that could provide useful 
and general climatic indices. For example, the 
Department of Environmental Sciences at the 
University of Virginia has records of cyclone 
frequencies since 1885 and 500 to 1000 mb 
thickness levels for all LTER sites. 


An index that seemed to have wide application 
for intersite comparisons emanated from air mass 
climatology. Climate at a place is dependent on 
exposure to a characteristic pattern of air mass 
types which integrate many climatic elements such 
as temperature, precipitation, and humidity. 
Wendland and Bryson (1981) refined the concept of 
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frequency of air mass climatology by using 
streamline analysis to map airstream regions. The 
regions are defined by the boundaries between 
airstreams from different global source areas. 
Wendland and Bryson traced the source of these air 
streams by mapping monthly surface level 
streamlines (j.e., lines of resultant winds along 
which air actually flows at any given moment). 
Every LTER site experiences periods during the 
year when there is a shift between being in the 
region of one airstream and being under the 
influence of air from another. An index for 
comparing LTER sites might be the number of months 
duration in different airstream regions. The time 
pattern of airstream regions could also explain 
the seasonal distribution of precipitation and 
strong site contrasts such as the extreme between 
Jornada and Andrews. Variation might increase 
with distance of a site from the source of the 
airstream. 


Wendland, who pioneered this work, has since 
examined air mass frequency data for all LTER 
sites (Table 1). These frequencies can be refined 
to ensure thet the boundaries for the air mass 
regions ere based on data representative for each 
LTER site. For example, the elevated Niwot site 
is not expected to be in the same air mass as the 
Plains site, both based here on Denver data. 

Table 1 indicates the duration of each air mass 
from various source regions and suggests the 
climate for the 1948-1963 period. In another time 
period, the air mass frequencies might change, 
especially at sites near the confluences of 
airstreams. Thus, this data form may provide 
evidence of moves from one climatic episode to 
another. 


The Committee thus recommends that sites, 
singly and as a network, investigate new and non- 
standard climatic indices to supplement the 
information obtained from standard climatic 
observations and summaries. Our goal is to define 
and refine relationships between climatic 
variation and ecosystem response. 


SIMILARITY AND DISSIMILARITY 


Outwardly, LTER sites appear so different that 
useful comparisons are either obvious or else 
impossible. A benefit of having LTER sites in 
very different biomes is that broad-scale 
comparisons, not often available to ecologists, 
can be made which should give valuable insights 
into ecosystem function and processes. This was 
demonstrated during the Workshop when similarities 
and dissimilarities between sites were examined. 


Many sites have not yet identified clear or 
obvious ecosystem responses to slow climate trends 
or even to events of mid-scale severity. But most 
sites have experienced major responses to a severe 
weather event. The Hubbard Brook ecosystem, for 
example, was not markedly disturbed by the 
droughts of the 1960s but still shows the effect 
of a single hurricane in 1938. This may be yet 
another exemple of our inability to perceive long- 
term changes. Tree blowdown has been a repeated 
catastrophic wind event at several LTER sites and, 





Table 1.--Mean number of months per year of 
domination by five air mass types for LTER sites, 
1948-1963. 





Site North North Ohio Arctic — High 

Pacific Atlantic Valley Plains 
and Gulf 

Andrews, OR 12 

Sevilleta, NM YW 1 

Konza Prairie, KS 5 5 2 

Hubbard Brook, NH 3 6 3 

Bonanza Creek, AK 2 7 5 

Jornada, NM 2 7 3 

Virginia Coast, VA I § 3 

Niwot Ridge, CO 3 5 2 2 

Central Plains, CO 3 5 2 2 

Okefenokee, GA 9 3 

North Inlet, SC 9 3 

Coweeta Lab, NC 1 i 4 

Luquillo For., PR 12 

Cedar Creek, HN 3 3 2 3 1 

Illinois Rivers, IL 3 4 2 j 2 

Northern Lakes, WI 3 4 1 2 2 

Kellogg Sta., HI 2 4 2 2 2 

Harvard For., MA ! 5 4 j 

Arctic Tundra, Ak 4 8 


since the Workshop, hurricane damage has 
significantly altered both the North Inlet and 
Luquillo sites. Many ecological responses are due 
to secondary effects of atmospheric events, such 
as flooding or landslides. For example, the 
redistribution of sediment by an intense rainstorm 
on the otherwise dry Jornada site has marked 
consequences on the biota either by burying them 
or by providing new microhabitats. 


Several sites reported possible time 
coincidence for deviations in climate variables. 
The years of climatic change suggested by shifts 
in freezing dates of Lake Mendota, WI, in 1880, 
1940, and possibly 1980, were noted as times of 
change at some other sites and also in general 
climate data. LTER sites may benefit from 
examining their own records for common break 
points in data sets. Data at most sites, as well 
demonstrated by the Central Plains Experimental 
Range, follow hemispheric, or at least regional, 


trends in temperature and precipitation. This 
augurs well for the extrapolation of results from 
the LTER network to larger areas. Yet, unique or 
isolated sites such as Niwot will not display the 
same spatial and temporal trends as adjacent 
dissimilar areas. 


At first the Kellogg Biological site was 
believed to be functionally different from other 
sites because of its emphasis on monoculture of 
agricultural crops and the attention given to 
short time-scale investigations emphasizing 
specific times of the year. These seasonal 
studies include winter impact on the life cycle of 
insect pests, spring weather affecting 
germination, and climatic influences on 
pollination and seed set. The lesson is that 
weather events are marked by phenological events, 
a phenomenon equally true at other LTER sites. 
The fact that the Kellogg ecosystem defines 
shorter time scales is another demonstration of 
the importance of recognizing all time scales as 
was discussed earlier. 


Discussion revealed that many LTER sites had 
considerably more data than simply monthly means 
and totals of temperature and precipitation 
values. In many cases, high-quality data for 
climate and ecosystem variables coexist. 
Opportunities were recognized for episodic studies 
on daily and other time scales in intersite LTER 
studies. 


In summary, several fertile areas for further 
research can capitalize on the similarities and 
dissimilarities of climate variability and 
ecosystem response across LTER sites. These 
include an investigation of (1) the importance of 
catastrophic events in relation to slower trends 
and cycles, (2) the time coincidence of certain 
major climatic breakpoints which appear to exist 
at several sites and the effects on ecosystems as 
they shift from one episode to another, (3) the 
relationship of climate to phenological studies 
across the LTER network. 


CONCLUSION 


Climatic variability and ecosystem response is 
clearly a topic having all the intricacies of a 
Gordian knot. Deliberations of the LTER Climate 
Committee have indicated some important starting 
points at which the knot may be unravelled. 

First, we must be very conscious of our 
terminology. Loose usage of terms may well hinder 
our conceptualization of reality. Second, we must 
put considerations of scale at the beginning of 
our investigations instead of making a priori 
assumptions about them. There is a tendency, of 
which we must be cautious, to impose human-— 
oriented concepts of scale on our real systems 
instead of letting the functions of the ecosystems 
themselves define our scale for us. Third, we 
have identified some exciting ways by which we can 
go beyond the use of simple temperature and 
precipitation values to relate to ecosystem 
functions or define breakpoints between climatic 
episodes. Finally, insights gained by comparing 
similarity and dissimilarity between the LTER 
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sites will improve understanding of on-site 
ecosystems as well as explain intersite variation. 


None of these ideas are new; but within the 
context of climate variability and ecosystem 
response at LTER sites, they take on a new 
Significance. The highly disparate nature of LTER 
sites allows the Committee to search for indices 
like air mass frequency that go beyond information 
restrained to local observations alone. This 
opportunity can lead to a broader search for new 
concepts and techniques in ecosystem science as a 
whole. The LTER Climate Committee Workshop 
generated ideas and concepts that should 
facilitate notable progress in understanding 
climate variability and ecosystem response. 
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